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ABSTRACT 

A  study  of  the  role  of  remote  sensing  for  geologic  reconnaissance  for  tunnel- 
site  selection  wr.s  commenced.  For  this  study,  remote  sensing  was  defined 
as  ultraviolet  to  thermal  infrared  multispectral  sensing,  S-  and  L-band 
synthetic  aperture  radar,  and  aerial  photography.  Data  from  these  sensors  were 
processed  and  evaluated  in  terms  of  their  complementary  use.  In  addition  to  gen¬ 
eral  geologic  interpretation  and  evaluation  of  remote  sensing,  specific  work  was 
performed  to  improve  the  state  of  the  art  and  thereby  contribute  to  better  geologic 
interpretation.  Our  areas  of  study  included:  (1)  image  ratioing  to  indicate  SiO,  and 
iron  oxide  content,  representing  the  final  step  toward  compositional  (lithologic  are 
remote  sensing  in  geology;  (2)  the  analysis  and  processing  of  remote  sensing  data 
from  a  glacial  drift  region  for  which  level  slicing,  ratioing,  and  multichannel  sta¬ 
tistical  pattern  recognition  of  multispectral  data  were  employed;  also  the  correla¬ 
tion  of  thermal  variation  studies,  aerial  photographs,  and  radar  imagery  to  obtain 
information  on  soils,  geology,  and  hydrology;  and  (3)  the  discussion  of  two-Jre- 
quency  multipolarized  radar  imagery  in  terms  of  its  capability  to  delineate  geo¬ 
logic  features.  This  report  can  ,j  used  tutorially  on  the  data  processing  and  basic 
instrumentation  of  conventional  remote  sensing.  Future  research  directions  are 
suggested,  and  the  extension  of  remote  sensing  to  include  airborne  passive  micro- 
wave  sensor  systems,  magnetometry,  gamma-ray  sensors,  gravimetry,  and  air¬ 
borne  electromagnetic  sounding  systems  is  discussed. 


V 


WILLOW  RUN  LABORATORIES 


CONTENTS 

Foreword . iii 

Abstract .  v 

List  of  Figures . lx 

List  of  Tables . x 

1 .  Summary .  1 

2.  Remote  Sensing  for  Tunnel-Site  Selection .  3 

2.1.  Introduction  3 

2.2.  Imaging  Remote  Sensors  3 

2.3.  The  Role  of  Remote  Sensing  in  Tunnel-Site  Selection  4 

2.4.  Related  Literature  Sources  7 

2.5.  Work  Performed  Under  This  Contract  8 

3.  Multispectral  Remote  Sensing . 10 

3.1.  Multispectral  Data  10 

3.2.  The  Multispectral  Concept  10 

3.3.  A  Multispectral  Scanning  System  11 

3.4.  Automatic  Processing  14 

3.4.1.  Single-Channel  Processing  16 

3.4.2.  Two-Channel  Processing  17 

3.4.3.  Multichannel  Processing  17 

3.5.  Data-Collectlon  Parameters  18 

3.6.  Applications  of  Automatic  Processing  of  Multispectral 

Data  19 

3.6.1.  Single-Channel  Applications  20 

3.6.2.  Two-Channel  Applications  26 

3.6.3.  Multichannel  Processing  Techniques  for  the 

Classification  of  Land  Uses  35 

4.  Compositional  Remote  Sensing  of  Rock  Outcrops  by  Image - 

Ratlolng  Techniques . 41 

4.1.  Imagery  Optimization  Analysis  41 

4.2.  Improvement  of  the  Thermal  Infrared  Ratio  Technique  52 

4.3.  Description  and  Improvement  of  the  Visible -Reflective 

Infrared  Ratio  Technique  64 

4.4.  Relevance  of  These  Methods  to  Tunnel-Site  Selection 

and  Other  Bureau  of  Mines  Problems  67 

5.  Applications  of  Sldelooklng  Radar  to  Geologic  Investigations . 68 

5.1.  Introduction  68 

5.2.  Simulation  of  Radar  Imagery  by  Means  of  Low  Sun-Angle 

Photography  70 

5.3.  Stereoscopic  SLAR  Imagery  and  Areal  Coverage  70 

5.4.  Radar  Return  Parameters  71 

5.5.  Use  of  Wide  Dynamic  Range  of  Radar  76 

5.6.  Optical  Processing  for  Enhancement  of  Geologic  Features  77 

6.  Future  Research  . 80 

6.1.  Recommendations  80 

6.1.1.  Field  Test  Cases  80 


vll 


- WILLOW  RUN  LABORATORIES  - 

6.1.2.  Continued  Development  of  Multispectral  Scanning 

and  Radar  Systemd  80 

6.1.3.  The  Addition  of  Other  Kinds  of  Sensors  For  a 

Complete  Remote  Sensing  System  82 

6.2.  Rationale  of  Approach  to  Future  Research  84 

References . 86 

Distribution  List . 90 


viii 


WILLOV.'  RUN  LABORATORIES 


FIGURES 

1.  The  Electromagnetic  Spectrum  and  Applicable  Sens'ng  Systems . 12 

2.  Multispectral  Data  Collection . 13 

3.  Comparison  of  Two  Video  Images  and  a  Level-Sliced  Image  for  the 

Discrimination  of  Standing  Water . 21 

4.  1.5-  to  1.8 -pm  Image  Mosaic  of  Oakville  Representative  Basin, 

Ontario,  for  Discrimination  of  Surface  Water . 22 

5.  High-Altitude  Photomosaic  of  Oakville  Representative  Basin, 

Ontario . 23 

6.  Comparison  of  Standard  Thermal  Video  Image  with  Calibrated  Level - 

Sliced  Thermal  Image;  Data  Collected  at  1200  hours . 25 

7.  Comparison  of  Standard  Thermal  Video  Image  with  Calibrated  Level- 

Sliced  Thermal  Image;  Data  Collected  at  0400  hours . 27 

8.  Thermal  Infrared  Images  of  Portion  A  of  No.  6  Highway  Between 

Hamilton  and  Guelph,  Ontario . 28 

9.  Thermal  Infrared  Images  of  Portion  B  of  No.  6  Highway  Between 

Hamilton  and  Guelph,  Ontario . 29 

10.  Comparison  of  Thermal  Images  Produced  from  Data  Collected  at 

0400  hours  and  at  1200  hours . 31 

11.  Comparison  of  Three  Video  Images  for  an  Area  in  Southeastern 

Pennsylvania . . 

12.  Comparison  of  Three  Ratioed  Images  for  an  Area  in  Southeastern 

Pennsylvania . 34 

13.  Aerial  Photographs  and  Video  Images  of  Oakville  Creek 

Basin  Area . 36 

14.  Recognition  Images  of  Six  Different  Land-Use  Classes  in  the 

Oakville  Creek  Basin  Area . 39 

15.  Ratio  Between  Radiances  in  Channel  1  and  Channel  2  Versus  the 

Percentage  of  Weight  of  SiOg  in  25  Rock  Samples . 43 

16.  Discrimination  of  Acidic  Silicates  Near  Mill  Creek,  Oklahoma 

(Sand  Quarry) . 44 

17.  Analog  Infrared  Images  of  Flight  Line  1,  Section  A,  Near 

Pisgah  Crater,  California . 45 

18.  Analog  Infrared  Images  of  Flight  Line  2,  Section  A,  Near 

Pisgah  Crater,  California . 46 

19.  Analog  Infrared  Images  of  Flight  Line  2,  Section  B,  Near 

Pisgah  Crater,  Californi-; . 47 

20.  Analog  Visible  and  Nec.r -Infrared  Images  of  Flight  Line  2, 

Section  A,  Near  Pi  jgah  Crater,  California . 50 

21.  Analog  Visible  apj  Near-Infrared  Images  of  Flight  Line  1, 

Section  A,  Nea.-  Pisgah  Crater,  California . 51 

ix 


- WILLOW  RUN  LABORATORIES  _ 

22.  Parabolic  Reflectometer  Optical  Schematic . . 54 

23.  Infrared  Spectra  of  Igneous  Silicate  Rocks . 56 

24.  Spectral  Variations  Resulting  from  Measurements  Taken  at 

Different  Locations  on  the  Samples . 63 

25.  Comparison  of  Direct  Spectral  Emittance  Measurements  (Solid 

Line)  with  Spectral  Emittance  Curves  Derived  from  Kirchoff's 

Law  and  Spectral  Reflectance  Measurements  (Dashed  Line) . 63 

26.  Spectral  Emittance  of  Iron  Rust . 63 

27.  Temperature-Corrected  Ratio  Recognition  Map  of  Silicate  Rocks 

Near  Pisgah  Crater,  California . 65 

28.  Mean  Reflectance  of  Three  Lavic  Phases  of  Pisgah  Crater  Basalts  ...  66 

TABLES 

1.  Spectral  Bands  of  Willow  Run  Laboratories'  Multispectral 

Scanner  System . . 

2.  Six  Spectral  Channels  Used  for  Recognition  of  Land-Use 

Features . . 

3.  Area  and  Percentage  of  Total  Area  for  Land-Use  Classes  in  the 

Region  Near  Milton,  Ontario . . . 40 

4.  Chemical  Composition  of  Igneous  Silicate  Rock  Samples . 53 

5.  Calculated  Thermal  Infrared  Ratios  by  a  Detector  at  an  Above- 

Ground  Altitude  of  1  km  on  a  Clear,  Dry  Summer  Day . 60 


X 


WILLOW  RUN  LABORATORIES 


TUNNEL-SITE  SELECTION  BY  REMOTE  SENSING  TECHNIQUES 

Final  Report 

1  April  1971  Through  4  March  1972 


1 

SUMMARY 

The  purpose  of  this  contract  is  to  evaluate  signal-enhancement,  data -reduction  techniques 
for  an  airborne  system  which  would  utilize  simultaneously  aerial  photographs,  synthetic  aper¬ 
ture  radar,  and  multispectral  scanning  systems  (the  latter  operating  in  the  0.32-  to  IS.S-pm 
spectral  region).  The  data  acquired  from  such  a  system  are  to  be  utilized  in  the  construction 
of  both  surface  and  subsurface  geologic  maps  which  could  be  used  to  predict  subsurface  con¬ 
ditions  in  possible  tunneling  sites.  Ground-truth  investigations  are  to  be  conducted  to  verify 
analyses  of  the  remote  sensing  data. 

For  these  evaluations,  primarily  we  are  using  recently  acquired  remote  sensing  imagery 
collected  by  Willow  Run  Laboratories.  We  are  presently  improving  and  refining  our  techniques 
for  collecting  and  analyzing  remote  sensing  data.  Recent  advancements  in  the  area  include: 

(1)  the  development  of  Vincent's  image  ratioing  method,  which  utilizes  multispectral  imagery 
to  detect  differences  in  the  SiOg  and  iron  content  of  geologic  formations;  and  (2)  development 
of  high-resolution  L-band  radar,  which  allows  us  to  acquire  like-  and  cross -polarized  images 
of  both  X-  and  L-band.  These  developments  in  remote  sensing  techniques  appear  to  be  applica¬ 
ble  to  the  identification  of  geologic  features.  Such  new  techniques  may,  in  fact,  be  required 
before  the  use  of  remote  sensing  in  tunne’  -site  selection  is  feasible. 

A  literature  search  was  performed  on  the  application  of  airborne  radar  and  multispectral 
scanning  systems  to  geologic  studies  in  order  to  avoid  duplication  of  effort,  to  build  upon  work 
performed  elsewhere,  and  to  acquaint  ourselves  generally  with  the  problems  involved  in  using 
remote  sensing  for  geologic  purposes.  We  have  now  acquainted  ourselves  with  much  of  the 
significant  work  in  the  geologic  applications  of  remote  sensing  and  with  most  sources  of  current 
contributions. 

At  Willow  Run  Laboratories,  compositional  remote  sensing  (the  determination  of  differences 
in  lithology)  has  been  advance  1  by  two  new  techniques  of  image  ratioing  performed  upon  multi¬ 
spectral  data.  These  techniques  allow  the  discrimination  between  silicate  and  nonsilicate 
geologic  formations  and  the  determination  of  the  amount  of  ferric  oxide  in  given  geologic  fea¬ 
tures.  Reflectance  and  emittance  calibrations  have  been  made  on  rock  samples  to  refine  and 
quantify  these  methods  further,  and  the  laboratory  tests  for  SiOg  have  confirmed  the  results 
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obtained  in  the  field.  Surficial  ferric  oxide  deposits  do  not  appreciably  alter  the  spectra  of 
underlying  rock  in  the  7-  to  14 -pm  range. 

The  geology  and  hydrology  of  a  glacial  drift  region  were  studied  through  the  data -processing 
techniques  of  level  slicing,  two-channel  ratioing,  and  multichannel  statistical  pattern  recogni¬ 
tion.  In  addition,  studies  of  thermal  variations,  aerial  photography,  and  radar  imagery  were 
used  for  this  study. 

Several  coherent  optical  processing  techniques  were  applied  to  available  radar  imagery  in 
an  attempt  to  reduce  the  amount  of  total  information  and  enhance  those  features  applicable  to 
the  identificatioii  and  Interpretation  of  geologic  features.  The  most  successful  technique  in¬ 
volves  an  image -enhancement  method  which  removes  the  low  spatial  frequency  portion  of  the 
image's  spectrum  by  means  of  a  spatial  filter.  This  method  resulted  in  the  capability  to  delin¬ 
eate  more  easily  the  structural  edges  and  lines  separating  differently  reflecting  terrain  regions. 

During  this  period,  we  have  also  interpreted  high-resolution  L-band  radar  imagery  along 
with  high-resolution  X-band  imagery.*  With  like-  and  cross -polarized  returns  in  each  band, 
we  have  the  capability  for  both  dual-spectral  and  dual -polarity  returns.  Interpretation  of  these 
radar  images  from  one  geologic  region  revealed  differences  in  the  imagery  related  to  the  texture 
of  the  surface  and  to  the  substratum  under  a  thin  surface  layer.  The  differences  shown  in  the 

L-band  and  X-band  imagery  designate  those  regions  for  which  ground-truth  investigations  should 
be  initiated. 

We  have  accurately  quantified  our  application  of  compositional  remote  sensing.  With  the 
use  of  compositional  remote  sensing  and  the  addition  of  multispectral  and  multipolar  radar, 
we  presently  anticipate  that  remote  sensing  will  be  an  economically  useful  technique  to  guide 
field  investigations  of  potential  tunnel  sites. 

No  remote  sensing  data  were  gathered  under  this  contract.  Although  both  radar  and  multi¬ 
spectral  data  were  available  for  several  sites,  no  site  was  suitable  for  a  case  study  of  remote 
sensing  for  geologic  areas  with  crystalline  rocks.  Fortunately,  data  from  a  glacial  drift  area 
were  recently  made  available  to  us.  Several  kinds  of  processing  were  performed  on  these  data 
to  aid  in  geologic  and  hydrologic  reconnaissance  for  tunneling. 

In  the  subsequent  contract  period,  it  is  anticipated  that  multispectral  data  will  be  obtained 
from  a  region  containing  a  wide  variety  of  crystalline  and  other  rock  outcrops.  A  test  site  in 
the  Black  Hills  was  chosen  for  a  study  of  tunneling  problems.  Radar  imagery  of  this  area  is 
available  from  other  sources.  We  hope  to  provide  a  test  case  to  evaluate  more  clearly  the  con¬ 
tribution  of  remote  sensing  to  tunnel -site  selection. 


^Specific  examples  of  the  radar  imagery  cannot  be  discussed  in  this  report  because  of 
security  and  proprietary  reasons. 
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2 

REMOTE  SENSING  FOR  TUNNEL-SITE  SELECTION 

2.1.  INTRODUCTION 

Multispectral  scanners  and  radar  systems  can  provide  geologic  and  terrain  information 
useful  for  exploration  and  planning  purposes.  The  potential  of  remote  sensing  systems  results 
from  the  new  kinds  of  information  collected,  the  large  amounts  of  data  involved,  and  their  com¬ 
patibility  with  modern  processing  facilities.  This  report  describes  the  results  of  a  one-year 
effort  to  process  and  analyze  existing  multispectral  and  radar  data  to  obtain  information  useful 
in  the  selection  of  tunneling  sites  and  in  planning  for  tunnel  construction.  Modern  remote 
sensing  imagery  is  used  to  determine  local  and  regional  geologic  and  hydrological  conditions, 
and  examples  of  high-altitude  aerial  photography  are  included  and  discussed. 

2.2.  IMAGING  REMOTE  SENSORS 

Imaging  sensor  systems  (cameras,  multispectral  scanners,  and  radars)  are  capable  of 
recording  only  the  surface  or  near  surface  of  the  terrain.  Thus,  the  problem  of  remote  sensing 
of  geologic  conditions  for  tunnel-site  selection  is  one  of  recording  and  interpreting  the  occur¬ 
rence  and  nature  of  surficial  materials  and  conditions  as  an  indicator  of  the  subsurface  struc¬ 
tures.  In  many  cases,  highly  accurate  predictions  of  the  geology  and  hydrology  of  an  area 
can  be  made  through  the  use  of  remote  sensor  imagery.  Structural  faults  and  fissures  have 
distinctive  appearances  on  imagery,  and  rock  types  are  often  identifiable  by  their  texture  and 
tone.  In  other  cases,  however,  because  of  dense  mantles  of  sediments  and/or  vegetation,  rela¬ 
tively  little  can  be  learned  about  subsurface  conditions  from  surface  images.  For  most  geo¬ 
graphical  areas,  remote  sensor  data  can  supply  much  useful  information  for  planning  tunnel 
construction,  but  seldom  all  that  is  lequired.  In  no  sense  can  remote  sensing  techniques  re¬ 
place  the  ground-survey  programs  currently  carried  out  in  advance  of  decisions  on  tunnel 
location.  Remote  sensing  can,  however,  reduce  ground-surveying  efforts  by  providing  early 
indications  of  potential  problem  areas  and  n.onitoring  surficial  changes  during  and  after  con¬ 
struction. 

To  obtain  surface  information  concerning  subsurface  conditions,  remote  sensor  data  can 
be  analyzed  in  two  ways:  (1)  by  an  examination  of  the  characteristics  and  occurrence  of  rock 
outcroppings,  wherein  rock  types  and  their  location  provide  direct  lithological  information 
concerning  the  site;  and  (2)  by  a  study  of  the  weathered  zone  of  the  earth's  mantle  (e.g.,  the 
vegetation  patterns,  soil  and  sediment  trends,  physiographic  features,  hydrologic  conditions, 
and  land-use  practices).  This  report  incorporates  recent  research  on  remote  detection  and 
mapping  of  rock  outcrops  (sponsored  under  this  and  related  contracts)  and  much  of  the  current 
state  of  the  art  in  the  remote  recording  and  mapping  of  nonlithologic  terrain  features. 
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2.3.  THE  ROLE  OF  REMOTE  SENSING  IN  TUNNEL-SITE  SELECTION 

In  general,  the  reliability  of  remotely  sensed  information  concerning  probable  subsurface 
conditions  decreases  with  increasing  depth  below  the  terrain  surface.  Necescary  information 
for  tunnei-site  selection  includes  knowledge  of  surface  and  ground  water  flow,  surficl.'il  sedi¬ 
ments,  and  even  land  use  in  the  form  of  building  locations  and  existing  surface  transportation 
routes.  Critical  sections  in  both  near -surface  and  deep  tunnel  construction  are  often  around 
portals.  Here  detailed  knowledge  of  slope  materials,  their  stratification,  engineering  proper¬ 
ties,  and  drainage  are  essential.  Knowledge  concerning  the  location  and  transportation  of 
materials  for  actual  tunnel  construction  can  also  be  useful. 

Szechy  (Ref.  1)  provides  the  following  list  of  rock  properties  of  interest  in  geologic  re¬ 
connaissance  for  tunnel  sites.  The  items  with  asterisks  are  geologic  or  lithologic  properties 
for  which  remote  sensing  can  contribute  information. 

*(1)  The  orientation  of  rock  stratification  (whether  horizontal,  sheet-like,  moderately 
inclined,  steeply  sloping,  reserved,  or  overfolded) 

♦(2)  The  thickness  of  individual  layers,  the  regularity  of  the  sequence  rock  layers,  or 
changes  in  mountain  types 

*(3)  Mineraiogical  composition  (detrimental  components) 

(4)  The  crystalline  structure  of  rocks  (uniformly  grained,  porphyrltlc,  etc.) 

(5)  The  bonds  between  the  individual  grains  (strong,  weak,  direct,  indirect) 

*(6)  The  hardness  of  the  rocks,  whether  they  would  allow  for  drilling  or  blasting 

♦(7)  The  structural  form  of  rocks  (massive,  stratified,  or  shaly) 

*(8)  Internal  structure  (whether  solid  or  porous,  ’-Ith  closed  or  open  voids) 

♦(9)  Deformations  suffered  during  the  origin  process  (cleavages,  crushed  zones,  faults) 
or  other  effects  (such  as  weathering,  mylonitizatlon,  or  kaolinlzation) 

(10)  The  probable  bearing  and  tensile  strength  of  the  mountain  at  various  tunnel  sections 

(11)  The  stability  of  the  mountain,  the  character  and  magnitude  of  probable  rock  pressure 

(12)  The  bulk  densities  and  dead  weights  of  component  rocks 

(13)  The  anticipated  durability  of  various  rock  types  to  be  penetrated,  the  length  of  en¬ 
trance  sections  to  be  lined  with  regard  to  the  danger  of  frost  effects 

(14)  The  depth  and  composition  of  cover  above  each  point  of  the  tunnel  for  each  rock 
constituent 

(15)  Temperature  conditions  affecting  the  mountain 

•(16)  Hydraulic  condltlono  at  tiie  construction  site  and  its  environment^ 


^Springs  and  watercourses  in  the  vicinity  should  be  observed  continuously  for  years  In 
order  to  obtain  evidence  concerning  the  amount  of  water  flowing  into  the  site  and  that  draining 
off.  The  water  should  also  be  tested  for  harmful  constituents. 
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(17)  Tho  possibility  of  the  occurrence  of  harmful  gases 

(18)  The  susceptibility  of  structures  to  earthquakes  and  artificial  vibrations 

*(19)  Surface  formations 

(20)  Safety  against  escaping  air  in  anticipation  of  operations  involving  compressed  air 

(21)  Hazards  to  structures  and  especially  to  entrance  portals  by  forces  of  mture  (e.g., 
slides,  rock  falls,  avalanches,  or  slumping) 

Of  these  twenty -one  geologic  features  applicable  to  the  selection  of  potential  tunneling  sites, 
remote  sensing  techniques  can  contribute  data  In  at  least  ten  categories. 

The  importance  of  the  information  that  remote  sensing  systems  can  provide  is  emphasized 
by  the  safety  hazards,  time  del.nys,  and  cost  Increases  caused  by  geologic  surprises.  Because 
basic  geologic  Information  was  Incomplete  before  tunneling  was  begun,  the  cost  of  some  tunnels 
has  been  multiplied  many  times  and  completion  delayed  for  years.  Although  the  Increasing 
sophistication  of  conventional  geologic  and  geophysical  reconnaissance  h.as  reduced  the  problem 
somewhat,  it  has  not  been  eliminated.  The  recent  difficulties  encountered  In  constructing  the 
Straight  Creek  Tunnel  on  Interstate  70  in  Colorado,  for  example,  have  resulted  In  a  substantial 
lncre.ise  In  cost  and  a  significant  delay  In  anticipated  completion. 

Aerial  photography  is  one  type  of  remote  sensing  that  has  been  used  conventionally  In 
geologic  reconnaissance  for  about  50  years.  Geochemistry,  measurements  of  radioactivity, 
and  geobotany  have  also  been  introduced  into  conventional  geologic  reconnaissance  (as  distinct 
from  geophysical  reconnaissance) during  this  century.  Multlspectral  scanning  and  sldelooklng 
airborne  radar  systems  are  two  of  the  newer  remote  sensing  systems  applicable  to  the  col¬ 
lection  of  geologic  data,  and  the  applications  of  these  systems  *0  geologic  investigations  are 
still  being  analyzed  and  evaluated. 

The  economic  feasibility  of  utilizing  these  newer  remote  sensing  systems  in  geol'?<^lc 
reconnaissance  has  been  questioned,  and  no  precise  figures  are  as  yet  available.  In  this  regard, 
however,  one  must  balance  the  Impact  of  the  additional  information  that  such  systems  can  pro¬ 
vide  against  the  safety  hazards,  the  scheduling  extensions,  anr  the  cost  Increases  Incurred  In 
construction  as  a  result  of  inadequate  geologic  information  about  the  tunnel  site.  The  following 
points  should  be  considered: 

(1)  The  amount  and  variety  of  remote  sensing  data  presently  avall.able  to  U.  S.  govern¬ 
mental  agencies  without  charge  is  large  and  growing  larger.  The  Strategic  Air  Command  h.as 
Imaged  approximately  15%  of  the  United  States  by  sidelooking  airborne  radar,  and  the  data 
are  available  from  the  Goodyear  Aerospace  Comp.iny  for  the  asking.  NASA  also  maintains 
three  remote  sensing  aircraft  (RD57E,  NP3A.  and  NC130B)  with  visible.  Infrared,  passive 
microwave,  and  .active  microwave  sensing  systems.  These  aircraft  are  devoted  to  »>itherlng 
remote  sensing  da»a  for  governmental  or  institutional  users,  and  NASA  maintains  a  continuously 
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growing  library  of  remote  sensing  data.  In  addition,  the  first  of  the  Earth  Resources  Technology 
Satellites  (ERTS-I)  was  launched  In  July  1972.  Carrying  remote  sensors,  these  satel¬ 
lites  are  coordinated  with  underflights  for  specific  research  projectn  and  are  gathering 
earth  resources  data  for  the  public  domain.  The  Earth  Resources  Observation  System  (EROS) 
program  will  maintain  a  center  In  Sioux  City,  Iowa,  to  disseminate  these  data. 

(2)  Modern  structural  geology  requires  a  large  regional  overview,  since  geologic  inferences 
are  made  in  the  context  of  regional  structures.  Remote  sensing  is  synoptic  and  is  thus  the  most 
economic  means  of  obtaining  regional  surface  features. 

(3)  Remote  sensing  Is  a  reconnaissance  tool.  The  basis  for  reconnaissance  is  the  identifi¬ 
cation  of  areas  for  intensive  study.  The  geologic  data  presently  collected  by  remote  sensing 
systems  are  largely  indicative  rather  than  definitive,  and  if  remote  sensing  were  the  only  in¬ 
vestigative  method  available,  the  lack  of  complete  definitiveness  would  be  a  serious  drawback. 
However,  when  used  as  a  complement  to  other  methods,  these  data  can  be  highly  Informative. 
Remote  sensing  systems  display  every  portion  of  the  surface  that  is  overflown  and  serve  as  a 
guide  for  the  field  geologist.  The  designation  of  a  distinctive,  although  undefined,  feature  Is 
highly  useful  In  reconnaissance.  The  fact  that  the  data  show  an  outcrop,  a  dike,  an  unexplained 
relief,  or  a  texture  Is  valuable  Information,  since  subsequent  mapping  expeditions  can  be  di¬ 
rected  to  pertinent  features.  The  time  required  for  geologic  mapping  is  reduced,  and,  In  addition, 
features  may  be  discovered  which  would  otherwise  be  overlooked. 

(4)  Remote  sensing  collection  and  data  processing  techniques  are  constantly  being  de¬ 
veloped  to  enable  us  to  define  geologic  features  more  clearly,  with  NASA  ?.nd  other  govern¬ 
mental  agencies  devoting  much  time  and  effort  toward  this  goal.  Remote  sensing  instrumenta¬ 
tion  and  data  analysis  methods  are  also  being  developed  for  application  to  other  disciplines, 
such  as  agriculture.  Ice,  soil,  and  land -use  studies,  and  the  new  instrumentation  and  methods 
that  have  evolved  are  being  applied  to  geologic  problems.  The  three  basic  contributions  to  the 
state  of  the  art  discussed  In  this  report  are  examples  of  the  continuing  advancement  toward 
more  varied  as  well  as  more  definitive  data. 

(5)  Although  geologic  reconnaissance  Is  a  small  portion  of  the  expense  of  tunneling,  the 
validity  of  the  available  geologic  Information  often  has  a  catastrophic  impact  upon  the  expense 
and  schedule  of  the  tunnel.  Since  more  complete  geologic  reconnaissance  can  potentially  elim¬ 
inate  unexpected  cost  increases,  safety  hazards,  and  scheduling  extensions,  every  possible 
geologic  interpretative  method  must  be  used,  Including  information  provided  by  remote  sensing 
techniques.  Remote  sensing  has  been  shown  to  be  useful  In  detecting  fault  zones  and  surface 
drainage  patterns,  both  of  which  ran  impede  tunnel  construction  and  produce  additional  ex¬ 
pense. 
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In  KUinniary.  tht>  nowor  ronioto  koiksIhj;  ttn  lmlqnt's  (l.o..  ninlltK|u>rlial  N(Mi\r,lnK  !>nd  nyn- 
Ihi'ltc  nporluro  radar  KyHlomn)  aro  rapidly  tuvoiulii)'  oporatlonal  and  ran  provido  Information 
nsrful  to  Ihr  Kolrrllon  of  tnnnolint;  hIIoh.  Tho  potruital  of  llirsr  Ky8trn\8  rrnnitn  from  llu*  i\r\v 
kinds  of  Information  rollrrlod.  Ihr  larno  amomits  ol  data  Involved,  and  llieir  rompalllilllly  with 
modern  proresnhij;  farlllllen.  Th,>  ainonni  of  remote  sensing  data  presently  available  at  no 
cost  l8  lar^e  and  niwInK'  larger,  and  extensive  libraries  of  sorb  data  already  exist  and  ronlliuie 
to  expand.  S>irli  data  are  the  most  eronon\lr  means  i»f  providing  an  overview  of  l  eglonal  geo- 
loglr  feal>ires.  and  berause  *'f  widespread  development  t)f  lnstnimentatli<n  and  analysis  metliods 
In  a  inimbei’  of  disripllnes,  the  appllratlons  of  riMiiote  sensing  d.tia  are  rapidly  being  extended. 
Alllunigb  the  remote  sensing  data  often  do  not  definitively  Identify  a  parlirnlar  geologir  sirurinre. 
they  do  serve  an  Imporlani  reronnals.sanee  fonellon  by  designating  dislinrilve  features  whieb 
might  otbei  wlse  t)e  overlookt'd.  Onre  the  f.'atores  are  loratt'd.  they  ran  serve  as  a  guide  to  the 
field  geologist  and  redure  the  time  reqnirt'd  for  geologir  mapping.  Hemole  sensing  data  ran 
ronirllnile  Information  that  may  elli\)lnate  nnexperted  rost  Increases,  safi'ly  hazards,  and  srhed- 
)dlng  extensions.  1  he  additional  rosi  of  iittllz.tng  remoti*  sensing  terhniqnes  In  i*o'  selertlon  of 
l)innel  sites  Is  low  when  rompar«Hl  to  the  rost  of  ovmromlng  unexpeded  geologir  diffirnllles 
after  Innneling  has  begun.  Hemole  sensing  bas  been  shown  to  be  parlirnlarly  useful  In  delerling 
fault  zones  and  s)irfaee  drainage  rharartri Islirs,  tH»tli  Imporlani  geologic  features  In  lonnel 
ronstnirllon. 


2.d.  HKLATKI)  IdTKKATUHi:  SOUIK'li' 

In  the  past  two  decades  there  has  grown  a  large  Isnly  t)f  llteralure  ronrerned  with  gi'ologir 
and  g»H)grapblral  Interpretations  from  aerial  |diolographs.  r\irrenl  efforts  to  obtain  geologir 
Information  from  new  types  of  remote  sensors  most  be  reg-.irdi'd  as  an  outgrowth  and  I'xtenslon 
of  lerhniques  developed  In  phologeology .  and  some  familiarity  widi  the  literature  Is  nxpilred. 
Kxlensive  bibliographies  of  phologeology  references  are  ronlalmx'  In  Talor  el  at.  |2|  .md 
Mlllei'  and  Srhumm  |3|.  Included  here  Is  a  listing  of  some  less  well-known  Information  sources 
which  relate  to  remote  seiusor  systems  other  than  or  In  addlllon  to  camera  sysimns. 

IMulotdiledly  the  single  most  comprehensive  compendinm  of  reporls  coiwernlng  remote 
sensor  systems  (Imaging  and  noi\lmaglng)  and  Iheir  nonmilllary  appllcallons  are  found  In  the 
I’loceedlngs  of  the  International  Synqtosla  on  Hemole  Sensing  of  the  Knvironment.  riu'se  sym- 
|H)sla  have  been  beld  al  IH-mon  Intervals  In  Ann  ArUtr,  Michigan,  since  Febniai-y  l!Ui2.  While 
the  first  four  .symposia  were  primarily  concerned  with  the  rapidly  developing  remole  sensor 
buimology.  since  April  lUltH.  Ilu'  d.tle  of  the  Fifth  Symposium,  an  Incre.ising  mimbei-  of  papers 
have  d'acussed  appllcalloi\s  of  remote  sensing  to  the  earth  sciences.  The  riweedlngs  of  (he 
Seventh  lntern.itlon.it  SynqH)sium  contain  7d  papers  and  are  com|H)sixl  of  three  volumes  of  7S0 
pages  each.  Al  least  half  of  Inese  p.ipers  present  the  re.sults  of  research  appileatde  1 1  (he 
earth  sciences  (o  g«'olog> .  i  cdology.  hydrology,  and  geography. 
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A  second,  less  comprehensive,  source  o  information  on  current  capabilities  of  remote 
sensors  can  be  found  in  the  record  of  the  Annual  Earth  Resources  Program  Review  of  the 
National  Aeronautics  and  Space  Administration  (NASA).  Since  NASA,  in  cooperation  with  a 
number  of  other  government  agencies,  sponsors  much  of  the  research  In  remote  sensing,  the 
program  reviews  provide  a  useful  guide  to  trends  In  current  applications.  The  latest  review, 
the  fourth  is  composed  of  five  volumes,  including  one  on  geology  and  geography. 

Since  1970,  the  quarterly  Journal  for  Remote  Sensing  of  Environment  (American  Elsevier 
Publishing)  has  provided  a  source  for  publication  of  current  remote  sensing  research.  Also, 
the  monthly  journal  of  the  American  Society  of  Photogrammetry,  Photogrammetrlc  Engineering, 
frequently  contains  articles  on  current  remote  sensing  applications,  although  it  includes  a 
large  range  of  subjects  other  than  remote  sensing. 

A  number  of  sources  which  treat  specific  aspects  of  remote  sensor  applications  relevant 
to  tunneling  will  be  referenced  where  appropriate  throughout  this  report.  Many  of  these  are 
limited-edition  reports  submitted  to  sponsoring  agencies  by  Willow  Run  Laboratories. 

2.5.  WORK  PERFORMED  UNDER  THIS  CONTRACT 

This  investigation  represents  a  first  step  in  expanding  the  use  of  remote  sensing  in  geo¬ 
logic  reconnaissance  for  potential  tunnel  sites  beyond  the  traditional  and  operational  photo- 
geology,  which  is  now  an  established  discipline  (4).  Our  investigation  has  shown  that  multi- 
spectral  scanninr  and  sidelooking  airborne  radar  systems  not  only  can  provide  more  geologic 
and  hydrologic  information  than  would  be  otherwise  available,  but  also  can  make  the  available 
photogeologic  data  more  definitive.  Listed  below  are  five  areas  in  which  new  geologic  infor¬ 
mation  can  be  obtained  through  remote  sensing  techniques.  The  last  three  of  these  areas  have 
been  under  intensive  investigation  during  this  contract  period. 

(1)  Extent  of  Fault  Zones.  The  application  of  remote  sensing  techniques  to  the  identification 
of  faults  is  Important,  since  weathered,  brecciated,  and  water -saturated  fault  zones  are  one  of 
the  most  difficult  and  expensive  problems  encountered  In  tunneling.  Rowan  et  al.  [5]  have  de¬ 
lineated  such  zones  in  the  Arbucklc  Mountains  of  Oklahoma  by  means  of  predawn  far-infrared 
imagery,  which  reflected  a  drop  in  temperature  caused  by  evaporation  of  water,  with  a  greater 
amount  of  moisture  in  the  fault  zone  than  in  luo  surrounding  area.  Infrared  Imagery  collected 

in  the  late  morning  and  afternoon  did  not  display  these  f.ault  zones,  nor  did  aerial  photography. 

(2)  Crystalline  and  Sedimentary  Rock  Outcroppings.  When  the  thermal  inertia  of  different 
types  of  rocks  (such  as  granite,  dolomite,  and  limestone)  was  investigated,  it  was  found  that 
each  rock  type  h.’s  a  different  thermal  inertia,  and  therefore  each  responds  differently  to  diurnal 
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temperature  variations.  Rock  outcroppings  can,  as  a  result,  be  mapped  and  identified  by 
suring  thermal  inertia  at  several  different  times  during  the  diurnal  cycle  (Watson,  Ref.  6). 

(3)  Hydrology,  Geomorphology,  and  Surficial  Soils.  Thermal  time  vari.ations,  single¬ 
channel  processing,  two-channel  ratioing,  and  multichannel  statistical  pattern  recognition  have 
been  used  to  analyze  multispectral  data.  In  Section  3  of  this  report,  Thomas  Wagner  discusses 
three  types  of  auromatic  processing  of  multispectral  dabi  useful  in  the  delineation  of  geologic 
features;  these  include  single-channel,  two-channel,  and  multichannel  processing,  the  applica¬ 
tions  of  which  are  ea'h  discussed.  It  was  found  th.at  multispectral  data  provided  more  Informa¬ 
tion  about  the  hydrology,  geomorphology,  and  soils  in  a  glacial  region  with  relatively  shallow 
bedrock  than  could  be  determined  from  aerial  photographs.  Single -channel  processing  allowed 
the  identification  of  standing  water  areas,  which  indicate  subsurface  hydrologic  conditions 
relevant  to  the  choice  of  tunneling  site;  such  areas  were  not  identifiable  on  an  aerial  photograph 
of  the  area.  The  multispectral  d.ata  also  delineated  surficial  deposits,  residual  soils,  crop 
vegetation,  and  various  classes  of  both  urban  and  rural  hand  use — all  of  which  are  Important 
considerations  in  tunneling. 

(4)  Rock  Composition.  In  Section  4,  Robert  Vincent  describes  compost. ional  remote  sens¬ 
ing  in  the  field  and  presents  data  confirming  his  Laboratory  calculations.  By  uslngthe  two-channel 
ratioing  method  of  processing  multispectral  data,  he  has  been  able  to  discriminate  between 
silicate  and  nonsilicate  geologic  formations  and  to  determine  the  amount  of  ferric  oxide  in  a 
given  geologic  feature.  This  ratio  method  is  appllc.able  to  the  Identification  of  surface  outcrops, 
the  tracing  of  past  landslides,  and  the  determination  of  the  number,  yield,  and  pH  value  of  exist¬ 
ing  springs  and  water  courses. 

(5)  Lineations,  Surface  Topogr.aphy,  and  Textures.  Den  Drake  reviews  the  applications 
of  sidelooking  radar  to  geologic  investigations  in  Section  5,  and  he  considers  the  parameters 
that  most  significantly  affect  radar  returns.  As  shown  in  the  literature  (Dcllwlg  e.  al..  Rowan 
and  C.annon,  Refs.  7  and  8),  faults,  joints,  and  contacts  between  different  lithologic  units  com¬ 
monly  have  relatively  shallow  depressions  along  them  which  are  enhanced  by  the  shadowing 
effect  of  SLAR,  and  regional  fold  and  fracture  p.itterns  commonly  are  more  e.asily  identified 
on  SLAR  imagery  than  on  aerial  photographs.  Radar  Imagery  particularly  well  displays  the 
relationship  between  the  topography  and  the  underlying  geologic  structure.  By  means  of  a  new 
multiplexed  radar  system.  Willow  Run  L.iboratories  will  soon  be  able  to  collect  imagery  simul¬ 
taneously  in  both  X-band  (3  cm)  and  L-band  (30  cm).  The  simultaneous  collection  of  data  in 
the  two  wavelength  regions  wiU  allow  .automatic  processing  of  data  .and  the  r.atioing  o^  the  two 
bands,  a  capability  that  should  provide  more  total  geologic  information  about  a  given  region, 
particularly  in  regard  to  surface  textures. 
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In  the  work  performed  under  this  contract,  we  have  assumed  that  the  remote  sensors  Im¬ 
portant  for  geolof^ic  reconnaissance  were  aerial  photography,  multispectral  scanners,  and 
sidelooking  airborne  radar  systems.  This  assumption  ignores,  however,  other  airborne  sens¬ 
ing  techniques,  such  as  passive  microwave  sensor  systems,  gamma-ray  sensors,  magnetometers, 
gravimeters,  and  electromagnetic  sounding  systems.  In  Section  6  of  this  report,  Philip  Jackson 
discusses  the  advantages  to  be  gained  b>'  including  thes  3  additional  remote  sensors  in  geologic 
investigations. 


3 

MULTISPECTRAL  REMOTE  SENSING 

3.1.  MULTISPECTRAL  DATA 

Up-to-date  references  are  currently  available  concerning  .'.erial  camera  systems,  and  to 
a  lesser  extent,  radar  systems.  There  is,  however,  a  paucity  of  available  information  con¬ 
cerning  modern  multispectral  scanning  systems  and  their  applications.  This  section  presents 
a  brief  introduction  to  the  concepts  and  technology  of  multispectral  data  collection  and  process¬ 
ing.  For  more  details  the  reader  is  referred  to  Shay  (1970)  [9]  and  several  state  of  the  art 
reports  prepared  by  Willow  Run  Laboratories;  in  particular,  see  Hasell  and  Larsen,  196i>  |10), 
Nalopka,  1970  (11),  Mallla  et  .al.,  1971  (12),  N.ilepka  et  al.,  1971  (13),  and  Erickson  and  Thomson, 
1971  (14). 

3.2.  THE  MULTISPECTRAL  CONCEPT 

The  nncept  of  multispectral  discrimination  is  relatively  simple.  The  reflectance  and 
emittance  of  radiant  energy  by  an  object  arc  wavelength  dependent  and  usually  specific  for  the 
materials  and  conditions  of  that  object  (15).  In  other  words,  different  objects  have  different 
reflectance  and  emittance  characteristics  which  vary  uniquely  over  the  electromagnetic  spec¬ 
trum.  By  collecting  information  in  a  number  of  spectral  bands  and  by  knowing  or  discovering 
the  unique  spectral  characteristics  of  objects  or  conditions,  it  is  possible  to  identify  these  and 
to  distinguish  one  from  another. 

An  image  formed  from  a  broad  spectral  r.ange  may  be  adequate  to  distinguish  objects 
which  have  grossly  different  radiant  characteristics  in  that  range.  For  example,  most 


‘The  term  "multispectral"  is  often  confused  with  the  term  "multiband"  in  remote  sensing 
literature.  We  regard  "multispectral"  d.ata  as  data  recorded  in  a  number  of  wavelength  ranges 
through  a  single  aperture.  Consequently,  all  of  the  data  are  in  exact  spatial  registration. 
"Multiband"  data  are  usually  the  product  of  multiple  aperture  camera  systems,  which  do  not 
assure  the  simultaneity  and  registration  of  data. 


10 


WILLOW  RUN  LABORATORIES 


photO|;raphic  films  arc  scns:(ivc  to  all  of  the  visible  spectrum,  l.c.,  to  the  0.4-  to  0.7-/im 
wavelength  range.  However,  by  comparing  two  or  more  images  of  the  same  object  made  in 
different  regions  of  the  spectrum,  the  chances  of  being  able  to  discriminate  and  identify  that 
object  arc  increased.  As  the  number  of  spectral  bands  is  increased,  additional  information  is 
available  about  the  object  of  interest.  Also,  the  potential  for  discriminating  between  different 
objects  Is  increased  with  the  addition  of  spectral  information.  If  the  spectral  characteristics 
of  objects  arc  not  specifically  known  (generally  the  ease).  It  is  an  advantage  to  have  a  number 
of  spectral  images,  because  one  or  more  will  prob.ably  show  tone  differences  between  two  or 
more  objects  better  than  others.  Thus,  the  inherent  ability  to  distinguish  between  similar  ob¬ 
jects  is  enhanced  with  an  increase  in  the  range  and  amount  of  spectral  information.  The  b.aslc 
strategy  in  multlspectral  sensing  calls  for  the  accurate  detection  of  spectnal  characteristics 
produced  by  objects  of  interest  in  a  scene  (with  use  of  as  many  spectral  bands  as  possible)  and 
the  correlation  of  these  spectral  characteristics  with  the  properties  or  conditions  of  the  objects. 

3.3.  A  MULTISPEC TRAL  SCANNING  SYSTEM 

To  take  advantage  of  the  spectral  information  of  a  scene,  the  reflected  and  emitted  radi¬ 
ation  must  be  sampled  in  a  number  of  discrete  bands  In  order  to  get  accurate  representation 
of  the  energy  distribution.  For  this  purpose,  an  airborne  optical -mechanical  sririi;cr  system, 
using  reflective  optics,  has  certain  advantages  over  multiple  aperture  systems.  First,  the 
spectral  information  is  recorded  in  electrical  form  on  computer-coPipatiblc  magnetic  tape, 
and,  therefore,  certain  limitations  in  data  handling  and  processing  Inherent  in  photographic 
processes  are  avoided.  Second,  the  detectors  in  scanner  systems  are  not  limited  to  the  rather 
narrow  spectral  range  over  which  conventional  films  are  sensitive.  Modern  scanning  systems 
collect  data  over  a  0.32-  to  14.0 -pm  spectral  range  (Fig.  1).  However,  because  of  atmospheric 
absorption,  surface  radiation  cannot  be  recorded  In  the  2.6-  to  3.5-pm  and  the  5.5-  to  8.0-pm 
wavelength  ranges. 

Basically,  an  optical-mechanical  scanner  consists  of  a  telescope,  a  rotating  mirror,  and 
some  sort  of  radiation -detector  configuration.  The  telescope  focuses  the  radiation  from  a 
small  portion  of  the  scene  onto  a  detector.  Figure  2  shows  a  schematic  diagram  of  such  a 
scanning  system.  The  rotating  mirror  causes  the  scanner's  field  of  view  to  sweep  rapidly 
across  the  flight  path.  The  direction  of  the  sweep,  or  scan,  is  perpendicular  to  the  direction 
of  flight  of  the  aircraft.  The  forward  velocity  of  the  plane  and  the  rotation  of  the  mirror  are 
such  that  a  complete  record  of  the  scene  is  sequentially  recorded  in  a  series  of  contiguous 
scans.  The  resulting  scan  pattern  is  similar  to  a  television  raster,  except  that  it  is  continuous 
for  the  entire  flightiinc.  The  rotation  rate  of  our  scanning  mirror  is  60  revolutions  sec.  At 
an  altitude  of  3000  ft,  for  example,  each  scanline  represents  a  ground  swath  10  ft  wide  and 
almost  1  mile  long  (00°  field  of  view).  In  this  manner,  the  spectral  data  of  a  scene,  composed 
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FIGURE  2.  MULTISPECTRAL  DATA  COLLECTION 
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of  many  parallel  scanlines,  are  obtained  in  the  length  of  time  it  takes  for  the  aircraft  to  fly 
over  the  area. 

Modern  multispectral  systems  employ  radiation-dispersing  prisms  behind  the  single  aper¬ 
ture.  After  the  prism  disperses  the  radiation  into  separate  wavelength  regions,  it  is  trans¬ 
mitted  by  an  array  of  relay  optics  to  a  number  of  radiation  detectors.  Photodetectors  are  used 
in  the  reflective  wavelengths  and  special  cooled  metal  detectors  in  the  infrared  wavelengths. 
These  detectors  convert  the  incoming  radiant  energy  to  electrical  signals,  which  are  then  amp¬ 
lified  and  recorded  as  voltage  levels  on  multitrack  magnetic  tape.  Radiation  calibration  is 
automatically  recorded  coincident  with  the  scene  radiation.  Calibration  information  includes 
synchronization  pulses,  upward-looking  sun-sensor  data,  two  reference  lamps,  and  two  thermal 
reference  plates.  The  fact  that  up  to  twelve  spectral  bands  are  recorded  simultaneously  for  a 
single  scene  is  important.  The  spectral  information  of  each  of  the  bands  is  in  exact  registration 
with  the  others,  a  fact  which  makes  possible  the  subsequent  computer  processing  of  the  data. 

Willow  Run  Laboratories'  current  multispectral  scanner  system  (M-7)  went  into  operation 
in  June  1971  (replacing  the  M-5  system,  which  had  been  operational  since  April  1966).  The 
spectral  bands  which  can  be  recorded  with  this  system  are  listed  in  Table  1.  Because  of  the 
magnetic  tape's  channel  limitations,  it  is  not  possible  to  record  more  than  12  of  these  spectral 
bands  simultaneously. 

3.4.  AUTOMATIC  PROCESSING 

The  large  amount  of  information  collected  by  the  multispectral  scanner  system  requires 
that  automatic  processing  techniques  be  used.  There  r  e  540  resolution  elements  per  scanline, 
and  60  scanlines  are  recorded  each  second,  which  amounts  to  some  32,400  recorded,  quantitative 
radiation  samples  per  second.  A  typical  15-mile  flightline  will  be  covered  in  6.6  min,  resulting 
in  over  300  million  samples,  each  sample  having  12  voltage  levels  representing  the  effective 
radiance  in  12  recorded  spectral  bands.  The  conventional  tools  of  image  interpretation  are  not 
capable  of  handling  this  amount  of  data  in  a  consistent  and  quantitative  manner.  Human  exam¬ 
ination  and  interpretation  of  sensor  Imagery,  using  accumulated  familiarity  with  specific  areas 
and  training  with  special  keys,  require  much  more  time  than  that  needed  to  obtain  the  original 
data  and  can  not  approach  the  information  potential  inherent  in  modern  sensor  data.  These  tra¬ 
ditional  techniques  are  neither  adequate  for  analyzing  large  amounts  of  data  in  a  short  time 
nor  capable  of  providing  results  that  are  not  biased  by  the  relative  skills  and  patience  of  the 
image  interpreter.  What  is  required  is  a  process  to  reduce  the  information  input  rate  re¬ 
quired  of  the  individual  interpreter,  i.e.,  to  insert  a  filter  between  the  sensor  system  and  the 
person  who  analyzes  the  data.  The  filter  should  be  capable  of  rejecting  much  of  the  undesired 
information  ("noise")  so  that  it  is  not  presented  to  the  interpreter.  This  data -filtering  process 
must  be  able  to  discriminate  automatically  between  objects  and  materials  of  potential  interest 
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TABLE  1.  SPECTRAL  BANDS  OF  WILLOW  RUN  LABORATORIES' 
MULTISPECTRAL  SCANNER  SYSTEM* 


50%  Response  Points 
(Mm) 

Spectral 

Range 

0.32-0.38 

Ultraviolet 

0.41-0.48 

Violet 

0.46-0.49 

Blue 

0.48-0,52 

Blue-  Green 

0.52-0.54 

Green 

0.53-0.57 

Yellow- Green 

0.55-0.60 

Yellow 

0.58-0.64 

Orange 

0.62-0.69 

Red 

0.67-0.91 

Reflected  Infrared 

1.0-1.4 

Reflected  Infrared 

1.5-1. 8 

Reflected  Infrarc  d 

2.0-2. 6 

Reflected  Infrared 

4.5-5. 5 

Thermal  Infrared 

8.2-10.9 

Thermal  Infrared 

8.2-13.5 

Thermal  Infrared 

9.4-12.1 

Thermal  Infrared 

♦Additional  bands  may  be  recorded  through  the  introduction  of 
different  filter  systems  in  the  0.32-  to  13,5-;im  range. 
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and  those  potentially  not  of  interest  and  present  the  results  of  this  discrimination  in  an  easily 
interpretable  form. 

Digital  and  analog  computer  systems  providt  such  a  capability  for  processing  multispectral 
data.  The  analog  computer  (or  processor)  can  process  the  data  at  a  near  real-time  rate,  that 
IS,  at  almost  the  same  rate  that  the  scanner  collecte,!  the  original  data.  The  same  filtering 
operation  can  be  implemented  on  digital  compute.,.^,  but  at  a  sacrifice  of  speed  and  often  of  the 
amount  of  data  which  can  be  economically  handled.  Digital  computers  have,  however,  the  ad¬ 
vantage  of  providing  precise  quantitative  analysis  of  the  data. 

Three  types  of  processing  (filtering)  of  multispectral  data  can  be  implemented  by  com¬ 
puter:  (1)  single -channel  quantization  and  level  slicing;  (2)  two-channel  ratioing;  and  (3)  multi¬ 
ple-channel  statistical  pattern  recognition.  These  three  processing  techniques  are  described 
below  in  order  of  increasing  sophistication;  examples  and  applications  of  these  three  techniques 
are  discussed  in  subsequent  sections. 

3.4.1.  SINGLE -CHANNEL  PROCESSING 

Since  a  single  spectral  region  contains  scene i-adianoe  information  distinct  from  all  other 
regions,  it  is  often  useful  to  select  one  band  for  automatic  processing.  A  knowledge  of  the 
potential  information  in  individual  spectral  bands  allows  their  use  in  extracting  certain  kinds 
of  information  automatically.  The  objective,  then,  is  to  isolate,  consistently  and  objectively, 
radiance  differences  which  are  indicative  of  conditions  or  objects  of  interest.  Optimally,  the 
data  will  be  presented  ir.  a  format  readily  interpretable  to  the  earth  scientist. 

This  objective  is  accomplished  by  playing  the  recorded  signals  from  a  single  channel 
through  a  data  processor.  In  this  case,  the  processor  automatically  quantizes  the  signal  levels 
into  a  preselected  number  of  discrete  ranges.  Signal  voltages  are  subsequently  recorded  as 
belonging  to  one  of  a  number  of  voltage  ranges.  The  total  voltage  range  may  include  that  range 
recorded  on  tape  for  the  entire  scene  or  some  subset  of  that  range.  When  the  result  of  this 
quantization  process  is  printed,  a  scene  image  may  be  portrayed  by  a  number  of  distinct  tones 
(density  levels)  on  a  film  transparency,  or  each  discrete  range  may  be  separately  printed  as  a 
black  tone  on  an  otherwise  clear  film  transparency.  The  latter  type  of  image  display  is  known 
as  a  level  slice  or  separation.  Subsequently  a  number  of  quantized  separations  may  be  color- 
coded  and  overlaid,  one  with  another,  to  provide  a  color-coded  composite  image. 

The  result  of  this  type  of  processing  is  similar  to  standard  techniques  of  film  densitometry, 
except  that  the  data  are  taken  directly  from  magnetic  tape  and  the  results  are  portrayed  in 
image  form.  Some  of  the  difficulties  associated  with  working  with  nonlinear  film  and  film 
processing  variations  are  thus  avoided.  Also,  since  sensor  calibration  sources  are  also  quan¬ 
tized,  an  accurate  measure  of  radiation  levels  is  available.  One  can  use  this  technique  effectively 
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for  quantitative  temperature  mapping,  using  a  thermal  infrared  spectral  channel,  and  for 
surface-water  mapping  by  means  of  reflective  infrared  data.  (Examples  of  thermal  mapping 
and  surface  hydrology  are  included  in  Section  3.6.I.I.) 

3.4.2.  TWO-CHANNEL  PROCESSING 

With  two  channels  of  multispectral  data,  it  is  frequently  useful  to  divide  the  signals  in  one 
channel  by  those  of  another  channel  to  produce  a  new  ratio  image.  The  shorter  wavelength  band 
is  divided  into  the  longer  wavelength  band,  and  the  new  image  displays  differences  between  the 
two  original  images.  A  ratioed  image  essentially  eliminates  the  similarities  between  two  bands 
of  spectral  data.  Also,  the  resulting  ratio  (voltage  level)  is  a  measure  of  the  average  slope  of 
a  spectral  curve  which  spans  the  two  wavelength  ranges.  This  technique  is  found  to  be  useful 
for  obtaining  certain  types  of  information  from  both  reflective  and  thermal  multispectral  data. 
The  theory  and  applications  of  this  technique  were  included  in  the  report  by  Drake  et  al.,  1971 
[16],  and  further  examples  are  contained  in  Section  4  of  this  report. 

A  ratio  of  two  bands  of  multispectral  data  is  produced  by  means  of  an  electronic  processor, 
so  that  for  every  scene  resolution  element  the  signal  level  recorded  in  one  channel  is  divided 
by  the  signal  level  recorded  in  another.  Normalization  of  the  two  bands,  to  adjust  for  optimal 
electronic  signals  on  the  ratioed  image,  is  accomplished  by  means  of  reference  standards  re¬ 
corded  coincident  with  the  original  data.  The  product  is  a  film  transparency  in  which  film 
density  is  proportional  to  the  ratio  (quotient)  of  the  signals  recorded  on  the  original  two  chan¬ 
nels.  Usually  a  neutral  grey  tone  is  selected  to  represent  a  ratio  value  of  unity.  Departures 
from  unity  are  represented  as  lighter  tones  for  higher  values  (greater  than  one)  and  darker 
tones  for  values  less  than  one.  If  the  ratioed  data  are  electronically  recorded  directly  onto 
magnetic  tape,  they  can  be  further  processed  as  a  single  channel. 

Ratioed  images  from  data  collected  in  the  reflective  wavelengths  (<3.0  pm)  result  in  in¬ 
formation  concerning  differences  in  the  spectral  reflectance  of  objects.  Ratioed  images  in  the 
thermal  wavelengths  (>3.0  pm)  display  useful  data  concerning  the  emittance  characteristics 
of  objects  or  materials,  i.e.,  the  degree  to  which  their  thermal  radiation  characteristics  depart 
from  those  of  a  perfect  emitter  (a  so-called  blackbody). 

3.4.3.  MULTICHANNEL  PROCESSING 

Multichannel  processing  is  the  most  sophisticated,  and  potentially  the  most  effective  type 
of  multispectral  processing.  Essentially  the  technique  utilizes  computer -implemented  pattern- 
recognition  schemes  to  analyze  and  discriminate  objects  or  materials  within  a  scene.  If  one 
has  a  sample  of  object  classes  of  interest  within  the  scene,  he  needs  no  other  a  priori  spectral 
information  to  use  this  technique.  Both  analog  and  digital  computers  can  be  used,  and  often 
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extensive  statistical  analyses  of  the  spectral  c  haracteristics  of  a  scene  are  made  before  the 
results  arc  printed  out  in  image  form . 

There  are  two  phases  to  the  implementation  of  pattern- recognition  schemes.  One  is  the 
training  phase,  in  which  the  data  are  calibrated  and  corrected  for  illumination  and  other  types 
of  signal  variations.  Samples  of  objects  of  interest  are  located  on  a  visual  display  of  the  scene, 
and  these  areas  are  then  "gated"  in  such  a  way  that  only  signals  from  the  sample  areas  are  re¬ 
corded  in  the  computer's  memory.  These  signals  are  the  statistical  means  and  standard  devia¬ 
tions  for  the  sample  objects.  Th?  combined  mean  and  variance  for  all  channels  is  known  as  the 
spectral  signature  for  a  given  sample.  If  that  signature  represents  the  recorded  spectral  char¬ 
acteristics  of  a  sample,  it  is  assumed  that  the  signature  indicates  the  spectral  characteristics 
for  all  similar  members  of  that  class  of  objects.  Thus,  once  signatures  are  programmed  into 
the  computer,  the  operational  phase  consists  of  having  the  computer  analyze  the  spectral  char¬ 
acteristics  of  each  resolution  element  from  a  set  of  unknown  data  and  determine  which  of  the 
elements  is  sufficiently  similar  to  any  of  the  signatures  to  be  recognized  as  a  member  of  a 
known  class. 

For  the  computer,  similarity  is  defined  by  a  mathematical  decision  rule.  The  rule  most 
commonly  used  is  the  likelihood  ratio.  In  likelihood-ratio  processing,  each  resolution  element 
of  the  data  is  classified  as  "target"  or  "not  target"  by  noting  if  the  likelihood  ratio,  L,  is 
greater  than  or  less  than  some  threshold  value,  T.  In  simple  form,  the  likelihood  ratio  is; 

P^(T)P(T/S)  >  T  =  target 

^  <  T  =  no  target 

where  P^(T1  =  a  priori  target  probability 

P(T/S)  =  probability  of  targ'  •,  given  a  data  sample 
=  a  priori  probability  of  background  n 

P(Bn/®)  “  probability  of  background  n,  given  a  data  sample 

The  result  of  this  yes/no  type  of  decision  rule  is  a  computer  recognition  map.  This  map 
is  an  image  display  wherein  all  areas  recognized  as  being  similar  to  a  given  signature  are 
displayed  as  a  particular  symbol  or  color.  Areas  of  the  test  site  that  are  not  similar  to  any  of 
the  signatures  are  left  blank,  and  thus,  the  recognition  map  displays  only  materials  of  interest. 
An  example  of  this  type  of  processing  is  given  in  Section  3.6.3. 

3.5.  DATA-COLLECTION  PARAMETERS 

Important  considerations  for  planning  the  collection  of  multispectral  data  involve  the 
spectral  bands  to  be  recorded,  the  time  of  day  or  night,  and  the  altitude  and  direction  of  the 
flight  or  flights.  The  attention  given  to  these  considerations  often  determines  the  relative 
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success  of  the  remote  sensing  mission  and  the  subsequent  value  of  the  data.  A  'ourth  import¬ 
ant  consideration,  often  neglected,  is  the  extent  and  nature  of  supportin.t  ground-truth  investi¬ 
gations. 

Multispectral  data  collected  at  night  are  restricted  to  the  thermal  wavelengths  (^3  ^xm). 
Because  of  the  diurnal  cycle  of  surface  temperatures,  it  is  usually  advisable  not  to  collect 
thermal  data  until  several  hours  after  sunset.  Prior  to  this  time,  the  differential  rates  of 
cooling  of  terrain  surfaces  obscure  subtle  temparature  effects  which  are  diagnostic  of  certain 
terrain  conditions  (such  as  mlcrorellef,  vegetation  distribution,  or  coil  moisture).  From  an 
operational  point  of  view,  dawn  provides  a  convenient  time  for  thermal  data  collection,  in  that 
there  is  usually  sufficient  light  for  aircraft  pilots  to  navigate  without  ground- marker  lights. 
Ground  fog,  however,  can  seriously  degrade  predawn  and  dawn  Imagery,  and  flights  under  foggy 
conditions  should  be  avoided. 

Multispectral  missions  for  the  collection  of  data  in  the  visible  and  near-IR  wavelengths 
are  most  successful  during  mid-day  hours.  For  most  purposes,  low  sun  angles  reduce  the 
quality  of  the  imagery  by  Introducing  pronounced  terrain  and  vegetation  shadows  and  by  causing 
radiance  variations  resulting  from  bidirectional  reflectance  differences  for  objects  seen  from 
different  scan  angles.  For  the  latter  reason,  flighMines  should  be  oriented  as  much  as  possible 
to  insure  that  the  aircraft  flies  parallel  to  the  sun's  rays,  l.e.,  into  or  180°  from  the  direction 
of  the  solar  azimuth.  Bidirectional  reflectance  is  not  a  serious  problem  when  two-channel  ratio 
processing  of  the  data  is  anticipated,  but  it  can  present  proble.ps  for  multichannel  discrimina¬ 
tion  techniques.  Daytime  thermal  data  are  often  best  collected  during  the  mid- afternoon  hours, 
at  which  time  the  terrain  has  reached  a  relative  state  of  thermal  equilibrium. 

Flight  altitude  is  determined  by  the  trade-off  between  the  degree  of  Im.age  resolution  and 
the  total  area  of  coverage  desired.  Both  the  resolution-element  dimensions  and  the  total  area 
scanned  increase  linearly  with  alUtude.  In  Willow  Run  Laboratories'  scanner,  the  instantaneous 
field  of  view  (IFOV)  is  2  mrad;  i.e.,  a  ground  patch  2  ft  sq  will  be  recorded  at  1000- ft  elevation. 
The  scanner  also  has  an  effective  90°  scan,  which  means  that  at  1000 -ft  altitude  a  swath  some 
2000  ft  wide  will  be  recorded,  and  at  3000  ft  the  swath  width  will  be  more  than  1  mile.  How¬ 
ever,  local  variations  in  terrain  influence  somewhat  the  amount  of  area  imaged. 

3.6.  APPLICATIONS  OF  AUTOMATIC  PROCESSING  OF  MULTISPECTRAL  DATA 

This  section  presents  results  of  the  application  of  the  three  types  of  automatic  multi¬ 
spectral  data  processing  discussed  in  Section  3.4.  Examples  of  both  processed  and  unprocessed 
imagery  were  selected  from  existing  data  files  to  provide  information  thought  to  be  applicable 
to  tunnel- site  selection.  These  results  represent  but  a  portion  of  the  remote  sensing  techniques 
currently  being  developed  for  earth  resources. 
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3.6.1.  SINGLE-CHANNEL  APPLICATIONS 

Quantizing  of  a  single  ch.innel  of  multispectral  data  clearly  delineates  often  subtle  con¬ 
trasts  inherent  for  the  scene  in  the  selected  v.'avelength  band.  Single- channel  processing  has 
been  useful  in  the  discrimin.'ition  of  surface  water  and  in  thermal  mapping. 

3. 6. 1.1.  Surface  Hydrology 

For  delineation  of  surface  water,  a  single  channel  in  the  near-infrared  spectral  region 
was  used.  Standing  water  areas  were  discriminated  from  all  other  (background)  materials  by 
their  low  signal  levels  in  the  1.5-  to  1 .8- pm  wavelength  range.  This  data  channel  was  played 
through  a  digital  processor,  and  all  signals  (voltage  levels)  falling  below  a  selected  threshold 
limit  were  imaged  on  a  cathode  ray  tube  (CRT)  and  subsequently  recorded  by  a  filmstrip  printer. 
All  recorded  signals  were  tallied  by  a  digital  counter  linked  with  the  system. 

Computer- implemented  discrimination  of  standing  water  consists  of  setting  a  threshold 
level  for  separating  the  low  sign-J  levels  of  standing  water  from  those  of  other  scene  elements. 
All  signals  below  the  threshold  are  categorized  as  standing  water.  The  initial  threshold  value 
is  empirically  determined  by  adjusting  the  threshold  level  to  the  point  that  known  areas  of 
standing  water  are  adequately  imaged  on  a  CRT  display,  and  background  elements  are  not  re¬ 
corded. 

Figure  3  compares  three  images  from  a  flightline  which  traverses  areas  of  extensive 
wetlands — lakes,  bogs,  swamps,  and  artificial  impoundments.  As  the  figure  shows,  the  proc¬ 
essed  image  of  standing  water  provides  the  location  and  areal  extent  of  all  surface  water, 
while  the  0.52-  to  0.57-/im  image  (the  spectral  range  to  which  the  human  eye  is  most  sensitive) 
provides  little  contrast  for  discriminating  these  surface-water  areas. 

The  imagery  in  Fig.  3  is  from  a  portion  of  a  flightline  which  crosses  such  Pleistocene 
features  as  outwash  gravel,  till  moraines,  kames,  and  drumlins.  A  porous  dolomite  formation 
underlies  this  area,  often  at  depths  of  less  than  50  ft.  Surface  drainage  is  impeded  by  the  rolling, 
hummocky  glacial  topography,  and  there  are  no  surface  streams  evident  on  the  imagery.  As 
indicated  by  the  abundance  of  surface  bodies  of  water  and  lack  of  integrated  drainage,  the  prob¬ 
able  occurrence  of  subsurface  high-yielding  aquifers,  either  within  the  surficial  sediments  or 
the  underlying  bedrock,  would  present  serious  problems  for  rapid  tunnel  construction. 

DistriuuMon  of  surface  water  as  an  indication  of  subsurface  hydrologic  conditions  is 
further  illustrated  in  Fig.  4.  A  mosaic  of  surface-water  discrimination  images  shows  a  160- 
sq-mlle  area  in  southern  Ontario.  (For  comparison,  a  high-altitude  photomosaic  of  approxi¬ 
mately  the  same  area  is  shown  in  Fig.  5.)  The  watershed  outlined  in  the  center  of  the  image 
is  the  90- sq- mile  drainage  basin  of  East  and  Middle  Oakville  Creek,  which  flows  into  Lake 
Ontario  at  the  bottom  of  the  image.  Surface-water  patterns  show  a  marked  difference  for  the 
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(a)  Video  Image;  0.52  to  0.57  ftm 


fti!  Video  Image;  1.5  to  1.8  /im 


(c)  I  evel -Sliced  Image;  1.5  to  1.8  pm 


FIGURE  3.  COMPARISON  OF  TWO  VIDEO  IMAGES  AND  A  LEVEL-SLICED  IMAGE  FOR  THE 
DISCRIMINATION  OF  STANDING  WATER.  Imagery  collected  over  No.  6  Hlghu.ay,  between 

Hamilton  and  Guelph,  Ontario. 
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FIGURE  5.  HIGH-AL'  ITUDE  PHOTOMOSAIC  OF 
OAKVILLE  REPRESENTATIVE  BASIN.  ONTARIO. 
Appro.xIm.itc  are.i:  160  sq  miles. 
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upper  and  lower  portions  of  this  mosaic.  In  the  upper  portion,  numerous  lakes,  swamps,  and 
bogs  indicate  the  location  of  glacial  sediments  u'«derlain  at  shallow  depths  by  the  Anabcl 
dolomite,  rncntio'.ied  above.  Most  of  the  central  and  lower  portion  of  the  image  show  scattered 
small  ponds  and  the  meandering  pattern  of  Oakville  Creek.  This  area  is  composed  of  well* 
drained  clay  loam  tills,  underlain  by  soft  shales  and  sandstones  of  the  Ordovician  Queenston  for¬ 
mation.  Thu  winding  course  of  the  creek  indicates  that  drainage  is  slow  over  this  area,  which 
slopes  gently  to  the  Lake  Ontario  shore.  Most  ol  the  small  ponds  are  associated  with  farms. 
Digital  analysis  of  these  two  types  of  terrain  indicates  that  approximately  2.1%  of  the  surface 
underlain  by  the  massive  Anabel  dolomite  Is  standing  water,  whereas  the  sloping  shale  plain 
of  the  Queenston  formation  has  0.55%  surface  water.  Except  for  an  occasional  sand  lens,  one 
would  not  expect  to  encounter  tunneling  problems  associated  with  cither  rapid  surface  recharge 
to  the  watcrtable  or  high-yielding  aquifers.  Indeed,  water  supply  is  known  to  be  somewhat  of 
a  problem  in  this  area  as  a  result  of  the  general  lack  of  aquifers  and  the  high  rates  of  surface 
evaporation  associated  with  land  cultivation  (Chapman  and  Putnam,  1966,  Ref.  17). 

3.6.1. 2.  Thermal  Mapping 

For  a  number  of  years,  useful  hydrological  and  geologic  information  has  been  obtained 
from  single- channel  thermal  Images  (Fisher  et  al.,  1964;  McLerran  and  Morgan,  1965; 

Friedman  and  Williams,  196b;  and  Quade  et  al,,  1970 — Refs.  18-21).  To  some  extent,  recent 
applications  of  thermal  sensing  have  been  hindered  by  limitations  of  conventional  qualitative 
approaches  to  the  interpretation  of  such  Imagery.  Simplistic  concepts,  such  as  hot-spot  de¬ 
tection,  have  limited  utility  in  defining  natural  processes. 

Two  sets  of  thermal  images  (6.0  to  13.5  pm)  of  a  morainal  area,  the  Blue  Springs  Basin, 
in  southern  Ontario  illusti  KC  the  kind  of  general  land-use  Information  that  can  be  obtained 
from  single- channel  processing  of  thermal  imagery. 

Data  for  Fig.  6  werp  collected  at  alcut  noon  from  an  experimental  watershed  near  Acton, 
Ontario.  One  image  shows  the  thermal  da'a  as  continuous  tone  variations,  with  the  warm 
areas  light  in  tone  and  the  cooler  areas  relatively  darker.  The  second  image,  derived  from 
the  same  data,  is  a  calibrated  thermal  composite  in  which  four  thermal  regions  are  designated 
by  color  code.  The  grey- scale  image  reveals  considerable  detail  concerning  the  thermal  com¬ 
plexity  of  the  scene  but  no  quantitative  information.  In  many  respects,  this  image  appears 
similar  to  a  conventional  panchromatic  photograph  of  the  area.  Both  field  patterns  and  the 
cool  water  areas  are  evident.  The  color-coded  image  provides  a  somewhat  simplified  view  of 
this  area,  with  the  total  surface-temperature  range  for  the  scene  somewhat  in  excess  of  30°F. 
Some  of  the  fields,  coded  as  green  and  brown,  are  characterized  by  higher  temperatures  than 
the  rest  of  the  area  and  are  known  to  be  bore  soil;  the  two  other  thermal  ranges  correspond  to 
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Color  Code 
Blue  <  65.3°F 

Light  Blue  *  65.3°  to  TO.OOp 
Green  t  70.0°  to  75.6°F 
Brown  >  75.6°F 


(a)  Video  Image 


(b)  Calibrated  Thermal  Composite  Image 


FIGURE  6.  COMPARISON  OF  STANDARD  THERMAL  VIDEO  IMAGE 
WITH  CALIBRATED  LEVEL-SLICED  THERMAL  IMAGE;  DATA  COL¬ 
LECTED  AT  1200  HOURS.  Blue  Springs  Basin,  II  May  1971. 
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densities  of  vegetation — the  cooler  the  temperature,  the  denser  the  vegetation  (perhaps  a  re¬ 
sult  of  the  effects  of  transpiration). 

The  data  *or  Fig.  7  were  collected  at  dawn  (0400  hours)  on  the  same  day  as  those  for 
Fig.  6.  The  differential  heating  and  cooling  of  the  surface  has  produced  considerably  different 
thermal  contrasts  between  the  mid-day  and  dawn  imagery.  The  continuous  tone  image  shows 
many  subtle  differences  between  fields;  the  lighter  toned  water  areas  provide  the  only  marked 
contrasts.  On  the  calibrated  composite,  the  water  areas  have  changed  very  little  In  terms  of 
actual  temperature,  but  the  terrain  surface  is  considerably  cooler.  The  factors  which  affect 
nocturnal  surface  temperatures  are  known  to  be  vegetation,  moisture,  and  the  movement  of 
cool  air  from  topographically  higher  areas  Into  lower  areas  or  depressions.  The  warm  areas 
In  the  dawn  Image  are  either  standing  water  (white)  or  the  upper  slopes  of  drumlins  or  eskers. 
The  morainal  areas  are  distinguished  from  the  outwash  plains  and  glacial  spillways  by  the  com¬ 
plexity  of  their  thermal  patterns.  The  hummocky  upland  topography  of  the  moraines  provides  a 
distinctive  thermal  contrast  to  the  relatively  cool,  level  plain  areas.  In  other  words,  nocturnal 
thermal  Imagery  can  be  expected  to  provide  information  that  Is  influenced  by  relief,  by  Pleisto¬ 
cene  topography  In  this  case. 

Thermal  Ini.agery  has  also  recently  been  used  to  locate  sinkholes  In  Karst  topography. 

In  a  study  by  A.  E.  Coker  et  al.  (22| ,  temperatures  related  to  water  stressed  and  unstressed 
vegetation  were  used  to  Indicate  areas  of  rapid  soll-molsture  depletion  (sink  holes). 

At  one  time  It  was  hoped  that  surflclal  thermal  imagery  would  provide  indications  of 
deep  subsurface  conditions,  but  recent  results  (23)  indicate  that  such  Information  Is  present 
only  when  there  are  large  thermal  gradients  within  the  overburden,  such  as  would  result  from 
underground  mine  fires. 

3.6.2.  TWO-CHANNEL  APPUCATIONS 

Section  4  of  this  report  Is  devoted  to  a  study  of  the  lithologic  properties  of  rock  outcrops 
which  provide  the  physical  basis  lor  their  discrimination  by  means  of  a  two-channel  ratlolng 
technique.  This  section  presents  examples  of  applications  of  the  ratlolng  technique  In  regions 
devoid  of  rock  outcrops.  Specifically,  the  two-channel  ratio  Is  applied  to  three  types  of  surface 
materials — surflclal  deposits,  residual  soils,  and  natural  and  crop  vegetation. 

3. 6. 2.1.  Surflclal  Deposits 

Figures  8  and  9  show  the  results  of  ratlolng  two  channels  of  thermal  data — 8.2  to  1C  9  pm 
and  9.4  to  12.1  pm — for  an  area  of  southern  Ontario.  Clearly,  agricultural  development  has 
had  a  significant  effect  upon  the  appearance  of  this  area  of  glacial  tills  and  outwash  sediments. 
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Color  Code 

Orange  <  38.2°F 
Green  =  38,2°  to  40.8of 
Red  =  40.80  to  44.70f 
White  >  44.7OF 


Old  Glacial 
Spillway 


Galt 

Moraine 


Paris 

Moraine 


Outwash 

Plain 


(a)  Video  Image 


(b)  Calibrated  Thermal  Composite  Image 


FIGURE  7.  COMPARISON  OF  STANDARD  THERMAL  VIDEO  IMAGE 
WITH  CALIBRATED  LEVEL-SLICED  THERMAL  IMAGE;  DATA  COL¬ 
LECTED  AT  0400  HOURS.  Blue  Sprlnj^  Basin,  11  May  1971. 
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(c)  Thermal  Ratio  Image 


FIGURE  8.  THERMAL  INFRARED  IMAGES  OF  PORTION  A  OF  NO.  6  HIGHWAY  BETWEEN 
HAMILTON  AND  GUELPH,  ONTARIO.  Length:  5  miles. 
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(a)  Video  Image:  8.2  to  10.9  jum  (b)  Video  Image:  9.4  to  12.1  jim  (c)  Thermal  Ratio  Image 


FIGURE  9.  THERMAL  INFRARED  IMAGES  OF  PORTION  B  OF  NO.  6  HIGHWAY  BETWEEN 

HAMILTON  AND  GUELPH,  ONTARIO 
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On  the  ratioed  thermal  images,  all  of  the  bare  soil  areas  appear  dark,  indicating  the  presence 
of  SiOg.  Areas  covered  with  vegetation  and  those  consisting  of  water,  roads,  and  man-made 
structures  are  medium  in  tone,  indicating  no  great  thermal  emittance  difference  between  the 
two  thermal  bands,  i.e.,  little  or  no  Si02.  Therefore,  ratioed  thermal  imagery  has  the  capability 
of  accurately  discriminating  cultivated  fields  (those  bare  of  vegetation).  In  Fig,  8,  the  ratioed 
image  shows  several  bare  fields  which  would  not  have  been  categorized  as  such  from  interpre¬ 
tation  of  the  thermal  video  images  alone.  Also,  in  the  upper  portion  of  the  image,  a  sand  and 
gravel  pit  is  also  dark.  Adjacent  to  the  pit,  however,  is  a  bare  field  which  has  an  anomalous 
irregular  light  tone  within  it.  The  farmer  who  cultivates  this  field  indicated  that  this  light  area 
corresponds  to  an  eroded  surface  on  a  hiU  slope.  Evidently  loss  of  topsoil  exposed  the  under¬ 
lying  parent  materials  to  view.  The  ratio  of  the  till  parent  material  is  greater  (lighter)  than 
the  soil  cover.  These  lighter  areas  may  result  from  a  greater  percentage  of  silica  in  the  up¬ 
per  soil  level  than  in  the  o.  iginal  parent  material — a  natural  result  of  the  soil  forming  process 
of  podzolization. 

The  bare  field  patterns  in  Fig.  9c  also  show  light  areas.  In  particular,  the  middle  field 
in  a  row  of  three  bare  fields  adjacent  to  one  another  appears  very  light  in  contrast  to  the  others 
on  the  ratioed  image.  In  discussing  with  the  farmer  the  differences  in  these  fields  at  the  time 
of  the  flight,  we  found  no  expi  'ation  other  than  that  the  middle  field  had  been  green  manured 
the  previous  fall.  Perhaps,  in  this  field,  accumulation  of  organic  matter  was  sufficient  to  elim¬ 
inate  the  effect  of  the  SiOg  on  the  thermal  ratio.  One  additional  field  shows  the  effect  of  the 
removal  of  topsoil,  xopsoil  from  the  northern  half  of  this  field  had  been  removed  by  a  scraper 
shortly  prior  to  the  time  of  the  flight.  The  remaining  portion  of  the  field  was  bare  of  vegetation, 
but  still  retained  its  soil  cover. 

It  was  originally  hoped  that  a  two- channel  thermal  ratio  technique  would  provide  informa¬ 
tion  on  the  relative  distribution  of  Pleistocene  deposits  on  the  basis  of  differences  in  their 
silica  content.  In  retrospect,  it  appears  that  cultivation  practices  induce  variations  which  par¬ 
tially  obscure  many  of  the  textural  differences  of  the  natural  terrain. 

The  ability  to  discriminate  areas  containing  silicates  by  means  of  a  thermal  ratio  does 
not  vary  as  a  result  of  the  time  of  day  that  the  data  were  collected.  Figure  10  compares  two 
thermal  ratio  images  for  the  same  area.  One  was  produced  from  data  collected  at  dawn 
(0400  hours),  and  the  other  from  data  collected  at  mid-day  (1200  hours).  The  patterns  indicat¬ 
ing  the  location  of  bare  soil  are  the  same  for  the  two  images.  This  lack  of  variance  in  the  ther¬ 
mal  ratio  image  with  time  of  day  is  in  sharp  contrast  to  the  results  obtained  in  the  single¬ 
channel  thermal  images.  (See  Section  3. 6.1. 2.) 
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(a)  Video  Image;  8.2-10.9  /jm.  (b)  Video  Image:  9.4-12.1  jum,  (c)  Ratioed  Image:  0490  hours 


0400  hours  0400  hours 


(d)  Video  Image:  8.2-10.9  /jm,  (e)  Video  Image:  9.2-12.1  jim,  (f)  Ratioed  Image:  1200  hours 

1200  hours  1200  hours  [ 

FiniTRE  10.  COMPARISON  OF  THERMAL  IMAGES  PRODUCED  FROM  DATA  COLLECTED  AT 

0400  HOURS  AND  AT  1200  HOURS 


I 
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3. 6. 2. 2.  Surflclal  Patterns 

Surface  soil  and  vegetation  patterns  are  often  indicative  of  subsurface  geologic  conditions. 
In  particular,  soils  derived  from  residual  parent  rocks  often  reflect  differences  related  to 
their  underlying  lithology.  Rates  of  erosion,  chemical  decomposition,  and  patterns  of  weakness 
vary  for  different  types  of  rocks,  giving  rise  to  characteristic  variations  in  topography  and  the 
appearances  of  soils.  Variations  in  color  are  often  indicative  of  microrelief,  drainage,  and 
mineral  composition.  The  purpose  of  this  task  was  to  use  ratios  of  reflective  spectral  bands 
to  enhance  surficial  patterns  on  multispectral  Imagery. 

Figure  11  shows  three  video  Images  for  an  area  of  southeastern  Pennsylvania.’*'  The 
area  is  intensively  cultivated,  except  for  the  steep,  sloping  ridge  areas  which  are  left  in  na¬ 
tural  upland  hardwood  forest.  The  pattern  in  Fig.  11  is  one  of  contrasting  light  field  areas 
and  areas  of  dark  natural  vegetation  or  crops  (data  collected  in  mid- May).  Barely  discernible 
are  tonal  patterns  within  the  light  fields  and  the  dark  vegetated  areas.  A  red- green  ratio  was 
used  for  enhancing  pattern  differences  within  the  bare  soil  areas,  and  an  infrared -green  ratio 
was  used  for  vegetative  patterns. 

Figure  12a  shows  the  results  of  ratlolng  the  0.66-  to  0.72- fun  band  and  the  0.50-  to 
0.52-Mni  band.  The  ratioed  image  shows  almost  all  of  the  vegetated  areas  as  dark,  but  within 
bare  fields  are  contrasts  of  light  and  dark  related  to  slight  differences  in  the  redness  of  the 
soils.  The  dark  mottled  areas  correspond  to  either  eroded  knobs  or  to  depressions  within 
these  fields.  Figure  12b  shows  a  ratioed  image  of  the  0.74-  to  0.85- jim  and  the  0.50-  to 
0.52-iLun  bands.  In  this  image,  all  bare  fields  are  uniformly  dark,  but  vegetation  patterns  are 
enhanced.  In  many  of  the  fields,  crop  patterns  show  the  effects  of  row  direction,  but  in  other 
cases,  soil  and  relief  variations  which  affect  the  crops  appear  as  light  or  dark  patterns.  Com¬ 
bining  the  information  contained  in  the  first  two  ratioes.  Figure  12c  shows  soil  and  vegetative 
patterns  for  the  entire  area.  In  this  ratio,  red- plus- infrared  divided  by  green,  all  contrast 
between  crops  and  bare  soil  is  lost,  but  surface  patterns  which  cross  soil  and  vegetated  areas 
can  be  discerned.  Also,  man-made  features,  such  as  roads,  buildings,  and  farm  ponds,  are 
clearly  identifiable  their  dark  tone. 

Further  applications  of  this  technique  will  attempt  to  enhance  surface- terrain  differences 
related  to  variations  in  the  mineral  composition  of  the  soil  and  perhaps  to  crop  types.  In  a 
recent  work  by  Karmanov  [24]  for  example,  a  ratio  technique  was  used  to  characterize  the 
spectral  properties  of  soils  as  indicative  of  the  content  of  organic  matter,  leaching,  podzollza- 
tion,  and  ferrallitlzatlon.  The  application  of  these  techniques  to  areas  in  which  soil  patterns 
reveal  underlying  lithology  are  thought  to  be  useful  for  purposes  of  tunnel- site  selection. 

*The  data  were  collected  for  the  Federal  Highway  Administration,  U.  S.  Department  of 
Transportation. 
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(a)  0.50-0.52  /im  (Green) 


(b)  0.66-0.72 /jm  (Red) 


(c)  0.74-0.84  ;jm  (Infrared) 

FIGURE  11.  COMPARISON  OF  THREE  VIDEO  IMAGES  FOR  AN  AREA  IN  SOUTHEASTERN 

PENNSYLVANIA 
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FIGURE  12.  COMPARISON  OF  THREE  RATIOED  IMAGES  FOR  AN  AREA  IN  SOUTHEASTERN 
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3.6.3.  MULTICHANNEL  PROCESSING  TECHNIQUES  FOR  THE  CLASSIFICATION  OF 
LAND  USES 

Recent  research  has  demonstrated  that  such  terrain  features  as  reck  types  [25],  soils 
[26],  and  urban  materials  [27]  (e.g.,  concrete,  asphalt,  lawn,  gravel,  buildings,  and  trees)  can 
be  discriminated  on  the  basis  of  their  spectral  reflectance  curves.  Multichannel  processing 
techniques  provide  for  the  development  of  diagnostic  signatures  from  multispectral  data,  allow¬ 
ing  the  distribution  of  these  features  to  be  mapped. 

Inherent  in  the  classification  of  land  uses  are  differing  hydrological  characteristics  as¬ 
sociated  with  each  class.  In  particular,  as  urban  development  (characterized  by  buildings,  hard- 
surface  roads,  and  parking  lots)  encroaches  on  rural  or  agricultural  arets,  surface  permeability 
and  evapotranspiration  are  reduced,  with  concomitant  changes  in  runoff  anl  microclimate. 

These  changes  in  hydrologic  characteristics  brought  about  changes  in  surface  condition  can 
be  measured  in  the  field  and  are  well  documented.  As  a  result,  rapid  land-use  determinations 
become  important  to  the  hydrologist  for  assessing  watershed  conditions  and  for  monitoring 
changes  over  time. 

Multispectral  data  of  a  1650-hectare  (6.4  sq  mile)  area  in  the  Oakville  Creek  Basin  were 
processed  on  a  special  purpose  analog  computer.  This  area  encompassed  the  town  of  Milton, 
Ontario,  and  neighboring  agricultural  land.  (Color  aerial  photographs  of  the  area  are  shown  in 
Fig.  13a  and  b.)  Five  primarily  rural  land-use  classes,  and  three  urban  land  uses  were  identi¬ 
fied.  These  are:  (1)  crop  vegetation,  (2)  marsh  vegetation,  (3)  forest,  (4)  bare  fields,  (5)  stand¬ 
ing  water,  (6)  asphalt  r?^ds  and  lots,  (7)  gravel  areas,  and  (8)  concrete  roads  and  buildings 
(roofs).  The  objective  oi  this  task  was  to  discriminate  areas  of  each  of  the  above  classes  and 
to  compute  total  areas  represented  by  each. 

Automatic  spectrum- matching  techniques  were  used.  An  analog  computer  identified  and 
printed  out  the  locations  of  terrain  classes  for  which  samples  were  available.  The  recognition 
criteria  were  based  on  distinctive  characteristics  recorded  for  each  sample  over  the  0.40-  to 
0.85- pm  spectral  range.  From  a  total  of  12  spectral  channels  synchronously  recorded  by  the 
multispectral  scanner,  a  subset  of  6  channels  provided  the  means  and  standard  deviations  for 
the  spectral  signatures.  These  channels  are  listed  in  Table  2. 

A  likelihood-ratio  decision  rule  was  Implemented  by  the  computer  to  determine  whether 
scene  elements  were  spectrally  similar  enough  to  be  recognized  as  a  given  terrain  class.  K 
the  spectral  characteristics  of  one  or  more  spectral  bands  for  a  land-use  class  of  interest  are 
different  from  the  characteristics  of  all  other  scene  elements,  that  class  is  automatically  dis¬ 
criminated  by  the  computer  and  recognized  on  play- back  imagery.  In  actual  practice,  multi¬ 
spectral  processing  techniques  are  used  to  establish  an  average  signature  for  a  class  of  objects, 
represented  by  a  number  of  contiguous  resolution  elements  known  as  a  training  set.  These 
average  signatures  have  a  mean  value,  variance,  and  covariance. 
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(a)  Color  Aerial  Photograph  (b)  Color  Infrared  Aerial  Photograph 


(c)  Video  Image:  0.43-0.48  pm  (d)  Video  Image:  0.58-0.65  pm  (e)  Video  image:  0.72-0.92  pm 

FIGURE  13.  AERIAL  PHOTOGRAPHS  AND  VIDEO  IMAGES  OF  OAKVILLE  CREEK 

BASIN  AREA 


36 


WILLOW  RUN  LABORATORIES 


TABLE  2.  SIX  SPECTRAL  CHANNELS  USED  FOR 
RECOGNITION  OF  LAND- USE  FEATURES 

Channel  %)ectral  Range  Color 

(4m) 

1  0.43-0.48  Violet 

2  0.46-0.49  Blue 

3  0.54-0.60  Yellow-Orange 

4  0.58-0.65  Red 

5  0.66-0.76  Dark  Red 

6  0.72-0.92  Reflective  Infrared 
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The  ability  to  establish  signatures  and  to  discriminate  terrain  classes  automatically  Is 
limited  by  the  res<jlutlon  of  the  system.  In  an  optical  mechanical  scanner  system,  the  ground 
patch  viewed  at  any  given  Ins'ant  by  the  telescope  Is  a  resolution  element  In  size.  The  resolu* 
tlon  of  the  Michigan  scanners  Is  approximately  2  mrad.  For  this  Investigation,  a  resolution 
element  represents  an  area  approximately  24  ft  (7.2  m)  on  a  side. 

Figures  13c,  d,  and  e  show  video  images  In  three  spectral  bands  for  the  Milton  area. 

(The  skew  In  the  Imagery  Is  caused  by  cross*wind  Induced  movement  of  the  data* collection 
aircraft.)  The  features  associated  with  different  land  uses  have  considerably  different  degrees 
of  contrast  In  the  three  Images.  This  consistent  change  In  contrast  (or  relative  radiance)  with 
wavelength  allows  the  Identification  and  mapping  of  the  various  land  use  classes  with  multi- 
spectral  data. 

Figure  14  contains  six  examples  from  twelve  different  recognition  Images  produced  for 
this  Milton  area.  All  white  areas  In  an  Image  are  those  designated  by  the  computer  as  being 
spectrally  similar  to  an  example  of  the  given  land-use  class.  In  each  class,  the  areas  of  rec¬ 
ognition  are  mutually  exclusive,  and  no  area  Is  identified  as  more  than  one  class.  The  Images 
for  the  eight  classes  are: 

(1)  Standing  Water.  Two  irregular  reservoir  areas  were  located  within  the  town  of 
Milton  and  one  rectangular  reservoir  near  the  highway  (401). 

(2)  Woodlots.  Five  areas  were  identified  north  and  west  of  Milton,  with  scattered  spots 
of  recognition  adjacent  to  portions  of  Oakville  Creek. 

(3)  Marsh  Vegetation.  This  class  was  principally  located  near  Oakville  Credc  before  It 
enters  the  reservoir  in  Milton. 

(4)  Green  Vegetation.  Recognition  is  comprised  of  three  Images  which  display  differences 
in  the  densities  of  crops  and  urban  lawn  vegetation. 

(5)  Gravel.  A  single  image  shows  numerous  areas  of  gravel  associated  with  parking 
lots  and  the  shoulders  of  roads  and  highways. 

(6)  Bare  Soil.  Two  Images  show  areas  of  plowed  fields  and  numerous  roads  within  the 
town  of  Milton.  One  Image  shows  portions  of  fields  mapped  as  moderately  drained, 
Chinguacousy  clay  loam,  and  the  other  represents  poorly  drained  areas  commonly 
mapped  as  Jeddo  clay  loam. 

(7)  Concrete.  Numerous  building  roofs  and  lots  and  a  single  lane  of  the  main  highway 
were  recognized  as  a  material  thought  to  be  concrete. 

(8)  Asphalt.  Parking  lot  areas  and  roads  not  recognized  as  dirt  or  concrete  were  dis¬ 
criminated  as  a  material  thought  to  be  asphalt. 

The  area  and  percentage  of  total  area  of  each  of  these  classes  are  listed  In  Table  3  below. 
Of  the  approxlm.it  >ly  1650  hectares  which  comprise  this  area,  1340  hectares  (81.2%)  are 
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(n)  Dare  Soil  (b)  Crop  and  Urban  Vegolatlon  (c)  Woodlols 


(d)  Concrete  (e)  Asphalt  (f)  Gravel 

FIGURE  14.  RECOGNITION  IMAGF.S  OF  SIX  DIFFERENT  LAND-USE  CLASSES  IN  THE 

OAKVILLE  CREEK  DASIN  AREA 
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TABLE  3.  AREA  AND  PERCENTAGE  OF  TOTAL  AREA  FOR  LAND- USE 
CLASSES  IN  THE  REQON  NEAR  MILTON,  ONTARIO 


Land-Use  Class 

Number  of  Points 

Percentage  of  Total 

Area 

(hectares) 

Standing  Water 

33,300 

1.8 

29.4 

Woodlots 

33,200 

1.8 

29.4 

Marsh 

12,300 

0.7 

10.9 

Green  Vegetation 

715,300 

38.3 

632.1 

Gravel 

40,200 

2.2 

35.5 

Bare  Soil 

387,500 

20.8 

342.5 

Concrete 

11,500 

0.6 

10.2 

Asphalt 

283,050 

15.2 

250.1 

Not  Classified 

350,000 

18.6 

309.4 

Total  Area 

1,866,350 

100.0 

1650.0 

40 


WILLOW  RUN  LABORATORIES 


classified.  Scene  elements  not  classified  are  those  which,  because  of  their  unique  spectral 
properties,  do  not  fall  Into  any  of  the  eight  classes.  These  •uiiclasslfled  scene  elements  might 
Include  such  factors  as  different  colored  roofs,  minor  variants  In  soil  and  vegetation,  and  differ¬ 
ing  surface  matei.al  within  the  urban  area.  Initial  inspection  of  the  Imagery  suggests  that 
areas  of  vegetation  comprise  most  of  the  unclassified  area. 

In  addition  to  some  areas  not  recognized.  It  is  probable  that  a  portion  of  the  areas  are 
classified  Incorrectly.  In  other  words.  It  la  llkfiy  that  the  spectral  characteristics  of  some 
scene  points  will  be  similar  to  terrain  claooco  t-'  which  they  do  not  belong.  Some  of  the  asphalt 
recognition  may  fall  Into  this  class.  Further  analysis  of  the  accuracy  of  the  recognition  results 
Is  proceeding. 


4 

COMPOSITIONAL  REMOTE  SENSING  OF  ROCK  OUTCROPS 
BY  IMAGE- RATIOING  TECHNIQUES 

According  to  K.  Sz^chy  [1| ,  the  first  step  In  geologic  exploration  prior  to  tunnel  con¬ 
struction  Is  a  survey  of  the  morphology,  petrography,  stratigraphy,  and  hydrology  of  the  en¬ 
vironment,  Including  thorough  field  reconnaissance  and  surface  explorations.  Among  the  Im¬ 
portant  data  needed  for  decisions  on  tunnel- site  selection  are  the  origin  and  actual  condition 
of  exposed  rocks,  geochemical  correlations  between  bedrock  and  topsoil,  stream  drainage,  and 
pH  values  of  springs  and  streams.  The  ratio  techniques  discussed  In  this  section  can  be  useful 
In  the  collcctlcm  of  these  and  other  pertinent  data,  once  the  techniques  are  understood  quantita¬ 
tively.  The  purpose  of  this  section  Is  to  explain  the  ratlolng  techniques,  to  report  what  has 
been  done  to  Improve  them,  and  to  suggest  ways  in  which  they  can  be  useful  for  tunnel- site 
s<  Ir.'k  tlon. 

The  thermal  Infrared  ratio  method,  originally  described  In  a  report  to  NASA  [28]  and  In 
the  Journal  of  Geophysical  Research  [29],  was  discussed  In  Section  3  of  the  semiannual  report 
for  this  contract  [16].  The  two  references  mentioned  above  show  that  the  ratio  magnitude  Is 
at  least  crudely  correlated  with  the  gross  percentage  of  SiOj  In  the  target;  therefore,  this  tech¬ 
nique  has  been  useful  for  discrimination  among  widely  different  silicate  rock  types  [31].  For 
completeness,  the  discussion  from  the  semiannual  report  Is  repeated  here  as  Section  4.1. 

4.1.  IMAGERY  OPTIMIZATION  ANALYSIS 

Geologic  remote  sensing  has  at  Its  disposal  two  new  techniques  of  Image  ratlolng  which 
are  applicable  (1)  to  the  discrimination  between  silicate  and  nonslllcate  geologic  formations 
and  (2]|  to  the  determination  of  the  amount  of  ferric  oxide  In  given  geologic  features.  The  first 
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technique  allows  the  discrimination  of  silicate  rocks  on  the  basis  of  a  ratio  of  the  radiances  in 
two  infrared  (IR)  channels.  Utilizing  25  rock  samples  selected  from  R.  J.  P.  Lyon  [30],  Fig.  15, 
Shows  the  ratio  between  channel  1  (8.2  to  10.9  iim)  and  channel  2  (9.4.  to  12.1  /im)  and  the  amount 
of  SiOg  in  percentage  of  weight,  with  the  ratio  symbolized  by  Rj  g-  The  crude  correlation  be¬ 
tween  Rj  2  and  the  percentage  of  SiOg  is  evidenced  by  the  straight  curve  on  the  graph.  The 
wavelengths  of  the  silicon- oxygen  reststrahlen  bands  are  monitored  by  R,  „.  The  more  fels  c 
rocks  (those  with  greater  SiOg  content),  generally  speaking,  have  reststrahlen  bands  (relatively 
large  departures  from  unit  emissivity)  at  shorter  wavelengths  than  do  the  mafic  rocks  (those 
with  lesser  SiOg  content),  and  therefore  Rj^  2  increases  with  decreasing  Si02. 

The  second  image- ratioing  technique  involves  a  ratio  of  a  channel  in  the  visible  green  re¬ 
gion  of  the  spectrum  (channel  5,  0.50  to  0.52  /im)  and  one  in  the  reflective  IR  (channel  7,  0.74 
to  0.85  /im).  This  ratio  can  show  the  rise  in  reflectance  of  ferric  oxides  from  the  green  to  the 

reflective  IR  wavelength.*?.  This  spectral  phenor.'enon  is  caused  by  electronic  transitions  of 

3+ 

the  Fe  ion  in  the  ferric  oxides  hematite  and  limonite. 

The  first  technique,  utilizing  the  ratio  of  two  IR  channels,  is  demonstrated  in  Figs.  16-19, 
which  display  infrared  Images  oi  specially  processed  scanner  data  gathered  by  The  University 
of  Michigan  aircraft  at  an  altitude  of  3000  ft  over  a  sand  quarry  at  Mill  Creek,  (Mclahoma,  and 
over  an  area  near  Pisgah  Crater  in  Southern  California  (also  at  3000-ft  altitude).  The  Mill 
Creek  data  were  gathered  on  June  25,  1970,  at  1000  hours  local  time,  while  the  data  from  Pisgah 
Crater  were  collected  at  0800  hours  local  time  on  October  30,  1970.  Each  of  the  figures  includes 
analog  infrared  images  of  two  single- wavelength  channels  of  data,  plus  an  image  of  the  ratio  of 
the  radiances  in  the  two  channels.  The  ratio  image  reflects  changes  in  the  chemical  composition 
of  the  rock  targets,  primarily  changes  in  their  Si02  content. 

Figure  16  shows  images  of  the  Mill  Creek  sand  quarry.  In  the  ratio  image  (Fig.  16c),  all 
features  that  appear  dark  (except  for  wavy  noise  lines)  are  either  exposed  quartz- sand  or  sand¬ 
stone.  The  outcroppings  around  the  lake  in  the  lower  part  of  the  images  are  sandstone.  Geolo¬ 
gists  from  both  the  United  States  Geological  Survey  and  The  University  of  Michigan  have  found 
almost  perfect  correlation  between  the  dark  regions  indicated  on  the  ratio  image  and  quartz- 
sand  or  sandstone.  This  type  of  discrimination  should  prove  useful  in  explorations  for  construc¬ 
tion  materials,  such  as  sand  and  gravel. 

Figure  17  is  the  first  third  of  a  north- to- south  flight  over  Pisgah  Crater,  located  in  the 
right-hand  portion  of  each  of  the  three  images.  The  right  half  of  these  images  consists  primar¬ 
ily  of  basaltic  lava  flows  which  erupted  from  Pisgah  Crater  in  the  late  Pleistocene  or  Recent 
epochs.  These  relatively  young  lava  flows  (approximately  10,000  to  100,000  years  old)  erupted 
in  three  phases.  Although  similar  in  chemical  composition,  they  have  quite  different  surface 
textures,  varying  from  ropy  to  blocky.  The  older  fanglomerate  in  the  left  half  of  the  images 
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©  For  T  =  3000k 
For  T  =  3200k 
For  T  =  280OK 

T  =  Rock  Temperature 
Liji  =  Total  Radiance  in  i-th 
Band  Received  from 
Rock  by  Detector  at  an 
Altitude  of  3000  ft  on  a 
Clear  Day. 

Line  Fit  by  Linear  Regression. 


FIGURE  15.  RATIO  BETWEEN  RADIANCES  IN  CHANNEL  1  AND  CHANNEL  2 
VERSUS  THE  PERCENTAGE  OF  WEIGHT  OF  SiOg  IN  25  ROCK  SAMPLES. 

Samples  taken  from  Ref,  [30] . 
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(b)  Channel  2:  9.4-12,1  jum 


Sandstone 

Outcroppings 


(c)  Ratio  Image  of  Channel  1  over  2 


FIGURE  16.  DISCRIMINATION  OF  ACIDIC  SILICATES 
NEAR  MILL  CREEK,  OKLAHOMA  (SAND  QUARRY) 


(a)  Channel  1:  8.2-10.9/1 
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(a)  Channel  1:  8.2-10.9  am 


(b)  Channel  2:  9.4-12.1  um 


Plsgah  Lava  (Basalt) 


I  Older  Fanglomerate^^J 
(c)  Ratio  Image  of  Channel  1  over  2 


FIGURE  17.  ANALOG  INFRARED  IMAGES  OF  FLIGHT  LINE  1.  SECTION  A,  NEAR  PISGAH 

CRATER,  CALIFORNIA 
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(b)  Channel  2;  9.4-12.1  uni 


(c)  Ratio  Image  of  Channel  1  over  2 


FIGURE  18.  ANALOG  INFRARED  IMAGES  OF  FLIGHT  LINE  2,  SECTION  A  NEAR  PISCAH 

CRATER,  CALIFORNIA 
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(a)  Channel  1:  8.2-10.9  |j 


(b)  Channel  2:  9.4-12.1  uni 


(c)  Ratio  Image  of  Channel  1  over  2 


FIGURE  19.  ANALOG  INFRARED  IMAGES  OF  FLIGHT  LINE  2,  SECTION  R.  NEAR  PISGAII 

CRATER,  CALIFORNIA 


I 
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consists  primarily  of  rock  fragments  from  both  the  basaltic  lava  and  the  felsic  mountains  on 
both  sides  of  this  region.  North  of  the  highway,  parts  of  the  second  eruptive  phase  of  the  Pisgah 
lava  are  partially  covered  by  alluvium. 

With  the  exception  of  Pisgah  Crater,  topographically  the  highest  feature  in  the  images, 
the  single- channel  images  show  very  little  temperature  variation  across  the  scene.  The  ratio 
image  (Fig.  17c),  however,  shows  considerable  variation  in  the  wavelengths  of  the  reststrahlen 
bands  over  the  same  region.  Lighter  and  darker  areas  on  the  ratio  image  indicate  reststrahlen 
bands  at  longer  and  shorter  wavelengths,  respectively.  The  mafic  rocks,  which  contain  lower 
amounts  of  SiOj,  should  appear  bright  (indicating  a  high  ratio),  and  the  felsic  rocks,  containing 
a  higher  SlOj  content,  should  appear  dark  (evidence  of  a  low  ratio).  In  the  ratio  image  (Fig.  17c), 
the  basaltic  lava  appears  bright  in  contrast  to  the  darker  tones  of  the  fanglomerate  and  gravel. 
The  basaltic  peninsula  ne.or  the  center  of  Fig.  17c  can  easily  be  discriminated  from  the  surround¬ 
ing  fanglomerate  and  gravel,  whereas  the  single- channel  images  (Fig.  17a  and  17b)  show  very 
low  thermal  contrast  for  this  feature.  In  addition,  the  wind-blown  sand  extending  over  the  lava 
just  south  of  the  peninsula  shows  up  readily  in  the  ratio  image,  but  is  not  distinguishable  in  the 
single- cliannel  Images.  In  the  ratio  image,  the  three  bright  spots  on  Pisgah  Crater  are  the  only 
experimental  artifacts;  these  resulted  from  the  much  higher  temperature  on  the  elevated  sun¬ 
ward  slopes,  which  caused  clipping  of  the  signal  in  one  of  the  channels.  Patches  of  exposed 
basaltic  lava  north  of  the  highway  r  jt  al.so  evident  in  Fig.  17c. 

Figure  ISshows  the  westernmost  halfofa  west-to-east  flight  line  over  a  region  south-south- 
east  of  Pisgah  Cr.iter.  From  ieft  to  right,  the  warm  region  (the  light  area  in  the  single-channel 
images)  in  the  extreme  western  part  is  a  dacitic  mountain;  the  adjacent  color  region  is  f.anglom- 
erate  and  gravel  (grading  eastwardly  to  smaller  fragments  and  sand);  and  the  broad,  warm  re¬ 
gion  is  the  basaltic  &inshinc  lava  flow.  Just  east  of  the  lava  flow  is  alluvium,  followed  by  playa 
deposits  (primarily  clay  and  carbonates)  and  the  southern  tip  of  the  basaltic  Pisgah  lava  flow. 

Once  again,  the  ratio  image  (Fig.  ISc)  shows  emissivity  variaMons  indicative  of  rock  type.  The 
felsic  mountains  and  fanglomerate  appear  darker  than  the  mafic  lava,  and  the  playa  material 
is  contrasted  sharply  against  the  alluvium.  Because  the  mountains  are  warmer  than  other  parts 
of  the  scene,  computational  temperature  corrections  should  heighten  the  contrast  Iwtween  dacite 
and  basal*.  t)>  lowering  the  ratio  Illustrated  in  the  ratio  image.  In  other  words.  Fig.  18c  shows 
a  higher  ratio  for  this  section  of  the  image  than  would  be  expected  on  the  basis  of  the  wave¬ 
lengths  of  the  reststrahlen  bands  alone.  In  contrast,  the  patchy  appearance  of  the  playa  is  caused 
solely  bv  emissivity  variations,  since  the  temperature  is  uniform  across  the  floor  of  Lavic  Lake 
of  which  this  section  is  a  part. 

Figure  19,  a  continuation  of  Fig.  18,  shows  the  western  half  of  the  flight  line  across  Lavic 
Lake.  From  west  to  east,  the  first  .apparent  fe.ature  is  the  Pisgah  basaltic  lava  with  playa  ma¬ 
terial  below  it,  then  alluvium,  a  basaltic  mountain,  .and  an  Intermittent  stream  bed,  now  dry,  which 
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has  eroded  down  through  the  surrounding  mountains  of  basalt.  Once  again,  on  the  ratio  Image 
(Fig.  19c),  the  southern  end  of  the  Plsgah  lava  flow  appears  similar  to  the  playa  material,  and  the 
playa-alluvlum  boundary  Is  sharply  demarcated.  On  the  northern  side  of  the  basaltic  mountain, 
the  small,  dark  regions  near  the  center  of  the  Image  have  been  confirmed  by  recent  field  trip  to 
be  felslc  outcroppings.  The  two  most  prominent  dark  areas  In  the  stream  bed  are  regions  of 
felslc  volcanic  tuff,  which  contain  minor  amounts  of  malachite,  a  copper  carbonate.  The  basaltic 
mountain  In  the  lower  right  corner  of  Fig.  19c  appears  bright,  as  did  the  basaltic  lava  flow  In 
Fig.  18c.  As  In  the  previous  ratio  Images,  In  Fig.  19c  the  felslc  rocks  and  rock  fragments  com¬ 
prising  the  alluvium  appear  to  have  lower  ratios  (and  therefore  appear  darker  In  the  ratio  Image) 
than  do  the  mafic  rocks. 

These  Images  show  th.'it  by  ratlolng  Images  In  the  two  IR  channels  8.2  to  10.9  pm  and  9.4  to 
12.1  pm,  one  can  enhance  emlsslvlty  variations  In  the  presence  of  temperature  gradients  on  the 
order  of  10°C  or  IS^C  without  computational  temperature  corrections.  Even  In  shadowed  re¬ 
gions,  compositional  Information  Is  still  present  In  the  ratio  Images.  Although  temperature- 
corrected  ratio  maps  are  available  In  digital  form,  they  do  not  show  the  contrasts  In  geologic 
formations  as  dramatically  as  the  analog  Infrared  ratio  Images  do. 

Figures  20  and  21  Illustrate  the  second  ratlolng  method— l.e.,  the  ratlolng  of  channels  In 
the  visible  green  and  the  reflective  IR  spectral  regions.  Figure  20  shows  the  same  area  as  that 
of  Fig.  18  (south -southeast  of  Plsgah  Crater),  except  that  the  data  have  been  recorded  In  channel 
5  (0.50  to  0.52  pm)  and  channel  7  (0.74  to  0.85  pm).  The  ratio  of  these  two  spectral  regions  Is 
sensitive  to  the  presence  of  Iron  oxides,  with  the  darker  areas  corresponding  to  greater  contents 
of  Iron  oxide.  In  the  ratio  Image  of  Fig.  20,  the  andesite  dikes  In  the  daclte  porphyry  mountain 
are  very  clearly  outlined,  while  they  are  hardly  discernible  In  the  single-channel  maps  of  the 
area.  The  greater  amount  of  Iron  oxide  In  the  basaltic  lava  (as  contrasted  to  that  In  the  alluvium) 

Is  also  clearly  Illustrated.  Figure  21,  winch  covers  the  ime  area  as  that  of  Fig.  17,  shows  that 
the  lava  flows  around  Plsgah  Crater  vary  In  Iron  oxide  content,  a  result  In  agreement  with  our 
analy.«=  s  of  samples  recently  collected  In  the  field  from  each  of  the  three  eruptive  phases. 

These  ratlolng  techniques  have  just  recently  been  developed*  and  are  presently  being  Im¬ 
proved.  Under  the  present  contract,  new  laboratory  data  are  being  measured,  and  the  results 
of  these  measurements  should  greatly  Improve  our  ability  to  relate  the  measured  Rj  j  h)  the 
chemical  parameters  of  geologic  targets.  Either  the  relationship  between  SlOj  and  Rj  ^  1" 

Fig.  15  will  be  Improved,  or  a  different  chemical  or  mlneraloglcal  parameter  better  related 
*^1,2  "'***  found.  In  regard  to  the  R^  ratio  technique,  several  samples  from  the  Plsgah 
Crater  test  site  are  presently  undergoing  laboratory  analyses  (both  spectral  and  chemical)  In 
an  effort  to  prove  quantitatively  the  relationship  between  spectral  reflectivity  features  and  the 


*As  of  September  1971. 
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(c)  Ratio  Image  of  Channel  5  ovc^  7 

FIGURE  20.  ANALOG  VISIBLE  AND  NEAR-INFRARED  IMAGES  OF  FLIGHT  LINE  2.  SECTION 

A.  NEAR  PISGAII  CRATER.  CALIFORNIA 
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(b)  Channel  7:  0.74-0.85  /im 


(r)  fbitio  Imace  of  Ch.mnel  5  over  7 


FIGUnE  21.  ANALOG  VISIBLE  AND  NEAR-INFRAHED  IMAGES  OF  FLIGHT  LINE  1.  SECTION 

A.  NEAR  PISGAII  CILXTER.  CALIFORNIA 
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Iron  oxide  content  of  exposed  rock  surfaces.  It  seems  likely  that  the  new  developments  In 
Imagc-ratlolnR  techniques  will  lead  to  Improved  methods  for  rapid  geologic  mapping  and  litho¬ 
logic  differentiation. 

4.2.  IMPROVEMENT  OF  THE  THERMAL  INFRARED  RATIO  TECHNIQUE 

As  indicated  In  the  previous  section,  the  thermal  Infrared  ratio  technique  has  been  limited 
by  lack  of  knowledge  about  the  precise  relationships  between  the  Infrared  spectral  features  of 
rocks  and  their  chemical  and  mlneralogieal  composition.  For  this  reason,  a  laboratory  program 
was  Initiated  to  provide  more  Information  about  this  relationship.  Twenty- six  Igneous  rock 
samples  from  various  ports  of  the  world  were  selected  as  a  test  set.  Table  4  lists  these  rocks 
and  gives  the  chemical  composition  of  those  for  which  analyses  have  already  been  completed. 

The  chemical  compositions  arc  currently  being  determined  for  the  remaining  specimens  and 
will  be  completed  by  July  1972.  The  U.  S.  Geological  Survey,  through  Dr.  Larry  Rowan  of  the 
Branch  of  Regional  Geophysics  In  Denver,  Colorado,  supplied  all  of  the  rocks  except  those 

from  Plsgah  Crater,  California,  and  also  Is  providing  the  chemical  analyses  at  no  cost  to  the 
Bureau  of  Mines. 

As  can  be  Inferred  from  Table  4,  the  silicate  rocks  In  this  group  vary  from  highly  felslc 
to  highly  mafic.  On  a  parabolic  renectometer  (xhown  In  Fig.  22)  owned  by  the  Martin-Marietta 
Corporation,  Denver  Division,  the  thermal  Infrared  spectra  In  the  approximate  wavelength  region 
from  7.5  to  14.0  pm  were  measured  tor  an  exposed*  surface  of  each  of  these  rocks.  Although 
the  renectometer  Ulumlnatcs  the  whole  hemisphere  of  the  sample,  the  rock  Is  viewed  In  a  smaU 
solid  angle  of  observation.  This  experimental  condition  Is  slmUar  to  that  of  a  total  Integrating 
sphere.  In  which  the  sample  Is  lUumlnated  within  a  smaH  solid  angle  of  Incidence  .and  the  ob¬ 
served  rcHccted  radiation  Is  collected  over  the  whole  hemisphere.  However,  from  a  design 
standpoint,  the  parabolic  renectometer  Is  superior  to  the  Integrating  sphere  In  the  thermal 
Infrared  wavelength  region.  First,  It  is  simpler  to  buUd  and  maintain  a  parabolic  aluminum 
mirror  than  an  Integrating  sphere  which  has  smooth  surfaces  compared  to  wavelengths  on  the 
order  of  10  pm.  Second,  the  parabolic  renectometer  offers  an  opportunity  to  observe  a  given 
s.ample  In  both  emission  or  renectlon,  with  almost  the  same  optical  train,  whereas  this  Is  not 
feasiUe  with  an  Integrating  sphere. 

For  renectlon  measurements,  a  heated  blackbody  cavity  at  600OC  acted  as  the  source  of 
infrared  lUumlnatlon.  The  source  radiation  was  chopped  to  minimize  the  contribution  of  radi¬ 
ation  cml»»«d  by  the  sample,  which  was  maintained  at  a  temperature  of  approxin;ately  300OK  by 
means  of  Ice  water  circulated  through  the  sample  holder  and  cooled  nitrogen  vapor  blown  across 

envlronr^'nt  surface  is  defined  as  one  exposed  to  air  In  the  sample’s  na.ural 
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FIGURE  22.  PARABOLIC  REFLECTOMETER  OPTICAL  SCHEMATIC 
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the  sample's  surface.  The  observation  angle  in  all  cases  was  15°  to  the  macroscopic  surface 
normal  to  the  sample.  A  portion  of  the  radiation  from  the  sample  was  collected  by  the  transfer 
optics  through  a  hole  In  the  second  parabolic  mirror  and  focused  (with  an  approximate  11°  field 
of  view)  by  mirrors  M3,  M4,  and  M5  onto  the  entrance  silt  of  a  Perkln-Elmer  Model  98  mono¬ 
chromator.  Dispersion  In  the  monochromator  was  accomplished  with  a  CsBr  prism,  which  had 
a  spectral  slit  width  that  varied  from  AX  «  o.51  /im  at  X  -  8  /un  to  AX  =  0.29  /im  at  X  =  14  /un. 
The  physical  slit  width  was  held  at  a  constant  0.5  mm  for  all  samples,  with  the  physical  slit 
length  12.0  mm.  The  signal  detector  was  a  high  sensitivity  Reeder  radiation  thermocouple, 
with  a  D*  (detectivity)  of  approximately  10®.  A  gold  specular  reflectance  standard  used  as  a 
reference  varied  In  reflectance  from  0.G921  at  X  ■  8  ;im  to  0.9930  at  X  =  14  /rni. 

The  resulting  Infrared  spectra  of  the  -6  rock  samples  are  shown  in  Fig.  23.  For  each 
sample,  the  spectral  emittance  has  been  plotted,  as  calculated  from  Klrchoff's  law 

^  ^  ■  ^X  (1) 


where  is  the  measured  directional  spectral  reflectance. 


The  connection  between  these  laboratory  spectra  and  the  infrared  ratios  measured  by  an 
airborne  scanner  wlU  be  shown  In  the  following  two  equations.  The  total  spectral  radiance  re¬ 
ceived  by  an  airborne  detector  from  the  direction  of  a  given  rock  target  was  calculated  from 
the  equation 


Where  -  tolal  spectral  radiance  of  wavelength  X  arriving  at  the  detector 

^X  "  ^  ■  ^X  “  spectral  emittance  taken  from  the  data  above  and  Eq.  (1) 

“  spectral  radiance  of  a  blackbody  at  same  temperature,  T,  ac  lock 

^sk)  “  spectral  radiance  from  sky  incident  on  diffuse  rock  surface 

T.  ■  spectr.ll  atmospheric  transmissivity 
X 

Ly  ■  spectral  radiance  from  emission  and  scattering  into  the  beam  by  the  at¬ 
mosphere  in  the  path  between  rock  and  detector 

The  atmospheric  parai.irlers  calculated  from  a  clear,  dry  summer 

X  A  X 

(visibility  of  20  miles)  atmospheric  model  developed  at  The  University  of  Michigan  [32j.  Typi¬ 
cally,  the  atmospheric  parameters  were  calculated  for  geological  targets  at  a  base  altitude  of 
0.76  km  (3000  ft)  and  at  an  aircraft  altitude  of  1  km  above  the  targets. 
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(a)  (b) 


(c) 

FIGURE  23.  INFRARED  SPECTRA  OF  IGNEOUS  SILICATE  ROCKS 


EMITTANCE 


2B5  Dacite 


(f) 

FIGURE  23.  INFRARED  SPECTRA  OF  IGNEOUS  SILICATE  ROCKS  (Continued) 
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For  two  selected  bands  and  an  assumed  rock  temperature,  Eq.  (2)  can  be  integrated  over 
the  respective  spectral  regions,  and  a  ratio  of  radiances  in  the  two  bands  for  each  of  the  26 
samples  can  be  calculated.  The  ratio  is  then  given  by 


“1,2™  ■  = 


(3) 


where  and  are  the  limits  of  band  1;  and  X^  are  the  limits  of  band  2;  and  fj(X)  and 
f2{X)  are  the  exact  spectral  responses  of  the  detector  in  the  two  respective  bands.  In  order 
that  the  theoretical  results  could  be  compared  with  an  existing  set  of  field  data,  the  two  spec¬ 
tral  channels  first  to  be  considered  corresponded  to  the  two  wavelength  regions  covered  by 
Willow  Run  Laboratories'  two-element  Honeywell  Hg;Cd:Te  detector;  i.e.,  AX^  is  8.2  to  10.9 
and  AXj  is  9.4  to  12.1  fun  at  the  10%  spectral  responsivity  points.  Actual  filter  functions  of 
these  two  channels  were  used.  The  samples  were  all  assumed  to  be  at  T  =  SOO^K.  With  the 
atmospheric  parameters  mentioned  above,  the  radiances  calculated  from  Eq.  (2)  were  sub¬ 
stituted  into  Eq.  (3)  to  yield  a  ratio  for  each  sample,  as  measured  by  this  specific  detector  at 
an  above-ground  altitude  of  1  km.  The  resulting  Rj  2(300OK)  ratios  are  given  in  column  1  of 
Table  5. 


In  column  2  of  Table  5  are  shown 


ratios  Rj  2(300°K)  for  no  atmosphere  =  1 


Lsky  =  column  3  are  ratios,  which  will  be  called 

Rf  jj(300°K),  for  two  hypothetical  channels  w^th  square-filter  functions  between  8.2  and  9.2  imi 
(channel  I)  and  10.3  and  11.3  (im  (channel  II);  the  same  experimental  conditions  as  were  used 
to  calculate  Rj  2(300°K)  were  assumed.  Finally,  in  column  4  are  the  ratios  rJ*  jj(300°K)  for 
these  same  two  hjT)othetical  channels  without  atmosphere.  Based  on  previous  studies  [28], 
channels  I  and  II  are  nearly  optimum  for  discrimination  on  the  basis  of  percentage  of  Si02. 


The  end  product  of  the  analysis  of  these  data  will  include  plots  of  the  ratios  listed  in 
Table  5  versus  percentage  of  Si02,  or  possibly  versus  some  better  chemical  or  mineralogical 
parameter(s).  Although  the  analysis  cannot  be  completed  until  the  next  contract  year,  three 
experimental  problems  merit  discussion,  and  three  useful  results  can  be  derived  from  the  data 
in  their  present  state.  The  three  experimental  problems  involve  the  small  spot  size  of  the 
reflectometer,  the  spectral  effects  of  surface  contaminants  (weathering  products)  on  a  few 
samples,  and,  for  some  samples,  the  lack  of  agreement  between  spectral  reflectance  and  direct 
spectral  emittance  measurements. 
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TABLE  5.  CALCULATED  THERMAL  INFRARED  RATIOS  BY  A  DETECTOR  AT  AN 
ABOVE-GROUND  ALTITUDE  OF  1  km  ON  A  CLEAR,  DRY  SUMMER  DAY 


Sample  Name  and  Number 

Rj  2(300°K) 

P®  2(300°K) 

Rj  jj(300°K) 

I^^(300°K) 

Granite  (A79) 

1.1761 

1.1541 

0.6756 

0.6383 

Granite  (A119) 

1.1946 

1.1812 

0.6977 

0.6868 

Granite  (A122) 

1.1987 

1.1828 

0.7030 

0.6918 

Rhyolite  (SA-49) 

1.1266 

1.0886 

0.6366 

0.5960 

Rhyolite  (Welded  Tuff) 

1.1741 

1.1555 

0.6852 

0.6731 

Syenite  (E83A) 

1.1784 

1.1540 

0.6960 

0.6846 

Trachyte  (Porphyritic) 

1.1442 

1.1053 

0.6561 

0.6211 

Nepheline  Syenite 

1.2094 

1.2058 

0.7030 

0.7016 

Granodiorite  (A117) 

1.2134 

1.2054 

0.7149 

0.7097 

Granodiorite  (A127) 

1.1723 

1.1461 

0.6834 

0.6636 

Dacite  (E12A) 

1.1724 

1.1512 

0.6715 

0.6485 

Dacite  (2A1) 

1.1436 

1.1062 

0.6566 

0.6203 

Dacite  (2B5) 

1.1574 

1.1264 

0.6737 

0.6497 

Dacite  (2B5,  Fresh  Surface) 

1.1686 

1.1457 

0.6704 

0.6476 

Diorite  (A129) 

1.1685 

1.1420 

0.6730 

0.6456 

Diorite  (E47C) 

1.1795 

1.1693 

0.6818 

0.6617 

Andesite  (E55A) 

1.1962 

1.1820 

0.7033 

0.6934 

Andesite  (2A3) 

1.1394 

1.1007 

0.6587 

0.6277 

Gabbro  (E80A) 

1.1810 

1.1567 

0.6993 

0.6870 

Basalt  (1A8) 

1.2154 

1.2087 

0.7184 

0.7151 

Basalt  (WI- 1-104) 

1.2204 

1.2180 

0.7188 

0.7190 

Basalt  (S95) 

1.2246 

1.2216 

0.7280 

0.7303 

Basalt  (5B4) 

1.2255 

1.2225 

0.7280 

0.7287 

Basalt  (1A3) 

1.2018 

1.1943 

0.7C47 

0.7026 

Anorthosite  (AN) 

1.2458 

1.2494 

0.7547 

0.7648 

Diabase 

1.2204 

1.2168 

0.7245 

0.7236 

Peridotite  (EP-1) 

1.3723 

1.4178 

0.9092 

0.9776 
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The  physical  area  observed  on  each  rock  during  the  spectral  measurements  was  approxi¬ 
mately  0.5  mm  X  12.0  mm  (the  slit  dimensions),  which  presented  a  problem  with  coarse-grained 
polymlnerallc  rocks.  The  greatest  variability  in  or  with  a  change  in  sample  position  was 
found  in  the  granite  samples,  as  expected.  Figure  24  shows  the  spectral  variations  caused  by 
taking  measurements  at  various  locations  on  three  of  the  coarsest  grained  samples:  granite, 
granodlorite,  and  peridotlte.  Whereas  such  variations  add  to  uncertainty  of  detection  via  the 
ratio  method,  this  problem  is  minor,  except  possibly  for  the  granites.  Monetary  constraints 
made  it  impossible  to  take  more  measurements  for  various  sample  positions  or  to  alter  the 
optical  train  to  enlarge  the  spot  size  without  losing  spectral  resolution. 

Surface  contaminants,  in  the  form  of  weathering  products  (both  authogenic  and  allogenic), 
constlti’led  a  second  problem,  although  one  that  was  deliberately  allowed  to  occur.  In  most 
cases,  the  surfaces  examined  on  these  samples  were  those  naturally  exposed  in  the  field,  be¬ 
cause  these  are  the  surfaces  that  a  remote  sensing  experiment  would  encounter.  Thus,  some 
idea  of  how  such  contamination  can  affect  this  type  of  experiment  can  be  gained  irom  the  more 
contaminated  samples,  such  as  the  1A3  basalt,  2B5  dacite,  and  the  nepheline  syenite.  Both  1A3 
and  2B5  are  Pisgah  Crater  area  samples  contaminated  wit.i  a  silt  that  seems  common  to  the 
playa  sediments  of  Lavlc  Lake,  which  is  within  a  2  or  3  mUe  radius  of  the  location  from  which 
each  of  these  rocks  was  collected.  The  1A3  basalt  sample  surface  was  quite  silty.  In  contrast 
to  the  cleaner  1A6  basalt  sample,  collected  within  100  yards  of  1A3,  the  1A3  sample  has  an 
emittance  minimum  at  about  9.5  /xm;  this  emittance  minimum  is  most  likely  caused  in  part  by 
the  silt,  in  agreement  with  R.  J.  P.  Lyon’s  recent  airborne  spectrometer  data  of  Lavlc  Lake 
playa  sediments  [33].  In  Fig.  23e,  curves  of  both  exposed  and  fresh  surfaces  of  the  2B5  dacite 
sample  are  shown.  They  demonstrate  that  the  9.5-Mm  emittance  minimum  in  the  weathered 
2B5  curve  is  caused  by  the  surface  contaminant,  which  appears  to  be  the  same  type  of  silt  as 
that  found  on  1A3,  i.e.,  a  light-colored  sUt  that  seems  to  be  composed  of  clay  particles. 

The  final  problem  encountered  Involved  an  experiment  that  compared  emittance  calculated 
from  reflectance  measurements  by  means  ot  Eq.  (1)  with  directly  measured  emittances.  A  45° 
folding  mirror  was  placed  in  front  of  the  source,  facing  down  into  a  container  of  liquid  nitrogen. 

A  hemisphere  of  approximately  100°K  radiation  (limited  by  the  emittance  of  the  mirrors)  was 
thus  effectively  produced  around  the  sample,  which  was  maintained  at  approximately  300°K  by 
means  of  warm  water  circulating  through  the  sample  holder,  l.ie  chopper  was  moved  to  a  new 
position,  in  front  of  the  monochromator  silt.  This  experiment  was  designed  to  determine  what 
spectral  emittance  differences  one  might  expect  when  the  sample  temperature  is  changed  from 
one  that  decreases  toward  the  interior  of  the  sample  (negative  temper  dure  gradient)  to  one 
that  increases  towaid  the  interior  (positive  temperature  gradient).  For  a  negative  temperature 
gradient,  the  emittance  is  calculated  from  reflectance  measurements,  while  for  a  positive 
gradient,  direct  emittance  measurements  are  used. 
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Figure  25  compares  the  spectral  emlttances  from  Eq.  (1)  and  from  direct  spectral  emit- 
tance  measurements  for  SA-49  rhyolite,  A79  granite,  A117  granodlorite,  ana  EP-1  perldotite. 
Other  samples  measured  in  this  way,  with  results  now  shown,  are  2A1,  2B5,  2A3,  E80A,  5B4, 
1A3,  1A8,  and  AN.  In  almost  all  of  the  direct  spectral  emittance  measurements,  a  local  emit- 
tance  minimum,  which  seems  independent  of  sample  composition,  occurs  in  the  8.3-  to  8.5- /im 
region.  This  result  is  thought  to  represent  an  experimental  artifact,  although  its  origin  has 
not  yet  been  deduced.  Besides  this  anomalous  minimum,  the  emlttances  measured  directly 
are  higher  than  those  derived  from  reflectance  measurements  and  Eq.  (1).  The  lack  of  quanti¬ 
tative  agreement  is  believed  to  result  from  the  different  temperature  gradient  in  the  two  experi¬ 
mental  cases  and  to  an  overestimate  of  the  sample's  surface  temperature  in  the  direct  emittance 
case.  The  measured  emittance  is  simply  the  detected  emitted  radiance  divided  by  a  biaitkbody 
radiance  at  the  sample’s  surface  temperature,  which  is  estimated  via  thermocouples.  Because 
of  these  olfficultles,  the  reflectance  measurements  were  chosen  as  the  principal  data  set.  The 
emlttances  from  Eq.  (1)  (plotted  in  Fig.  23)  are  for  samples  with  a  negative  temperature  gradient 
near  the  outermost  surface,  as  one  would  expect  for  rocks  in  the  day-time.  Figure  25,  with  the 
exception  of  that  artifact  in  the  8.3-  to  8.5-/im  region,  indicates  that  a  positive  temperature 
gradient,  which  would  be  expected  under  night-time  conditions,  should  not  produce  emlttances 
greatly  different  from  day-time  emlttances. 

The  three  useful  results  which  can  be  derived  from  these  data  in  their  present  state  in¬ 
volve  (1)  the  apparent  transparent  behavior  of  ferric  oxide  coatings  in  the  7-  to  14- /im  wave¬ 
length  region,  (2)  the  very  uniform  emittance  of  25  of  the  26  silicate  samples  in  the  11.75-  to 
13.75-/im  region,  and  (3)  the  excellent  agreement  between  ratios  calculated  from  the  Pisgah 
Crater  samples  in  this  data  set  and  those  actually  measured  in  scanner- produced  ratio  im¬ 
agery  from  Pisgah  Crater,  California. 

Sample  2A1  dacite  is  stained  oark  reddish-brown  by  ferric  oxide,  and  yet,  in  the  wave¬ 
length  region  from  7-14  jim,  it  has  a  distinctly  different  spectral  structure  from  the  ferric 
oxide  coated  basalt  samples  (including  1A3  and  1A8  from  Pisgah  Crater).  This  result  indicates 
that  in  this  wavelength  region  ferric  oxide  coatings  do  not  appreciably  alter  the  spectra  of  the 
underlying  rock  (for  coating  thicknesses  on  the  order  of  100  pm).  As  support  for  this  argument, 
Fig.  26  shows  an  emittance  spectrum  calculated  from  Eq.  (1)  for  a  thick  layer  of  rust  (a  milli¬ 
meter  or  more  in  depth)  on  a  piece  of  iron.  The  few  spectral  features  of  rust  are  not  observe.l  on 
the  2A1  dacite,  2 A3  andesite,  1A3  basalt,  or  1A8  basalt  samples. 

The  second  immediate  result  is  the  uniform  emittance  of  all  the  samples,  excluding 
peridotlte,  in  the  11.75-  to  13.75-/im  wavelength  region.  The  mean  emittance  for  those  samples 
in  that  spectral  region  is  0.956  ±  0.008,  while  peridotlte,  the  most  highly  mafic  sample,  has  a 
lower  mean  emittance,  about  0.895.  This  uniform  emittance  for  so  many  silicates  indicates 
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FIGURE  25.  COMPARISON  OF  DIRECT  SPEC¬ 
TRAL  EMITTANCE  MEASUREMENTS  (SOLID 
LINE)  WITH  SPECTRAL  EMITTANCE  CURVES 
DERIVED  FROM  KIRCHOFF’S  LAW  AND  SPEC¬ 
TRAL  REFLECTANCE  MEASUREMENTS 
(DASHED  LINE) 


FIGURE  24.  SPECTRAL  VARIATIONS  RESULT¬ 
ING  FROM  MEASUREMENTS  TAKEN  AT  DIFFER 
ENT  LOCATIONS  ON  THE  SAMPLES 
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FIGURE  26.  SPECTRAL  EMITTANCE  OF 
IRON  RUST 
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that  the  11.75-  to  13.75-/ini  region  is  an  excellent  one  in  which  to  determine  the  temperature 
cf  a  silicate  target.  Previously  that  assumption  was  made  to  facilitate  temperature  corrections 
of  the  Pir^ah  Crater  scanner  data  (34 j.  Now  it  is  (luantitatively  verified  to  be  a  valid  assumption. 
Such  knowledge  should  also  Improve  the  accuracy  of  absolute  temperature  measurements,  which 
arc  important  in  the  interpretation  of  thermal  variations  across  the  scene. 

Finally,  the  Rj  values  of  Plsgah  Crater  samples  from  Table  4,  when  multiplied 

by  a  temperature  correction  factor  Rgjj(300OK)  -  0.8429  [34],  agree  quite  well  with  the  mea¬ 
sured  values  from  scanner  data  collected  in  an  October  1970  flight  under  NASA  Contract  No. 
9-9784.  Figure  27  is  a  reprint  of  a  figure  from  Ref.  (35),  which  shows  a  digital  ratio  map  of 
part  of  flight  line  2,  section  A,  5  miles  southwest  of  Plsgah  Crater.  The  dacltc  in  Fig.  27 
marked  D  is  represented  by  sample  2A1,  which  has  a  temperature- corrected  ratio  of  0.964, 
correctly  represented  by  the  lightest  color  range  R  <  0.960,  where  R  »  Rg|^(300°K)Rj  ^(SOO^K)- 
Also  correctly  in  that  r.ange  Is  the  2A3  andesite  sample,  located  adjacent  to  2A1.  The  alluvium 
A  consists  prlm.arlly  of  rock  fragments  of  these  two  rock  types.  The  Plsgah  basalts  1A3  and 
1A8  have  temperature- corrected  ratios  of  1.013  and  1.025,  respectively,  which  would  make 
them  fall  Into  the  darkest  color  range  (R  >  0.992).  Accordingly,  the  basaltic  Sunshine  lava 
flows  marked  LA  (not  very  different  from  the  Plsgah  basalt  in  percentage  of  SiO,)  are  dark 
in  Fig.  27.  The  playa  material  marked  P  consists  mostly  of  clays  that  have  a  spectral  be¬ 
havior  similar  to  sample  1A3,  and  hence  patches  of  the  playa  appear  to  be  in  the  d.arkest  color 
range.  In  summary,  there  is  excellent  quantitative  agreement  between  laboratory  results  as 
shown  In  Table  4  and  field  results  (Fig.  27). 

4.3.  DESCRIPTION  AND  IMPROVEMENT  OF  THE  VISIBLE- RE  ELECTIVE  INFRARED  RATIO 

TECHNIQUE 

The  vlslble-reflectlve  Infrared  ratio  technique,  originally  described  in  the  NASA  report 
of  Ref.  [34],  is  discussed  in  Section  4.1  of  this  report,  excerpted  from  Ref.  [16].  In  short,  a 
ratio  of  nadiances  collected  In  channel  5  (0.50-0.52  ;xm)  to  those  collected  in  channel  7  (0.74- 
0.85  Mm)  was  hypothesized  as  useful  for  mapping  variations  of  Iron  oxides  on  exposed  rock 
surfaces.  That  technique  had  not,  however,  been  verified  by  laboratory  data.  The  first  step 
In  proving  its  validity  required  a  demonstration  that  the  scanner  data  represent  true  spectral 
variations  in  the  observed  rocks.  It  is  a  crucial  step,  since  good  quantitative  atmospheric 
data  have  not  previously  been  available  for  this  spectral  region. 

To  correct  this  deficiency,  spectral  reflectance  measurements  of  14  Plsgah  basalt  sam¬ 
ples  were  made  on  a  Beckmann  spectrophotometer  with  an  integrating  sphere  reflectance  at¬ 
tachment.  The  sample  was  illuminated  near  normal  incidence,  and  reflected  radiation  was 
gathered  over  the  whole  hemisphere.  The  samples  were  of  three  eruptive  phases  of  Plsgah 
basalt.  Including  sample  1A3  (phase-3  basalt)  and  1A8  (phase- 1  bas.alt).  Figure  28  shows  the 
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FIGURE  27.  TEMPERATURE -CORRECTED  RATIO  RECOGNITION  MAP  OF 
SILICATE  ROCKS  NEAR  PISGAIT  CRATER.  CALIFORNIA.  From  [34]. 
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FIGURE  28.  MEAN  REFLECTANCE  OF  THREE  LA  VIC  PHASES  OF  PISGAH  CRATER  RASALTS 
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mcMi  curves  for  the  oldest  phase- 1  lava  (averaged  over  12  curves),  pha8e-2  lava  (averaged 
over  10  curves),  and  the  youngest  phase- 3  lava  (averaged  over  12  curves).  The  Individual 
curves  were  taken  from  randomly  collected  samples  and  at  randomly  located  spots  on  the  sam¬ 
ples.  Some  samples  were  measured  at  more  than  one  spot.  Although  the  standard  deviation  is 
large  within  a  given  lava  phase,  the  mean  curves  indicate  that  the  visible- to- reflective  infrared 
ratios  should  increase  from  phase  1  to  phase  3.  The  average  reflectance  ratios  of  channel  5 
(0.50-0.52  pm)  to  channel  7  (0.74-0.85  pm)  for  lavlc  phases  1,  2,  and  3  are  0.457,  0.578,  and 
0.743,  respectively.  On  the  average,  then,  phase-1  lava  should  appear  the  darkest  of  the  three 
on  tho  Rg  ,j  ratio  image;  phase  3  should  appear  lightost;  and  phase  2  should  be  an  intermediate 
tone.  An  examination  of  Fig.  21,  which  shows  single  channel  and  R,  „  ratio  images  of  flight- 

D,  I 

line  1,  section  A,  confirms  that  the  relative  magnitude  of  the  R^  ,,  ratios  as  measured  by  the 

I 

airborne  multiL'pectral  scanner  agree  well  with  the  laboratory- measun'd  values.  A  more  com¬ 
plete  quantitative  analysis  is  not  possiUe  because  of  lack  of  accurate  atmospheric  parameters. 

Therefore,  the  Rg  ^  ratio  measured  by  a  scanner  can  be  used  to  monitor  spectral  vari¬ 
ations  of  exposed  rock  surfaces,  in  spite  of  the  fact  that  shadowing,  scan- angle  effects,  and 
target- geometry  variations  are  present  across  the  scene.  When  the  chemical  analyses  are 
completed,  an  attempt  will  be  made  to  correlate  the  Rg  ^  variations  with  Fc2®3  content  of  the 
basalt  samples. 

4.4.  RELEVANCE  OF  THESE  METHODS  TO  TUNNEL- SITE  SELECTION  AND  OTHER  BUREAU 

OF  MINES  PROBLEMS 

The  observation  of  surface  outcrops,  tracing  of  past  landslides,  and  the  determination  of 
the  number,  yield,  and  pH  value  of  existing  springs  and  water  courses  are  all  important  in 
tunnel- site  selection.  The  ratio  methods  presented  can  aid  in  each  of  these  areas.  For  in¬ 
stance,  initial  outcrop  rock-type  discrimination  can  be  done  quickly  over  large  areas  via  air¬ 
borne  scanners.  Although  only  igneous  rocks  have  been  considered  in  this  report,  these  tech¬ 
niques  should  also  bo  useful  with  sedimentary  and  metamorphlc  rocks.  Discrimination  between 
silicate  and  nonsilicate  rocks  is  relatively  easy,  as  has  been  shown  between  quartz  sandstone 
and  limestone  [35].  Image  ratiolng  is  also  a  good  method  for  detecting  outcrops,  even  if  they 
are  only  a  few  feet  in  diameter  (for  lower  altitude  flights)  and  in  small  clearings  between  the 
trees.  These  te  hniques  should  prove  superior  both  to  aerial  photography  and  field  observation 
for  the  location  of  outcrops,  except  in  very  dense  vegetation. 

These  techniques  should  prove  as  good  as  aerial  photographic  methods  for  tracing  past 
landslides  and  could  prove  to  be  even  better  if  the  rock  types  involved  are  appreciably  dif¬ 
ferent  in  Si02  content,  but  similar  in  visible  color.  The  infrared  technique  has  been  useful, 
for  instance,  in  distinguishing  volcanic  rhyolitic  ash  flown  from  less  acidic  detritus,  where 
the  two  materials  were  indistinguishable  in  aerial  photographs  [34]. 
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The  determination  ol  the  number,  yield,  and  pH  values  of  existing  springs  and  water 
courses  can  probably  also  be  assisted  via  these  techniques,  although  there  have  been  no  test 
cases  of  that  nature  as  yet.  Stream- bed  sedL’icnts  can  probably  be  correlated  with  their  parent 
materials,  if  the  whole  upstream  channel  is  flown.  Although  water  is  opaque  to  thi  rmal  in¬ 
frared  radiation,  the  water  edges  can  bo  identified.  From  these  data,  one  could  calculate  vol¬ 
ume  flow,  if  the  standard  channel  width  (including  the  whole  stream  bed)  to  channel  depth  model 
is  used,  along  with  the  measured  water  level.  Finally,  some  information  regarding  the  pH 
values  of  streams  might  be  obtained  by  m.apping  ferric  oxide  stains  on  the  banks  of  acidic 
streams  or  on  the  bottoms  of  intermittent  acidic  stream  beds.  Such  mapping  could  be  accom¬ 
plished  with  the  visible-reflective  infrared  ratio  technique.  Although  false-color  photography 
would  also  detect  the  ferric  oxide  stains,  the  sc.anner  ratio  method  is  insensitive  to  shadowing 
effects,  which  makes  it  superior  for  shadowed  stream  beds.  The  ferric  oxide  stain  mapping 
procedure  might  also  be  useful  for  monitoring  acidic  water  runoffs  from  either  surface  strip 
mines  or  abandoned  underground  coal  mines,  both  considerable  water- pollution  hazards. 

Finally,  these  techniques  should  eventually  be  useful  for  rapid  afisaylng  of  identifiable 
regions  of  surface  ore  deposits,  although  more  laboratory  measurements  similar  to  the  ones 
presented  in  this  section  would  be  needed,  so  that  correlations  could  be  made  between  observed 
ratio  values  and  bulk  content  of  the  mineral  resource  to  be  mined.  The  use  of  these  techniques 
should  be  especially  useful  when  rapid  decisions  must  be  made  by  the  Department  of  Interior 
concerning  the  designation  of  wilderness  areas. 


5 

APPLICATIONS  OF  SIDE  LOOKING  RADAR  TO  GEOLOGIC  INVESTIGATIONS 
5.1.  INTRODUCTION 

Sidelooking  airborne  radar  (SLAR)  has  imaging  capabilities  that  both  supplement  and 
complement  the  passive  remote  sensors.  Althoi^h  SLAR  effectively  collects  geologic  infor¬ 
mation  when  used  alone,  maximum  data  will  be  obtained  If  it  is  used  prior  to  ground  studies 
in  conjunction  with  other  sensors,  each  contributing  its  special  Information.  Sidelocking  radar 
has  the  same  limitation  as  multlspectral  scanning  and  aerial  photography  in  that  it  provides 
Information,  often  relatively  detailed,  about  surface  and  very  near- surface  geology,  but  no 
definite  Information  about  geologic  features  more  than  a  few  inches  beneath  the  surface,  since 
the  wavelengths  presently  being  used  (mainly  K-,  X-,  and  L-bands)  have  very  limited  penetra¬ 
tion  into  geologic  materials.  Radar  penetration  of  geologic  materials  increases  as  the  wave¬ 
lengths  get  longer.  The  use  of  much  longer  radar  wavelengths  (several  to  tens  of  meters) 
should  provide  useful  information,  depending  upon  the  terrain  conditions,  about  deeper  geologic 
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features  (such  as  burled  structures  and  rock  types),  but  the  interpretation  of  such  imagery  will 
be  more  difficult  and  more  subject  to  error  than  that  of  present  radar  imagery.  In  general,  the 
major  geologic  advantage  of  SLAR  is  that  It  detects  and  Identifies  surface  features  more  thor¬ 
oughly  and  therefore  aids  the  geologist  In  making  inferences  about  subsurface  features. 

Sldelooklng  radar  has  several  advan*  ges  over  multlspectral  scanning  and  aerial  photog¬ 
raphy  aside  from  the  fact  that  It  collects  information  from  a  different  part  of  the  electromag¬ 
netic  energy  spectrum.  As  an  active  sensor,  SLAR  provides  Its  own  source  of  illuminating 
energy  and  records  the  energy  reflected  back.  Therefore,  radar  Imagery  can  be  obtained  with¬ 
out  regard  to  sun  conditions.  Radar  penetrates  fog,  haze,  and  all  types  of  clouds  except  cumu- 
lonlmbjs  (thick,  dense  moisture- laden  clouds).  Although  rain  attenuates  the  radar  signal,  the 
extent  of  attenuation  depends  upon  the  radar  wavelength  ar.d  the  rate  of  rainfall;  some  radar 
wavelengths  will  not  penetrate  a  heavy  rain.  Thus,  SLAR  Is  essentially  an  all-weather,  day  or 
night  sensor  which  Is  not  affected  by  dally  variations  In  temperature  and  sun  illumination.  Ra¬ 
dar  Is  the  only  sensor  capable  of  produ'  ing  the  same  high  quality  Imagery  of  a  region  regard¬ 
less  of  the  time  of  day  or  season  of  year.  The  advantages  of  using  SLAR  to  gather  geologic 
data  in  regions  of  adverse  weather  and/or  sunlight  conditions  arc  oovlous.  Such  adverse  con¬ 
ditions  commonly  exist  in  the  polar  regions,  especially  during  winter.  Sidelooking  radar 
recently  provided  the  first  Imagery  of  the  entire  Darien  Province,  Panama,  which  Is  essentially 
cloud  covered  year  around.  In  the  past  20  years,  only  a  small  portion  of  the  province  was  ever 
photographed  because  of  the  cloud  conditions.  Berlin  and  Schaber  [36]  state  that  it  Is  readily 
becoming  apparent  that  SLAR  (opcr.itlng  at  K-band  and  longer  wavelengths)  Is  the  only  sensor 
that  can  be  used  with  a  high  degree  of  success  in  the  humid  tropics. 

Although  It  appears  possible  that  some  radar  wavelengths  penetrate  certain  types  of  vege¬ 
tation,  other  radar  wavelengths  will  not  penetrate  any  vegetative  cover.  Longer  wavelengths 
than  those  presently  being  used  very  likely  will  be  able  to  penetrate  vegetation  and  provide 
information  on  the  nature  of  the  ground  surface.  Much  more  research  Is  needed  to  determine 
the  penetrating  capabilities  of  various  radar  wavelengths,  to  ascertain  the  types  and  character¬ 
istics  of  vegetation  that  can  be  penetrated,  and  to  provide  information  on  the  extent  to  which  the 
vegetation  and  ground  surface  each  contribute  to  the  total  radar  return. 

Radar  shadows  are  created  when  the  unldlrectlonally  transmitted  energy  Is  blocked  by  a 
terrain  feature  which  prevents  the  energy  from  reaching  the  shadowed  area  behind  the  feature. 
Relative  relief  and  roughness  of  the  terrain  are  portrayed  on  SLAR  Imagery  mainly  by  the  dis¬ 
tribution  of  radar  shadows.  Rad.ar  imagery  Is  similar  to  a  shaded-rellef  topographic  map  in 
that  the  shadowing  creates  a  three-dimensional  effect.  Extensive  and/or  excessively  dark 
shadows,  particularly  in  high  relief  areas,  commonly  mask  features  .and  details  and  hinder 
geologic  interpretation.  By  careful  choice  of  the  look  direction  and  depression  angle,  r.idar 
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shadowing  can  eraphasize  any  preselected  direction  of  topography.  Radar  shadows  are  exten¬ 
sively  used  in  interpreting  structural,  geomorphlc,  and  local  lithologic  units  seen  on  the  im¬ 
agery.  Because  of  the  high  depression  angles  (45-60°),  shadowing  will  be  much  less  extensive 
on  radar  imagery  obtained  by  satellites  scheduled  to  be  launched  in  the  next  several  years. 
Although  the  reduced  amount  of  shadowing  will  somewhat  hinder  interpretation,  more  total  ter¬ 
rain  area  and  more  geologic  features  will  be  visible. 

5.2.  SIMULATION  OF  RADAR  IMAGERY  BY  MEANS  OF  LOW  SUN-ANGLE  PHOTOGRAPHY 
Lyon,  Mercado,  and  Campbell  [37]  have  shown  that  low  sun- angle  photography  (LSAP), 

taken  under  specified  conditions,  approximates  SLAR  imag''ry  made  at  comparatively  high  de¬ 
pression  angles.  They  state  that  by  using  a  low  sun- angle  (around  20-30°  above  the  horizon) 
combined  with  an  increase  in  the  y-contrast  of  the  prints  to  intensify  the  blackness  of  the  shad¬ 
ows,  one  can  produce  aerial  photography  with  some  of  the  shadowing  features  of  radar  imagery 
used  in  geologic  interpretation.  Low  sun-angle  photography  has  some  advantages  over  SLAR; 
it  is  less  expensive  over  limited  areas  than  radar  imaging  and  requires  no  electronic  gear. 
Also,  the  quality  of  the  imagery,  when  conditions  are  optimum,  appears  to  be  equal  to  that  of 
SLAR  imagery.  However,  in  addition  to  the  necessity  for  clear  weather,  LSAP  has  other  im¬ 
portant  disadvantages:  (1)  sun  angles  can  only  be  chosen  ircm  a  finite  set  of  elevation  and 
azimuth  angles;  (2)  a  particular  solar  elevation  and  azimuth  must  be  used  at  different  times 
of  the  year  for  a  specific  latitude  to  get  photographs  which  closely  simulate  radar  imagery. 
Lyon  et  al.  [37]  report  that  "one  should  plan  the  use  of  the  sun  in  azimuth,  as  well  as  in  ele¬ 
vation,  to  enhance  the  lineatlons.  If  north- south  linears  are  present,  one  should  use  the  dawn 
or  dusk  sunlight;  if  east-west,  the  winter  sun,  and  so  on.  Above  all,  one  should  plan  to  use 
several  Illumination  angles."  Normally  it  would  not  be  possible  to  wait  long  enough  to  satisfy 
the  conditions  necessary  for  LSAP  simulation  of  radar  imagery. 

It  should  be  remembered,  however,  that  low- sun  angle  photography  and  SLAR  do  not  pre¬ 
sent  similar  information  for  at  least  two  reasons.  First,  unlike  optical  wavelengths,  radar 
wavelengths  interact  with  surface  textures  with  roughness  on  the  order  of  radar  wavelengths. 
Second,  the  albedo  of  the  surface  does  not  correspond  to  the  dielectric  constant  at  microwave 
frequencies. 

5.3.  STEREOSCOPIC  SLAR  IMAGERY  AND  AREAL  COVERAGE 

It  is  presently  possible  to  generate  Ibnlted  amounts  of  stereoscopic  SLAR  imagery  by 
overlapping  adjacent  flight  lines.  Stereoscopic  SLAR  imagery  enables  one  to  interpret  the  im¬ 
agery  qualitatively  to  determine  if  the  sh.''dowing  is  caused  by  topography  or  by  the  differing 
reflective  properties  of  materials. 
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Without  stereoscopic  coverage,  techniques  which  aid  in  the  differentiation  of  areas  of  di¬ 
verse  relief  are  very  helpful.  As  indicated  earlier,  limited  shadowing  can  give  SLAR  imagery 
a  pseudT  three-dimensional  appearance.  Dellwig,  MacDonald,  and  Kirk  [38]  describe  a  simple, 
relatively  inexpensive  technique  which  imposes  a  pseudo  three-dimensional  effect  upon  mono- 
scopic  SLAR  imagery.  This  technique  is  easily  applicable  to  the  present  methods  employed 
for  monoscopic  imagery  viewing  and  interpretation  and  facilitates  the  interpretation  of  radar 
imagery  in  areas  of  low  relief,  such  as  the  Gulf  Coast.  Magnification  by  a  binocular  micro¬ 
scope  causes  the  three-dimensional  effect  to  become  even  more  pronounced,  allowing  an  inter¬ 
preter  to  identify  readily  features  that  formerly  were  on  the  threshold  of  recognition  [38].  The 
technique  is  known  only  to  exaggerate  the  relative  relief,  its  radargrammetric  value  not  having 
been  invest^ated. 

Sidelooking  radar  is  capable  of  providing  broad  areal  coverage,  on  either  a  single  film  or 
on  a  mosaic,  in  a  minimum  amount  of  fl’ght  time.  This  coverage  is  difficult  to  obtain  by  any 
other  sensor.  The  continuity  of  regional  structures  and  landforms  can  easily  be  traced  for 
long  distances  on  SLAR  single  images  and  mosaics,  and  the  regional,  commonly  deeper  extend¬ 
ing,  structural  pattern  can  be  determined.  This  identification  would  not  be  possible  with  aerial 
photography  without  extensive  use  of  mosaics,  which  could  very  possibly  hinder  the  tracing  of 
the  desired  structure  or  landform. 

5.4.  RADAR  RETURN  PARAMETERS 

Sidelooking  radar  imagery  portrays  the  reflectance  characteristics  of  the  terrain  for  cer¬ 
tain  wavelei^ths  of  the  electromagnetic  spectrum  that  are  much  longer  than  those  of  visible 
light.  Those  parameters  generally  agreed  upon  as  most  significantly  affecting  the  general  ra¬ 
dar  return  are  surface  roughness,  topography  (including  orientation  and  slope  of  surfaces), 
complex  dielectric  constants  of  soils  and  rocks,  object  geometry,  moisture  content  of  the  fea¬ 
ture  (especially  near  the  surface),  type  and  extent  of  vegetation,  moisture  content  of  low  vege¬ 
tation  (such  as  brush,  crops,  and  grasses),  and  the  incidence  angle,  wavelength,  and  polariza¬ 
tion  of  the  transmitted  radar  energy.  Rydstrom  [39]  also  considers  the  manner  of  illumination 
by  radar,  i.e.,  unidirectional  illumination,  to  be  a  fundamental  parameter.  The  relative  value 
of  the  different  parameters  commonly  is  difficult  to  assess,  even  after  a  comprehensive  study. 
Generally,  several  parameters  contribute  to  the  magnitude  of  the  radar  return,  but  one  param¬ 
eter  may  be  dominant.  The  relative  importance  of  each  parameter  depends  upon  the  feature 
and  region  being  imaged. 

A  surface  with  roughness  less  than  the  wavelength  of  the  incident  electromagnetic  energy 
appears  smooth  and  flat  to  the  energy  and  acts  as  a  specular  or  mirror-like  reflector.  The 
specular  reflection  may  be  either  a  strong  return  or  a  no- return,  depending  upon  the  angle  at 
which  the  energy  strikes  the  surface.  Rougher  surfaces  appear  multifaceted  to  the  incident 
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energy  and  act  as  diffuse  reflectors,  scattering  energy  unequally  in  all  directions.  Diffuse  re¬ 
flectors  return  only  a  fraction  of  the  energy  to  the  antenna.  Because  radar  wavelengths  are 
much  longer  than  those  of  visible  light,  many  more  surfaces  appear  smooth  to  radar  energy 
than  to  visible  light.  Thus,  radar  imagery  displays  a  greater  number  of  specular  reflections 
than  does  aerial  photography.  Depending  upon  the  particle  diameters  and  radar  wavelength, 
silt,  sand,  and  gravel  deposits  can  act  as  specular  reflectors.  Water  bodies  with  smooth  hori¬ 
zontal  surfaces  act  as  specular  reflectors  and  are  recorded  as  no- returns  on  the  imagery. 

The  complex  dielectric  constant  of  the  soil  or  rock  is  directly  determined  by  its  chemical 
and  mineralogical  composition,  moisture  content,  conductivity,  and  porosity.  What  part  each 
of  these  factors,  plus  possibly  others,  plays  in  determining  the  complex  dielectric  constant  of 
a  material  is  generally  hard  to  determine.  Gillerman  [40]  noted  a  difference  in  the  return  on 
simultaneously  obtained  like-  and  cross- polarized  K-band  radar  imagery  of  rhyolitic  to  basaltic 
Tertiary  volcanic  rocks  ’r.  the  western  United  States.  He  suggests  that  mineral  composition  an 
and  glass  content  may  play  a  part,  although  very  likely  a  subtle  one,  in  causing  the  difference 
in  return.  He  also  speculates  that  the  age  of  the  rocks,  as  it  might  be  related  to  compositional 
differences  and  degree  of  weathering,  may  be  a  factor  contributing  to  the  difference  in  the 
radar  return  of  the  lava  domes  and  craters  of  Mono  Craters,  California.  The  complex  dielectric 
constant  varies  directly  with  the  moisture  content  of  the  soil  or  rock,  and  high  moisture  con¬ 
tent,  especially  near  the  surface  of  the  feature,  causes  a  bright  or  very  light  toned  return  on 
the  radar  imagery.  Examination  of  imagery  during  the  course  of  this  contract  suggests  that  a 
high  moisture  content  of  soil  and  rock  may  so  dominate  the  other  radar- return  parameters  as 
to  be  essentially  the  only  parameter  influencing  the  return.  Accordingly,  radar  imagery  of  a 
region  generally  should  be  obtained  during  a  dry  period.  The  texture  of  the  rock,  including 
crystal  orientation  and  possibly  crystal  size,  influences  the  rock's  conductivity.  Cosgriff, 

Peak,  and  Taylor  [41]  discuss  the  variations  in  radar  return  from  materials  of  many  different 
compositions. 

The  angular  relationship  between  the  surfaces  of  a  geologic  feature,  i.e.,  the  object  geom- 

greatly  influences  the  strength  of  the  radar  return.  Hadar  energy  may  or  mav  not  be 
specularly  reflected  back  to  the  antenna  from  a  right-angled  dihedral  corner  reflector.  A 
trihedral  corner  reflector  in  which  the  three  surfaces  are  mutually  perpendicular  will  return 
a  strong  specular  reflection  to  the  antenna  regardless  of  the  look  direction.  Rydstrom  [39] 
ably  Ulustrates  the  relationship  between  dihedral  corner  reflectors  and  specular  reflection. 

Vegetation  causes  diffuse  radar  reflections.  Because  the  radar  eneriry  is  primarily  re¬ 
turned  from  the  first  reflecting  surface  that  is  encountered,  the  return  may  be  mainly  from, 
and  thus  more  indicative  of,  the  vegetative  cover  than  the  underlying  soil  or  rock.  Thus,  the 
type  and  extent  of  vegetation  are  very  important  radar  return  parameters.  Radar  Lmagery  of 
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areas  with  a  heavy  vegetative  cover  will  appear  much  different  than  that  of  areas  without  veg¬ 
etation  or  with  onlj'  sparse  vegetation.  The  differences  in  the  radar  images  of  various  types 
of  vegetation  may  indirectly  indicate  differences  in  the  underlying  soil  or  rock,  such  as  the 
thickness  of  the  soil  or  the  mineral  content  of  the  soil.  The  occurrence  and  distribution  of 
substantial  forested  areas  may  indirectly  indicate  the  topography.  The  radar  return  varies 
directly  with  the  moisture  content  of  low  vegetation;  h^h  moisture  content  causes  high-energy, 
light- toned  returns. 

The  incidence  angle,  wavelength,  aud  polarization  of  the  transmitted  radar  beam  are  the 
only  radar- return  parameters  which  can  be  chosen  and  directly  controlled.  The  incidence 
angle  on  a  particular  surface(s)  and  object  geometry  control  how  much  energy,  if  any,  is  spec¬ 
ularly  reflected  from  the  surface(s)  back  to  the  antenna.  Radar  shadowing,  radar  layover,  and 
radar  foreshortening  of  a  particular  terrain  are  functions  of  the  general  incidence  angle  of  the 
terrain,  v/hich  is  the  complement  of  the  depression  angle  of  the  radar  beam  (Refs.  42  and  43). 

The  same  feature  commonly  looks  different  on  imagery  of  different  wavelengths,  and  the 
interpretation  of  the  imagery  of  the  two  wavelengths  commonly  yields  different  information 
about  the  feature.  Thus,  it  is  desirable  to  image  the  same  area  with  different  radar  wave¬ 
lengths,  simultaneously  if  possible.  Not  only  are  the  penetration  characteristics  of  the  various 
wavelengths  different  (although  none  of  the  present  radar  wavelengths  can  penetrate  to  more 
than  a  very  shallow  depth),  but  the  speculai'-diffuse  reflection  boundary  in  terms  of  particle 
size  is  different.  Even  though  more  areas  are  being  imaged  by  two  different  radar  wavelengths, 
it  has  not  been  possible  previously  to  do  so  simultaneously.  The  inability  to  image  simultan¬ 
eously  has  obvious  disadvantages,  since  return  parameters  such  as  moisture  content  of  the 
soil  or  rock  and  vegetation  characteristics  (type,  extent,  moisture  content)  can  easily  change 
between  flights.  Willow  Run  Lalxjratories,  by  means  of  a  new  multiplexed  radar  system,  will 
soon  have  the  capability  of  simultaneously  imaging  an  area  with  two  different  radar  wavelengths. 

Afeature  commonly  appears  different  on  the  simultaneously  produced  like-  and  cr  o&s -polar¬ 
ized  radar  imagery  of  a  particular  wavelength,  and  the  difference  in  appearance  provides  inlcrmation 
about  the  feature.  Four  differently  polarized  images  generally  are  obtained  for  one  radar 
wavelength.  Radar  energy  is  transmitted  linearly  polarized  (either  vertical  or  horizontal)  and 
is  returned  to  the  antenna  either  linearly  polarized  to  produce  like-polarized  (HH  and  VV)  ra¬ 
dar  imagery  or  orthogonally  polarized  to  produce  cross-polarized  (HV  and  VH)  radar  imagery. 

Certain  return  parameters  are  more  important  than  others  in  determining  the  appearance 
of  geologic  features  on  radar  imagery.  The  most  important  parameters  generally  are  topog¬ 
raphy,  surface  rougimess,  object  geometry,  complex  dielectric  constants,  and  the  incidence 
angle,  wavelength,  and  polarization  of  the  transmitted  radar  energy.  Rydstrom  [39]  would  add 
the  unidirectional  illumination  of  the  terrain  by  radar  as  a  fundamental  return  parameter;  it 
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may  or  may  not  be,  but  certainly  it  is  a  large  factor  in  how  the  geologic  feature  lodes  on  the 
radar  imagery. 

At  the  present  time,  the  other  parameters  are  generally  considered  to  be  much  more  im¬ 
portant  than  the  complex  dielectric  constant  (composition)  of  the  soil  or  rock  in  determining 
the  radar  return  of  a  geologic  feature;  indeed  the  contribution  to  the  return  by  the  other  param¬ 
eters  masks  the  lesser  contribution  made  by  the  compKJsition  of  the  material.  However,  certain 
compositional  characteristics  of  rocks  may  be  inferred  from  their  physical  configurations, 
which  affect  the  radar  Linage  more  directly  [39].  Generally,  it  is  only  possible  to  differentiate 
between  crystalline  and  sedimentary  rocks  and  then  only  on  the  basis  of  the  inherent  character¬ 
istics  of  the  rock  classes  (bedding,  jointing,  object  geometry)  andjmt  on  the  basis  of  mineralogy 
or  chemistry. 

The  unidirectional  illumination  of  the  terrain  by  radar  energy  creates  radar  shadows 
which  give  the  imagery  a  pseudo  three-dimensional  appearance  and  indicate  the  topography 
and  relative  relief.  Radar  imagery  generally  shows  smaller  variations  in  surface  relief  than 
does  aerial  photography.  The  topography  and  relief,  discussed  later,  are  important  parameters 
used  in  identifying  structural  and  geomorphic  features,  especially  linear  ones. 

Look  direction  is  an  important  factor  in  whether  structural  lineaments  or  linear  geomorphic 
features  are  readily  recognized.  A  study  by  MacDonald,  Kirk,  Dellwig,  and  Lewis  [44]  showed 
that  lineaments  normal  to  or  at  a  high  angle  to  the  lock  direction  are  much  more  easily  seen 
than  those  that  are  not.  Lineaments  parallel  to  the  look  direction  commonly  are  not  seen  at 
all.  Obviously,  the  geologic  feature  may  be  enhanced  or  not  seen  depending  upon  its  orientation 
to  the  look  direction.  Lineaments  not  seen  or  only  faintly  seen  on  the  imagery  obtained  with 
one  look  direction  may  be  clearly  visible  on  the  imagery  obtained  from  another  look  direction. 
MacDonald  et  al.  [44]  state  that  the  maximum  geologic  information  from  radar  imagery  of 
an  unmapped  or  poorly  mapped  region  can  be  obtained  by  imaging  the  region  from  four  orthog¬ 
onal  look  directions,  but  that  imagery  taken  using  only  one  look  direction  can  still  contain 
significant  information. 

A  surface  composed  of  particles  with  diameters  less  than  the  wavelength  of  the  incident 
radar  energy  appears  smooth  and  flat  and  generates  specular  reflections.  Thus,  silt,  sand, 
gravel,  and  glacial  deposits  commonly  can  be  differentiated  as  to  the  type  of  deposit,  but  not 
as  to  composition,  by  determining  whether  they  act  as  a  specular  or  diffuse  reflector  for  a 
particular  wavelength.  This  differentiation  of  grain  size  is  best  done  by  imaging  the  deposit 
with  two  or  more  radar  wavelengths,  such  as  will  be  done  with  the  new  multiplexed  radar  sys¬ 
tem.  Surfaces  that  appear  rough  to  one  wavelength  may  appear  smooth  to  another.  Smooth, 
essentially  horizontal,  specular  surfaces  (such  as  playas)  are  recorded  as  no- returns  and 
generally  are  easily  recognized  in  arid  environments.  However,  areas  within  a  playa  of  high 

74 


WILLOW  RUN  LABORATORIES 


near- surface  moisture  content  and/or  vegetation  content  may  be  imaged  as  local  bright  spots. 
Using  shorter  wavelength  radars,  alluvial  fans  can  be  readily  distinguished  from  playas  by 
their  diffuse  reflections  as  well  as  by  their  characteristic  fan  shape.  Coarser  debris  such  as 
boulders  in  a  dry  stream  bed|  distributary  channels  of  coarse  material,  and  talus  slopes  gen¬ 
erally  will  show  as  high  returns  on  shorter  wavelength  radar  imagery,  a  result  which  can  be 
misleading.  Lava  flows  have  rough  surfaces  and  thus  generally  show  a  moderate  to  high  re¬ 
turn  on  radar  imagery. 

Topography  and  surface  roughness  are  particularly  important  parameters  in  determining 
the  radar-return  signal  that  is  especially  applicable  to  the  recognition  of  linear  features  and 
geologic  boundaries  (GUlerman,  Ref.  45).  There  commonly  is  a  correspondence  between  the 
topography  and  the  underlying  structure  or  geomorphic  feature,  and  this  relationship  is  par¬ 
ticularly  well  displayed  on  SLAR  imagery.  The  topography  can  reflect  either  different  resis¬ 
tances  of  rocks  to  erosion  or  structural  or  geomorphic  patterns.  MacDonald's  study  [42]  of  the 
radar  imagery  of  Darien  Province,  Panama,  showed  that  in  areas  of  substantial  forest  cover, 
where  the  radar  return  is  primarily  from  the  foliage,  the  imagery  indirectly  depicts  topography. 

Faults,  joints,  and  contacts  between  different  lithologic  units  commonly  have  relatively 
shallow  depressions  along  them  which  are  enhanced  by  the  shadowing  effect  of  SLAR.  Gen¬ 
erally,  many  more  of  these  structural  features  can  be  seen  on  SLAR  than  on  aerial  photography 
or  as  proven  by  many  studies,  such  as  multispectral  scanning,  Gillerman  [45].  Topographically 
expressed  folds  are  clearly  displayed  and  traced  on  SLAR  imagery.  Indeed,  the  regional  fold 
and  fracture  patterns  commonly  are  more  easily  delineated  on  SLAR  imagery,  with  its  synoptic 
view,  than  on  aerial  photographs. 

It  must  be  kept  in  mind  that  shapes  are  distorted  by  the  near- range  compression  of  slant 
range  radar  imagery.  The  elongated  distortion  of  shapes  parallel  to  the  ground  track  is  a 
problem  when  delineating  landforms  which  might  have  structural  significance  (Lewis  and 
MacDonald,  Ref.  43).  Wing  [46]  and  Wing,  Overbey,  and  Dellwig  [47]  point  out  that  is  is  im¬ 
portant  to  distinguish  on  slant  range  imagery,  with  its  nonuniform  distortion,  between  linea¬ 
ments  that  are  truly  straight  and  those  that  only  appear  to  be.  Straight  parallel  lines  on  the 
ground  do  not  appear  such  on  the  imagery,  unless  they  are  either  90°  or  180°  to  the  flight  line. 
Slightly  curved  lineaments,  approximately  parallel  to  the  flight  line,  may  appear  straighter 
in  the  near  range,  and  straight  lineaments  can  appear  curved  if  they  are  long  enough  and  not 
oriented  90°  to  the  flight  line.  Much  of  the  problem  of  lineament  distortion  is  compensated 
for  by  plotting  the  lineaments  on  topographic  maps,  provided  large  scale  coverage  is  avail¬ 
able  [46,  47]. 

The  object  geometiy  of  rocks  plays  a  large  part  in  their  appearance  on  radar  imagery. 

Rocks  with  surfaces  that  form  dihedral  or  trihedral  corner  reflectors  commonly  return  much 
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more  energy  to  the  pjitenna,  and  thus  they  appear  much  lighter  toned  on  the  imagery  than  they 
would  if  the  naturrJ  reflectors  were  not  present.  Trihedral  corner  reflectors  will  return  a  high 
amount  of  energy  regardless  of  the  look  direction,  but  the  highest  return  from  dihedral  corner 
reflectors  is  from  those  oriented  normal  to  the  look  direction.  Sedimentary  rocks,  such  as 
limestone  or  dolomite,  in  which  there  are  a  large  number  of  joint  planes  perpendicular  to  each 
other  and  to  the  bedding  planes  commonly  contain  many  dihedral  and  trihedral  corner  reflectors 
High  returns  can  be  expected  from  the  cliff  faces  of  such  rocks  and  from  valleys  cut  into  them. 
The  irregular  surfaces  of  lava  flows  commonly  form  corner  reflectors  so  that  the  return  from 
lava  flovs  is  both  specular  and  diffuse.  Plutonic  igneous  rocks  should  contain  few,  if  any,  cor¬ 
ner  reflector  surfaces,  while  the  number  of  corner  reflectors  to  be  found  on  metamorphic  rocks 
is  determined  by  the  amount  of  foliation  on  the  rock,  its  resistance  to  erosion,  and  the  degree 
of  jointing. 

5.5.  USE  OF  WIDE  DYNAMIC  RANGE  OF  RADAR 

Both  synthetic  and  real  aperture  radar  systems  are  capable  of  a  wide  dynamic  range  of  up 
to  approximately  60  db.  However,  images  used  for  interpretation  are  most  often  recorded  on 
photographic  paper  with  a  dynamic  range  of  perhaps  15-20  db.  Geologic  diagnostics  available 
from  subtle  tonal  gradations  are  lost  when  such  processing  methods  are  used.  Usable  data 
are  degraded  and  useful  diagnostic  information  destroyed.  The  degradation  of  the  imagery  is 
approximately  1/2^°  (1/1024)  in  energy  or  1/2®  (1/32)  in  amplitude  of  the  avaUable  tonal  grada¬ 
tion  (conservatively  based  on  50-db  recording  and  20-db  photographic  printing), 

Although  the  signal  film  of  the  synthetic  aperture  radar  will  only  record  a  dynamic  range  of 
approximately  20  db,  the  signal  is  spread  over  a  relatively  large  area  and  is  then  compressed 
to  a  small  area.  The  ratio  of  areas  is  the  factor  of  improvement  in  dynamic  range.  Recently 
N.  Massey  of  Willow  Run  Laboratories  developed  a  means  of  making  from  signal  film  a  holo¬ 
gram  which  can  be  viewed,  or  scanned,  in  a  simple  optical  viewer.  The  subtle  gradations  of 
tone  resulting  from  the  larger  dynamic  range  are  thus  accessible.  Also,  Willow  Run  Labora¬ 
tories  is  presently  installing  an  image  dissector  with  analog  to  digital  conversion  to  obtain  a 
wide  dynamic  range  and  to  enable  sophisticated  digital  processing.  During  this  contract,  how¬ 
ever,  we  did  not  have  holograms  of  suitable  geologic  features,  nor  was  the  image  dissector 
available  to  determine  the  usefulness  of  a  wide  dynamic  range  to  the  selection  of  tunnel  sites. 

A  wide  dynamic  range,  properly  used,  might  allow  for  the  successful  identification  of 
fault  zones,  an  important  feature  in  the  selection  of  tunnel  sites.  Without  such  a  wide  dynamic 
range,  our  ability  to  identify  such  zones  has  been  limited.  For  example,  in  the  Arbuckle 
Mountains  of  C&lahoma,  a  fault  zone  had  been  identified  by  means  of  predawn  far-inirared 
imagery  [5]  which  revealed  a  drop  in  temperature  caused  by  evaporation  of  water,  with  a 
slightly  greater  amount  of  moisture  in  the  fault  zone  than  in  the  environs. 
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We  examined  classified  radar  imagery  of  the  area  in  an  attempt  to  identify  this  fault  zone, 
assuming  that  radar's  sensitivity  to  moisture  would  help  in  the  delineation  of  the  feature.  Dif¬ 
ficulties  in  optical  processing  produced  unsubdued  speckle  patterns,  and  the  feature  could  not 
be  identified.  Later  a  technique  termed  noncoherent  integration  was  used  to  process  the  im¬ 
agery  to  bring  out  ground  painting.  In  the  processed  imagery,  the  fault  zone  appeared  slightly 
darker  in  tons,  and  the  tonal  difference  was  similar  to  that  found  in  some  published  works  il¬ 
lustrating  the  use  of  radar  imagery  in  identification  of  fault  zones.  This  tonal  difference  might 
iTave  been  caused  by  a  slight  difference  in  vegetation,  however,  and  microdensitometer  traces 
across  the  film  were  inconclusive.  Had  a  larger  dynamic  range  been  available,  the  subtle 
gradations  in  tone  might  have  allowed  precise  identification  of  this  fault  zone. 

During  the  period  of  the  ensuing  contract,  we  will  have  available  holograms  of  geologic 
features  and  an  image  dissector  with  analog  to  digital  conversion,  both  of  which  will  produce 
a  wide  dynamic  range.  With  the  image  dissector,  level  slicing  can  be  performed  on  the  digital 
output  of  the  optical  processor.  Either  the  output  of  the  optical  processor  or  the  holographic 
viewer  image  can  be  used  for  level  slicing.  With  the  larger  dynamic  range  from  which  to 
choose,  level  slicing  becomes  a  powerful  diagnostic  tooi,  which  may  well  be  used  effectively 
in  the  identification  of  geologic  features. 

5.6.  OPTICAL  PROCESSING  FOR  ENHANCEMENT  OF  GEOLOGIC  FEATURES* 

Coherent  radar  imagery  has  shown  considerable  usefulness  for  all-weather  mapping  of 
terrain.  The  complex  scattering  mechanism  in  radar  imaging  is  strongly  affected  by  geomor- 
phic  factors  and  tends  to  enhance  subtle  topography  and  vegetation  [39].  The  ability  of  a  radar 
image  to  reveal  small  tectonic  features  and  lineaments  has  been  recognized  by  others;  Rowan 
and  Cannon  [8]  have  utilized  radar  imagery  to  locate  both  previou.sly  mapped  and  newly  dis¬ 
covered  faults  in  the  Mill  Creek,  CMahoma  area. 

Although  radar  imagery  effectively  displays  a  large  number  of  geomorphic  features,  it  still 
absorbs  a  considerable  amount  of  a  geologist's  activity  to  locate  and  map  tectonic  features. 
Particularly  in  view  of  the  image  coverage  afforded  by  a  radar  sensor,  whereby  geologic  map¬ 
ping  can  be  extended  over  much  of  the  earth's  surface,  the  present  capacity  to  collect  radar 
imagery  is  far  beyond  that  of  geologists  to  examine,  analyze,  and  interpret  it  adequately. 

This  problem  could  be  alleviated  by  processing  the  imagery  to  reduce  the  amount  to  be 
examined,  allowing  the  geologist  to  use  his  talents  for  interpretation  of  the  structural  signatures. 
Since  the  images  are  two-dimensional,  it  is  natural  to  consider  a  form  of  optical  processing 
which  can  operate  upon  the  imagery  without  transformation. 

*Section  5.6,  authored  by  Ralph  Mitchel,  was  originally  published  in  the  semiannual  re¬ 
port  on  this  contract.  Ref.  [16]. 
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Coherent  optical  systems,  in  particular,  are  suitable  for  interpreting  radar  imagery. 
Along  with  the  capability  of  analyzing  each  point  of  the  image  simultaneously,  coherent  optical 
image- processing  systems  possess  a  physically  accessible  Fourier  transform  plane,  allowing 
the  spectrum  of  an  image  to  be  modified.  Thus,  linear  filtering  operations  can  be  performed 
upon  the  imagery  very  rapidly.  Examples  of  the  filtering  operations  commonly  performed  by 
a  coherent  optical- processing  system  are  directional  filtering,  matched  filtering,  and  image 
enhancement. 


A  series  of  experiments  was  performed  on  some  hig.i- resolution,  synthetic- apertime  radar 
imagery  collected  by  The  University  of  Michigan  to  evaluave  the  usefulness  of  coherent  optical 
processing  as  an  aid  to  the  geolc^ist.  The  area  mapped  was  near  Mill  Creek,  Oklahoma,  an 
area  underlain  by  gently  folded,  severely  faulted  pre- Pennsylvanian  rocks,  parts  of  which  are 
covered  by  Pennsylvanian  sandstone  and  conglomerates  of  tl.e  lower  part  of  the  Pontotoc  Group. 
The  major  faults  and  folds  in  these  pre- Pennsylvanian  rocks  run  northwest- southeast.  In  many 
regions,  the  conglomerates  lack  the  extensive  structures  which  might  revea!  offsetting.  Rubble 
along  a  fault  trace  appears  similar  to  the  parent  rock,  making  the  recognition  of  a  fault  diffi¬ 
cult,  even  with  meticulous  field  work. 


Optical  directional  filtering  experiments  were  attempted  on  this  imagery  with  little  si  cess. 
Narrow  angular  wedges  placed  in  the  optical  system  to  modify  the  spectrum  of  the  radar  imag¬ 
ery  produced  poor  discrimination.  Although  the  filtered  image  had  poor  contrast,  it  did  con¬ 
tain  many  streaks  identifying  lineaments.  However,  this  type  of  filtering  did  not  effectively 
reduce  the  information  to  be  interpreted.  Similar  results  were  obtained  when  a  two-dimensional 
Gaussian  slit  was  used  to  modify  the  spectrum  of  the  imagery. 

In  an  effort  to  obtain  better  discrimination,  optical  matched  filtering  was  attempted.  The 
mathematical  model  that  we  utilized  for  the  radar  energy  reflected  from  a  fault  tr.ace  and  re¬ 
corded  on  a  radar  image  was  as  follows; 


X  -  -i  (y) 
i(x,  y)  =  Rect-^^ 


Rect 


y-yp 

2a 


where  i  =  the  radar  im^ery  reflected  from  a  fault  trace  and  recorded  on  a  radar  image 
X  =  coordinate  across  width  of  fault  trace 
y  =  coordinate  along  length  of  fault  trace 

1  =  a  constant 
o 

2  A  =-■  width  of  the  fault  trace 
2a  =  length  of  the  fault  trace 
x^  =  path  taken  by  a  fault  trace 

x^(y)  =  center  of  fault  trace  as  a  function  of  the  location  along  its  length 
y^  =  center  of  fault  length 
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2A  and  are  written  as  functions  of  the  length  y,  since  the  trace  may  not  be  of  constant  width 
or  exactly  straight. 

We  defined  probability  measures  over  x^,  A,  and  a,  along  with  the  associated  probability 

density  functions  f  (x)  and  f  (x)  over  x,  as  well  as  f  (y)  over  y.  x  (y)  and  A(y)  are  random 
o  a  o 

processes;  i.e.,  for  each  y  they  are  random  variables.  We  assumed  that  f  (x)  is  symmetric 

*0 

around  some  value  x^(y)  =  M  (the  centroid  of  a  fault  trace  in  the  x-direction);  that  f^(x)  is 
symmetric  at  about  A(y)  =  X  (average  width  of  a  fault  trace);  and  that  ^^(y)  is  symmetric  at 
about  some  mean  value  a  =  a  (average  length  of  a  fault  trace). 

The  probability  i(x,  y)  =  is  the  probabUity  (P)  that  x^  -  A  s  x  s  +  A  and  -  a  s  y  <  a, 
which  can  be  computed  to  be; 

P[i(x,  y)  =  I^]  =|[f^(-x)  -  f  ,^(x)]  ♦  F^jx)|[F_Jy)  -  Fjy)] 

=  |[F.^(x)  -  F^(x)]  ♦  f^^(x)|[F_^(y)  -  FJy)] 

where  F  =  the  probability  distribution  function. 

A  particular  case  should,  however,  be  noted;  when  the  fault  trace  is  narrow,  having  a  con¬ 
stant  width  (2A  =  2e),  where  2e  is  on  the  order  of  the  radar  image's  resolution,  we  have: 

P[i(x,  y)  =  I^]  =  f^  (x)(2e)[F_^(y)  -  F^(y)] 

That  is,  the  probability  of  the  fault  trace's  existence  at  a  particular  point  is  defined  in  terms  of 
the  probability  density  function  of  the  random  variable  defining  its  path. 

For  the  case  of  uniform  noise,  a  coiierent  optical  correlator  using  a  matched  filter  in  the 
form  of  a  complex  spatial  filter  can  display  an  output  light  amplitude  distribution  proportional 
to  the  cross  correlation  between  a  reference  signal  (stored  on  the  filter)  and  the  input  image. 
The  results  of  the  experiment  utilizing  optical  matched  filtering  were  mixed.  Some  of  the  known 
faults  were  not  detected.  However,  with  the  filter  described  above,  we  were  able  to  locate  the 
contact  separating  the  Deese  and  Dornick  Hills  formations. 

The  last  experiment  performed  was  image  enhancement,  designed  to  eliminate  from  the 
radar  image  that  information  which  is  unnecessary  to  geologic  interpretation.  The  low  spatial 
frequency  portion  of  the  image's  spectrum  was  removed  by  means  of  a  spatial  filter.  In  the 
processed  image,  structural  edges  and  lines  dividing  differently  reflecting  terrain  regions 
were  emphasized.  As  a  result,  this  experiment  allowed  geologic  interpretation  with  much  less 
time  and  effort  than  that  required  by  the  unprocessed  pictures. 
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6 

FUTURE  RESEARCH 

6.1.  RECOMMENDATIONS 

Future  research  in  the  applications  of  remote  sensing  techniques  to  tunnel- site  selection 
should  be  concentrated  in  three  categories:  (1)  field  test  cases;  (2)  continued  development  of 
remote  sensing  systems  as  defined  under  this  contract,  i.e.,  multispectral  scanning,  sidelooking 
radar,  and  aerial  photography;  and  (3)  the  addition  of  other  kinds  of  sensors,  including  physical 
field  sensors,  for  a  complete  remote  sensing  system. 

6.1.1.  FIELfy  TEST  CASES 

Although  Wf  know  that  we  can  obtain  synoptic  remote  sensing  data  which  can  complement 
information  acquired  from  other  types  of  geologic  investigations,  these  data  must  be  evaluated 
experimentally.  For  such  an  evaluation,  we  need  a  test  area  comprised  of  crystalline  and  other 
rock  outcrops  in  which  folding,  faulting,  tilting,  jointing,  and  weathering  have  occurred. 

In  the  continuation  of  the  work  described  in  this  report,  we  have  chosen  as  a  test  site  a 
region  in  the  Black  HUls,  South  Dakota,  where  the  surface  exhibits  great  geologic  and  litho¬ 
logic  diversity.  Sidelooking  radar  imagery  of  this  region  has  been  received  from  the  Strategic 
Air  Command,  and  multispectral  scanning  and  concurrent  stereographic  aerial  photographic 
flights  were  obtained  in  May  1972.  Ground-truth  studies  of  the  test  site  were  made  at  the  time 
of  the  flights.  We  will  undertake  to  use  all  of  the  remote  sensing  data  to  infer  the  geology  and 
lithology  by  attempting  to  reconstruct,  extend,  and  possibly  correct  the  present  knowledge  of 
the  geology  of  the  region,  and  we  will  then  evaluate  our  efforts. 

Remote  sensing  tests  of  present  tunnel  sites  in  which  geologic  difficulties  have  been  en¬ 
countered  should  also  be  undertaken.  If  the  precise  geologic  features  which  caused  the  diffi¬ 
culties  could  be  identified  or  detected  for  subsequent  identification,  the  value  of  remote  sensing 
techniques  in  tunnel- site  selection  would  be  confirmed.  Unfortunately,  no  tests  of  this  kind 
are  presently  funded. 

6.1.2.  CONTINUED  DEVELOPMENT  OF  MULTISPECTRAL  SCANNING  AND  RADAR 
SYS"  MS 

Within  this  contract  period,  for  the  first  time  multispectral  scanning  systems  have  been 
shown  to  be  capable  of  identifying  the  chemical  content  of  rocks.  As  described  in  Section  4, 

SiOg  and  ferric  oxides  have  been  identified  by  means  of  a  ratioing  technique  which  takes  ad¬ 
vantage  of  the  fact  that  spectral  characteristics  of  radiation  depend  upon  the  material  which 
emits  or  reflects  that  radiation.  Although  Vincent  has  demonstrated  this  capability  with  only 
two  chemical  compositions,  the  basic  technique  could  be  used  for  many  different  geologic  ma¬ 
terials.  Investigations  into  the  applicability  of  this  technique  to  the  identification  of  other 
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geologic  materials  should  be  continued.  With  the  further  development  of  such  techniques,  more 
complex  lithology  could  be  classified  and  more  effective  geologic  reconnaissance  performed. 

Although  multichannel  processing  of  data  by  means  of  statistical  pattern  recognition  is 
quite  weU  developed,  its  application  to  geologic  investigations  has  been  limited.  Wagner  shows, 
however,  the  applicability  of  the  method  to  the  identification  and  classification  of  surface  fea¬ 
tures.  Although  much  remains  to  be  done  to  correlate  this  technique  with  other  data  processing 
methods,  without  doubt  multichannel  processing  not  only  utilizes  the  spectrum  in  a  more  sophis¬ 
ticated  way  than  does  ratioing  but  also  can  be  used  to  identify  materials  relevant  to  the  selection 
oi  tunneling  sites. 

The  structural  geologist  may  be  able  to  utilize  the  technique  at  optimal  effectiveness  by 
interacting  directly  with  the  computer,  since  geologic  inferences  require  considerable  knowl¬ 
edge  and  experience.  The  scene  could  be  presented  on  the  screen,  and  the  interpreter  could 
choose  the  spectral  bands  he  wished  to  use  and  the  defined  area  he  wished  to  investigate. 

With  a  double  cursor,  he  would  define  the  area,  and  the  computer  would  immediately  present 
on  the  screen  those  areas  which  were  within  a  preset  deviation  of  the  channels  found  within 
the  defined  area.  The  geologist  could  then  draw  upon  his  experience  to  ascertain  if  the  rock  in 
the  different  areas  was  of  similar  composition  or  possibly  from  the  same  formation.  Such 
interaction  with  the  computer  might  offer  more  effective  data  processing  for  geologic  investi¬ 
gations,  and  such  processing  techniques  should  be  developed  to  be  used  with  known  spectral 
signatures. 

The  emphasis  of  future  developments  in  sidelooking  radar  is  upon  multipolarity  and  multi¬ 
frequency  radiation.  Willow  Run  Laboratories  is  presently  capable  of  collecting  both  X-band 
(3  cm)  and  L-band  (30  cm)  radar  imagery;  however,  data  are  not  presently  collected  simulta- 
leously  in  the  two  bands.  On  one  pass  imagery  is  collected  in  one  band,  and  then  the  antenna 
is  changed  and  an  attempt  made  to  collect  imagery  of  the  same  area  and  the  same  flight  di¬ 
rection  in  the  other  band.  Since  the  imagery  from  the  two  flights  can  not  be  superimposed,  we 
have  not  been  able  to  process  both  bands  automatically  by  means  of  ratioing  or  image  sub¬ 
traction. 

Willow  Run  Laboratories  will  soon,  however,  install  a  new  multiplexed  radar  system  which 
will  allow  the  simultaneous  collection  of  both  X-  and  L-band  imagery  on  one  pass  over  an  area. 
Superimposible  imagery  will  allow  automatic  processing,  and  the  bands  can  be  ratioed  to  de¬ 
lineate  surface  textures.  Simultaneously  collected  radar  data  will  extend  the  possible  geologic 
applications  discussed  in  Section  5. 
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6.1.3.  THE  ADDITION  OF  OTHER  KINDS  OF  SENSORS  FOR  A  COMPLETE  REMOTE 
SENSING  SYSTEM 

If  remote  sensing  is  used  for  geologic  reconnaissance  (and  this  report  has  shown  the  ben¬ 
efits  to  be  gained  from  its  use),  all  types  of  airborne  sensors  should  be  included,  preferably 
aboard  a  single  aircraft.  The  latter  requirement  might  not  be  feasible,  however,  in  certain 
cases,  such  as  the  inclusion  of  both  a  gamma- ray  sensor  and  sidelooking  radar;  gamma- ray 
sensing  is  most  effective  at  very  low  altitudes,  while  radar  is  used  at  a  relatively  high  altitude. 

Even  if  different  aircraft  must  be  used,  however,  the  increase  in  the  amount  of  confirming 
data  that  could  be  gathered  by  the  additional  systems  offset  the  cost.  Thus  far,  we  have  been 
concerned  with  detecting,  identifying,  and  mapping  geologic  units,  and  we  have  emphasized 
surface  structural  details,  such  as  faults,  jomts,  soil  moisture,  drainage  patterns,  and  general 
geomorphology,  since  these  data  are  necessary  to  infer  geologic  features  at  depth.  To  in¬ 
crease  our  ability  to  define  and  map  surface  features,  a  passive  microwave  sensor  system  and 
a  gamma-ray  (nuclear  radiation)  sensing  system  should  be  added. 

Sensors  of  volume  characteristics  rather  than  of  surface  characteristics  should  also  be 
utilized,  particularly  since  the  geologic  formations  at  depths  of  a  mile  or  more  are  important 
in  tunneling.  The  direct  sensing  of  the  extent  of  buried  rock  masses  could  provide  valuable 
geologic  information.  Also,  just  as  inferences  about  structures  at  depth  can  be  made  from 
surface  data,  knowledge  of  surface  formations  can  be  rendered  more  complete  by  data  on  sub¬ 
surface  features.  The  most  complete  knowledge  of  the  geology  of  a  region  is  achieved  when 
both  kinds  of  data  are  used.  Three  volume  sensors  that  should  be  considered  are  the  magne¬ 
tometer,  the  gravimeter,  and  an  electromagnetic  sounding  system. 

The  two  additional  surface  sensors  and  the  three  volumetric  sensors  are  discussed  below. 

6. 1.3.1.  Additional  Surface  Sensors 

Passive  Microwave  Sensor  Sy, stems.  A  passive  microwave  system  is  similar  to  a  multi- 
spectral  scanning  system  in  that  it  depends  upon  the  sensing  of  natural  radiation  that  is  re¬ 
flected  and  emitted  from  objects.  The  power  received  by  the  radiometer  is  considered  in 
terms  of  an  equivalent  temperature  known  as  the  "apparent  temperature"  or  "antenna  temper¬ 
ature."  Two  types  of  passive  microwave  sensor  systems  have  been  developed — the  single¬ 
value  sensing  system  and  the  scanning  system.  The  former  establishes  a  single  value  for  any 
location  under  the  flight  path,  while  in  the  latter,  successive  scans  are  recorded  normal  to  the 
flight  path  and  a  two-dimensional  image  is  obtained.  The  scanning  system  is,  of  course,  con¬ 
siderably  more  expensive  than  the  single-value  sensing  system,  but  the  two-dimensional  im¬ 
age  that  it  produces  increases  the  potential  information  to  be  gained  and  allows  correlation 
with  other  types  of  remote  sensing  imagery. 
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Edgerton  et  al.  [48]  state  that  "microwave  brightness  temperatures  of  soils  and  sediments 
are  dependent  upon  .  .  .  moisture  content,  density,  surface  roughness  (including  microrelief, 
particle  or  fragment  size  and  shape,  vegetal  cover,  etc.),  thermometric  temperature,  and 
stratigraphy  (layering,  etc.)."  These  dependencies  have  been  demonstrated  experimentally  [48], 
and  since  data  indicating  each  of  these  physical  properties  can  be  used  to  locate  geologic  map¬ 
ping  units,  the  value  of  passive  microwave  systems  in  geologic  reconnaissance  is  obvious. 

Since  it  is  known  that  fault  zones  can  be  identified  by  far-infrared  detection  of  predawn  evap¬ 
orative  cooling  [5] ,  the  capability  for  detection  of  moisture  becomes  particularly  applicable  to 
tunnel- site  selection.  Passive  microwave  sensor  systems  might,  in  fact,  be  more  effective 
than  infrared  systems  in  the  detection  of  fault  zones,  since  the  underlying  volume,  particularly 
that  close  to  the  surface,  plays  a  part  in  microwave  emission. 

Gamma-Ray  Sensors.  Airborne  gamma-ray  surveys  conducted  in  Canada  (Darnley  and 
Fleet,  Ref.  49)  and  in  Puerto  Rico  (Sohwarzer  et  al..  Ref.  50)  have  shown  the  applicability  of 
nuclelr  remote  sensing  to  geologic  investigations.  In  these  surveys,  low- altitude  flights  were 
used  to  establish  gamma- ray  signatures  for  several  substances,  such  as  thorium,  uranium, 
potassium,  bismuth,  and  thallium.  Once  the  signatures  were  known,  rock  types  could  be  identi¬ 
fied  on  the  basis  of  the  amount  and  ratios  of  these  elements.  Schwarzer  et  al.  state,  "...  this 
study  suggests  that  it  is  possible  not  only  to  discriminate  between  various  bedrock  lithologies 
but  to  identify  them,  provided  there  is  adequate  rock  exposure  or  the  overlying  soU  cover  ac¬ 
curately  reflects  the  composition  of  the  underlying  parent  material - Gamma- ray  remote 

sensing  ...  can  be  an  effective  method  of  revealing  geologic  patterns  and  serving  as  a  guide  for 

minimum  investigation  on  the  ground." 

6. 1.3. 2.  Volumetric  Sensors 

Airborne  Magnetometer.  Airborne  magnetometers  have  been  employed  operationally  since 
World  War  II  in  m.ineral  and  oil  exploration  (Dobrin,  Ref.  51),  and  any  thorough  geologic  re¬ 
connaissance  program  should  incorporate  aeromagnetic  interpretation.  The  magnetometer 
can  easily  be  incorporated  with  the  other  sensors  in  a  single  aircraft.  Magnetometer  data 
should  also  be  correlated  with  other  remote  sensing  data,  although  research  into  the  possibilr- 
ties  of  combining  and  processing  the  data  is  necessary.  For  example,  meaningful  image  ra- 
tioing  might  be  performed  with  a  magnetic  intensity  map  and  other  remote  sensing  data. 

Airborne  Gravimetry.  Nettleton  et  al.  [52]  first  investigated  airborne  gravimetry,  the 
purpose  of  which  is  to  measure  the  variation  in  the  gravitational  field  by  detecting  small  dif¬ 
ferences  in  the  weight  of  a  constant  mass  at  different  points  on  the  earth's  surface.  The  lo¬ 
cation  and  composition  of  deeply  buried  features  can  be  inferred  from  knowledge  of  the  gravi¬ 
tational  field.  Unlike  the  airborne  magnetometer,  the  airborne  gravimeter  has  not  yet  been 
incorporated  into  operational  use  in  mineral  and  oU  exploration.  According  to  Garland  [53], 


83 


WILLOW  RUN  LABORATORIES 


"It  is  now  generally  agreed  that  the  greatest  limitation  to  gravity  measurements  from  moving 
vehicles  is  the  lack  of  navigation  systems  which  could  be  used  anywhere  on  earth  to  give  a 
continuous  record  of  the  vehicle's  velocity." 

Research  on  airborne  gravity  surveys  is  still  continuing,  however.  Hall  and  Townsend  [54] 
produced  Bouguer  anomaly  maps  with  data  gathered  from  a  helicopter  survey,  and  contour  lines 
differing  in  10  milligals  were  plotted  as  close  together  as  2  km.  In  Szabo  and  Anthony's  dis¬ 
cussion  [55]  of  aerial  gravity  measurements  (in  which  Doppler  radar,  astrotracker,  hypsometer, 
terrain  clearance  radar,  and  aerial  camera  data  were  used),  a  typical  result  was  a  total  5°  x  5° 
error  of  11.8  milligals.  In  the  continuing  research  on  airborne  gravimetry,  one  of  the  prob¬ 
lems  that  must  be  solved  is  the  determination  of  the  relative  accuracy  required  in  tunnel-site 
selection  in  comparison  to  that  required  for  mineral  exploration. 

Airborne  Electromagnetic  Sounding.  Airborne  electromagnetic  sounding  (AEM)  systems 
are  essentially  very  low-frequency,  intensity- measuring  nonranging  radars  (a  few  hertz  to  a 
few  thousand  hertz),  in  which  the  transmitter  and  receiver  are  separated,  usually  by  more  than 
60  m.  The  separation  is  achieved  either  by  towing  a  glider  or  by  using  two  aircraft.  According 
to  Dobrin  [51],  the  AEM  technique  can  be  used  to  a  maximum  depth  of  approximately  1500  ft. 

The  depth  of  penetration  is  an  inverse  function  of  frequency,  where  the  optimum  frequency,  f, 
for  greatest  penetration  [51]  is 


where  h  is  depth  of  penetration  in  meters  and  p  is  resistivity  in  ohm  •  cm. 

In  a  tutorial  article,  Ward  [56]  states,  "The  purpose  of  all  of  the  totally  airborne  [AEM] 
systems  is  to  locate,  rapidly,  good  conductors  of  electricity  within  the  top  several  hundred 
feet  of  the  earth's  surface.  Such  good  conductors  might  include  graphitic  schists,  carbonaceous 
sediments,  fault  and  shear  zones,  swamps,  etc.,  as  well  as  the  primary  target,  massive  sulfide 
mineralization."  Although  not  basically  used  for  maR)ii^,  airborne  electromagnetic  sounding 
is  particularly  suited  for  locating  fault  and  shear  zones  (probably  because  of  the  increased 
moisture  content  in  these  zones),  and  the  ability  to  identify  such  zones  is  an  obvious  benefit 
in  tunnel- site  selection. 

6.2.  RATIONALE  OF  APPROACH  TO  FUTURE  RESEARCH 

Geologic  reconnaissance  for  tunnel- site  selection  has  much  in  common  with  mineral  and 
oil  exploration  and  general  geologic  mapping  in  that  it  seeks  to  infer  features  that  cannot  be 
directly  seen.  For  such  a  purpose,  every  possible  confirming  method  should  be  used,  particu¬ 
larly  in  investigations  involving  tunneling,  where  the  cost  of  construction  is  so  much  greater 
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than  the  cost  of  a  complete  geologic  and  geophysical  survey  and  where  one  geologic  oversight 
can  cost  many  times  more  than  the  total  survey. 

In  future  research  into  geologic  applications  of  remote  sensing  data,  test  sites  must  be 
utilized  to  evaluate  present  techniques,  to  provide  a  framework  for  improvement  of  these  tech¬ 
niques,  and  to  develop  new  methods  of  acquiring  useful  data.  The  progress  that  we  have  made 
in  geologic  applications  of  image  ratioing,  statistical  pattern  recognition,  and  sidelooking  radar 
should  be  continued.  All  available  remote  sensing  systems  should  be  incorporated;  the  corre¬ 
lation  of  data  from  many  different  types  of  sensors  enables  us  to  make  more  reliable  geologic 
inferences  than  we  could  make  on  the  basis  of  any  one  type  of  remote  sensing  data  or  on  a 
limited  number  of  types. 

Because  of  the  increasing  availability  of  remote  sensing  data  (such  as  the  Strategic  Air 
Command's  sidelooking  radar  imagery  and  the  data  resulting  from  the  ERTS  projectl.  recon¬ 
naissance  geologists  will  in  the  near  future  be  required  to  use  such  data  effectively.  Research 
should  currently  be  directed  toward  the  determination  and  development  of  the  most  effective 
techniques  of  utilizing  these  data  in  geologic  investigations. 
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